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ABSTRACT 


23-6 metres of K~-meson track in nuclear emulsion have yielded unbiased 
samples of 815 absorptions at rest and 106 in flight. The corresponding 
imteraction mean free path is 21-5 nea em for K~-mesons of energy below 
90 Mey. Observations are reported on the disintegration stars themselves 
and on the charged hyperons produced. It is in particular noted that the 
angular correlation between hyperon and pion in two-pronged stars is 
closely that expected for a system without any nuclear interaction in the 
final state. It is shown that this and other aspects of the situation are most 
simply explained if the K~-meson is absorbed in the tenuous outer surface of 
the nucleus from a high Bohr orbit according to the model of Jones, and 
is not absorbed uniformly throughout the nuclear volume as is usually 
assumed. This surface absorption model allows of a strong hyperon— 
nucleon interaction such as appears to be indicated by a variety of other 
evidence and which would be inconsistent with the volume absorption 
model. 


§ 1. INTRODUCTION 


INVESTIGATIONS so far indicate that the K~-meson interacts with nuclear 
matter by two distinct processes : absorption by one nucleon with the 
production of a hyperon and a z-meson; and secondly, absorption by two 
nucleons which have high relative momentum owing to correlations in the 
nuclear ground state wave function. Such interactions provide a source 
of information about the fundamental K~-meson interaction, and may 
also be used to study the properties of the capturing nucleus. Hornbostel 
and Salant (1956) have outlined the characteristics to be expected of 
K--meson absorption, and they and many others (Alles e al. 1957a, 
Bacchella et al. 1958, Barkas et al. 1957, Fry et al. 1957, Gilbert et al. 1957, 
Haskin et al. 1956, Hornbostel and Zorn 1958, and Lohrmann ef al. 1958) 
have investigated the absorption of K~-mesons in nuclear emulsion. A 
body of work as yet unpublished has been presented at recent conferences 
(Dilworth 1957, Alles ef al. 1957b, Bhowmik et al. 1957, Hoang et al. 1957, 
and Webb et al. 1957). 
Unfortunately the absorption reaction appears to be so complicated 
that only qualitative deductions can be made, and these must be based 


y+ Communicated by the Authors. This work was begun at Cambridge and 
finished at Oxford. 
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on a particular model. It has, so far, usually been assumed that a slow 
K--meson is initially captured by the electrostatic field of the nucleus and 
that it then cascades by electromagnetic transitions through a series of 
Bohr orbits until finally absorption by the nucleus occurs. One of the 
present authors (Jones 1958) has investigated such a mechanism and has 
concluded that the absorption probability, deduced from the observed 
K--proton scattering and absorption cross sections, is so large that most 
of the K--mesons will be removed from rather high circular Bohr orbits, 
the absorption taking place in the surface of the nucleus in a region of 
low nucleon density. However, evidence has been presented recently 
(Gilbert and White 1958), which seems to favour an absorption 
probability more uniform throughout the nucleus. It is important to be 
able to choose the better of the two models, viz. surface or volume 
absorption. 

Owing to the scarcity of statistically significant data on K~-meson 
interactions, we felt it advisable to make a completely systematic study 
of K~-meson capture stars in nuclear emulsion in order to obtain the 
widest variety of information. By finding K~-meson tracks where they 
enter the emulsion stack and following them until they are ultimately 
absorbed we obtain an unbiased selection of stars, as well as a value for 
the mean free path of K~-mesons in emulsion. These stars are an 
abundant source of hyperons some of whose characteristics can be studied. 
The importance of the branching ratios and angular correlations of the 
decay products is already well known. 

In §2, the experimental procedures used to find and examine the 
K~-meson interactions are described. The determination of the mean 
free path for K~-mesons in nuclear emulsion, the characteristics of K- 
capture stars and observations on the interaction and decay of the 
x -hyperons produced are all described in §3. The results obtained are 
discussed in §4. In particular, it is shown that the observed © and 7 
production angular correlation in the laboratory system of coordinates for 
2-prong stars having a &+-hyperon, a 7+-meson and no visible nuclear 
excitation, provides strong evidence for the surface absorption model. 
The consistency of the model, and the form of the effective nucleon 
momentum distribution it predicts for K~-meson absorption by a heavy 
nucleus receive some discussion. The relative production of U--hyperons 


in the alsorption of K~-mesons by bound neutrons and protons is 
determined. 


§ 2. EXPERIMENTAL DETAILS 


A stack of 125 sheets of 400 micron thick Ilford G.5 nuclear photo- 
graphic emulsion, 15 x 15 cm, was exposed to the September 1956 K~ beam 
of the Bevatron at the Berkeley Radiation Laboratory of the University 
of California. This beam was almost completely separated from 7—-meson 
contamination by an energy degrader which reduced the K- momentum. 
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from 430 to 300Mev/c and left only «--mesons of the latter momentum, 
after a second momentum analysis. The integrated K~ and p- meson 
fluxes were 70 and 2 x 10° per cm? respectively, 

In order to make a completely unbiased selection of K--meson capture 
stars and interactions in flight, and to obtain a value for the mean free 
path in nuclear emulsion, tracks having a normalized blob density in the 
range 2-1 to 2-8 were picked up near the beam entrance edge of the stack 
and followed in until some interaction occurred. All K~-mesons coming 
to rest, including zero prong stars, were then identified by a comparison 
of range with initial grain density. Interactions in flight were identified 
by multiple scattering versus grain density measurements close to the 
star. In this way the energy of the K~-meson at the point of interaction 
was obtained. 

The stars were analysed by following to rest all prongs with normalized 
blob density > 2-0. This limit corresponds to a proton energy of 175 Mev 
and to a = hyperon energy of about 220Mev. The minimum plateau blob 
density in the stack was in the range 14-16 blobs per 100 microns only, 
and so it was difficult and in many cases impossible to follow minimum 
ionizing tracks ; further, such tracks are very likely to leave the stack. 
Consequently no systematic attempt was made to obtain a 7-meson 
energy spectrum or charge ratio. 

Each K~--meson ending was examined carefully for minimum ionizing 
tracks, as were the endings of those prongs followed. With such a low 
minimum blob density, it was impossible to estimate the scanning loss for 
fast z-mesons. The frequency of 7-meson emission from our K~-meson 
_ stars is close to the generally accepted value but the number of single prong 
stars in flight may be rather low. We have avoided estimating any 
quantities which require a knowledge of the efficiency for detecting 
minimum ionizing tracks. 

Systematic identification of zero prong & -hyperon capture stars by 
Auger electron emission was not possible because of a heavy slow electron 
contamination of the stack. Decays of &+-hyperons in flight to protons 
were ascertained to be kinematically consistent with the momentum of the 
decaying particle as determined by grain density measurement. 


§ 3. RESULTS 
3.1. The Mean Free Path and Related Matters 

The interaction mean free path of K~-mesons in nuclear emulsion, 
found from events occurring along 23-56 metres of K~-meson track 
scanned, is shown in table 1. The 110 interactions are divided into three 
groups according to the K~ energy at interaction. For this determination, 
all possible decays in flight (stars with one light prong and no sign of 
nuclear excitation), and interactions with free protons were excluded. 
Inelastic scatters on nuclei, however, were included in the reaction cross 
section. 

4M2 
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Table 1. The Mean Free Path of K~-mesons in Nuclear Emulsion 


is Shown as a Function of Energy 


Be 
| 


Energy interval 5-30 30-60 60-90 Mev 
Track length 430 1070 856 cm 
No. of absorptions 18 46 42 
No. of inelastic 

scatters 0 1 3 
Mean free paths Does 22:82 3° 19-0737 cm 


Overall mean free path 21-5773 cm 


In rejecting 11 interactions in flight with only one light prong, we have 
probably rejected some true interactions. Since 44% of all K~ capture 
stars have one light prong only, the number of true interactions rejected 
could well be 5 in the present case. However, the length of track scanned 
corresponds to a proper time of 1-9 x 10~’sec, and so we would expect to 
see about 14 decays in flight if we accept the lifetime of the K~-meson to 
be 1:3x10-8sec as reported at the Venice—Padua Conference. These 
numbers are not inconsistent but may suggest some scanning bias against 
finding the very fastest particles. This has little effect on the evaluation 
of the mean free path, and the resulting figure of 21-5773 cm is in good 
agreement with 23-4+0-7cm, a figure representing the weighted average 
of the results referred to in § 17. 

Only one elastic scatter, and no absorptions of fast K~-mesons, were 
seen on free protons. Such a result is not inconsistent with an average 
cross section for elastic scattering on protons of 40 mb (Alles et al. 1957a, 
Bhowmik et al. 1957, and Glasser et al. 1958), which would lead us to 
expect 3 such scatters. 

Data on the inelastic scattering of the K~-mesons are presented -in 
table 2. 


Table 2. Inelastic Scattering of K~-Mesons on Emulsion Nuclei. 
The Data are for the Laboratory System of Coordinates 


Incident K~-meson energy | Angle of scatter | Energy loss 


87 Mev 47 degrees 48 Mev 
84. iyi 24 
80 104 61 
ay) 97 8 


+ This mean free path is considerably shorter than the geometrical mean 
free path of about 30cm. This is possible because the K--meson can be 
absorbed on a single nucleon and so can be caught in the tenuous nucleon 
atmosphere outside the main body of the nucleus (cf. p). Contrast the situation 
for the 7-meson which must find two nucleons to be absorbed and which shows 
a mean free path close to the geometrical. 
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The fraction of inelastic scattering events is about 4°% of all interactions 
in flight, in agreement with the figure of 1-4°%, quoted by Goldhaber 
(1957 a) at the Venice—Padua Conference. The energy loss is in general 
large as was first pointed out by Alles et al. (1957 a) 


3.2. Characteristics of K~-Meson Absorption Stars 
As the mechanism for the absorption of a fast K~-meson by a nucleus 
is very different from that for absorption of a captured K--meson, we 
have divided the absorption stars found into these two classes. Their 
characteristics are presented in tables 3 and 4. 


Table 3. Characteristics of 815 K~-Meson Capture Stars are Shown 


All identified 


No. of No. of 
prongs stars at HE See SS oe ee 

0 140 0 0 Cetin aki patie 

1 141 36 1 te, iil ae <0) 

2 195 93 2 ie AW he IL 

3 161 51 7 G2) eon ea oiees 

4 82 25 10 Stays hen. SRL Ta) 

5 56 18 4 eee eee Ty 

6 28 5 3 PP rae 

7 12 2 1 Cpege0 Os) Ome) 

Totals 815 230 28 Bi RN Ty Et 
Percentages 100-0 28-2 Bde 14% 3-8) 14153-2501 ot 


The stars are divided into categories according to the total number of 
prongs. Columns 8 to 7 show the total numbers of 7+-mesons, hyper- 
fragments and &-hyperons in each category. Columns 8 and 9 show the 
total number of stars in each category which have both a &-hyperon and 
a 7-meson. 

In these tables, an identified hyperon is one with charge determined by 
its mode of decay or interaction at rest, or in some cases by conservation 
of charge if the same star gave a 7-meson which could be followed to rest. 
In a systematic study of 2-prong stars, Gilbert et al. (1957) found no cases 
in which both S-hyperon and 7-meson from such stars had the same 
charge. This result is confirmed by our data for cases in which it was 
possible to identify the hyperon by its decay or interaction. The columns 
labelled E+ contain the decays in flight, U+>n+7*, for which the sign of 
the hyperon is unknown. The column labelled X~ contains &~ capture 
stars with >1 prong only. 

Non-mesonic hyperfragments could not generally be identified because 
of their short range, but in some cases charge conservation is useful in 
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Table 4. Characteristics of 106 Fast K~-Meson Absorption Stars are 
Shown. The Form of the Table is the Same as that of Table 3 


All identified 


No. of No. of 7 

prongs stars aw HE See ae ee ee 

0 17 0 9) 0) i) 0 Q O 

i i O+ 0 0 1 0 ©) 0) 

2) 25 a 1 0 3 0 ] 2 

3 20 6 1 X 1 iL 0 

4 16 I 1 0 1 0 0 0) 

5 8 i 1 0 0 0 0 0 

6 4 1 0 ) ) 0 0 0) 

il 4 0) 1 0 0 0) 0 O 

| Totals 106 16 5 2 6 0 Y 2 
Percentages 100-0 15-0 4-7 1-9 --5°T - O-D NV ieoe ES 


+ See earlier comment in discussion of mean free path. 


distinguishing between X~-hyperon capture stars and true hyper- 
fragments. Figure 1 shows a range distribution of all possible hyper- 
fragments, and we conclude from this and the X~-hyperon energy spectrum 
of fig. 2 that no significant number of these ‘ hyperfragments ’ are actually 
=~ capture stars. 


Fig. 1 
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No. HYPERFRAGMENTS 


| 2 5 0 20 50 100 200 500 
RANGE (MICRONS - LOG SCALE) 


Range distribution of 28 hyperfragments emitted from 815 K--meson captures 
at rest. 


K~-Meson Interactions in Nuclear Emulsion 1199 


The probabilities for hyperon production which can be inferred from 
table 3 for capture stars, are 0-047 and 0-038 for identified X+- and 
~-hyperons respectively, rather lower than those usually quoted (see the 
references of §1). Such a situation may result from larger scanning loss, 
or more strict criteria for identification of events. In order to check this, 
we can determine the branching ratio of the S+-hyperon decay. The best 
estimate will be that determined from =+ decays at rest. The ratio 
(n7*)/(p7°)=13/10. This ratio is consistent with the figure of 1-15 + 0-29 
derived from the hydrogen bubble chamber data (Goldhaber 1957 b) and 
we feel that this is an indication that scanning losses will not invalidate 
conclusions drawn from the present investigationt. 

From table 3 it is seen that the ratio of the number of stars having both 
a &*-hyperon and az -meson to the total number of stars with ©+-hyperons 
is 26/38=0-68, whereas the proportion of all capture stars having 
7*-mesons is 0-28. An approximate comparison of these figures is possible 
if we assume that about 15% of all K~ absorptions occur through the 
two-nucleon mechanism (see later), and that the nucleus is self-conjugate 
and the interactions charge independent. Using the theorem of 
Watson (1952, eqns. 20-22), it is seen that the average escape probability 
for 7+-mesons from all capture stars is about 0-50. This figure includes the 
scanning loss and takes into account both X- and A°-hyperon production. 
The difference between these two figures, though hardly statistically 
significant, can be understood using either the surface or the volume 
absorption model if the absorptions leading to the Xz class of events occur 
further from the centre of the nucleus than in the average case. 

The energy spectra of the X-hyperons emitted from the K~ capture 
stars are shown in fig. 2. The energy spectra for &*+, &~ and &+-hyperons 
are shown in parts (a), (b) and (c) respectively. A fairly clear group of 
high energy hyperons appears above 60 Mev (see fig. 2), which can be 
taken as evidence for the capture of K~-mesons by two nucleons. This 
interpretation is also suggested by the fact that only one of these fast 
hyperons is accompanied by a w-meson. The lower energy hyperons 
will then result from the single nucleon absorption process with a small 
contribution from the two nucleon process owing to scattering of the 
hyperons by the residual nucleus. In the absence of knowledge of the 
hyperon-nucleon cross sections we can only assume that such a contri- 
bution is small; consequently we infer that the ratio of two-nucleon to 
single-nucleon processes producing escaping hyperons is 0-14 + 0-05. 

=+- and =--hyperons are scattered and absorbed equally in a self- 
conjugate nucleus if the interactions are charge symmetric. In this case, as 
pointed out by Gilbert e¢ al. (1957), the spectra of figs. 2 (a) and 2 (6) should 


+ No Dalitz electron—positron pairs were observed, although from 815 stars 
we would expect about 3 7°-mesons to have decayed in this way. However, 
no case of such a pair has been reported by other groups from a total of several 
thousand K- stars, indicating that the apparent non-occurrence of this mode 


may be significant. 
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be displaced in opposite senses by the Coulomb potential of the nucleus. 
The total shift of the two spectra will then be approximately twice the 
average Coulomb potential V,, effective in the region of capture. Such 
a shift is obvious in the figures, and leads to an estimate of this Coulomb 
potential of V..~6+4 Mev, in agreement with Gilbert et al. 


Fig. 2 


a 
S 
a 
ros 
8 

asl a 

() 
fs 
0 = a l a a L | 
20 40 60 80 100 


KINETIC ENERGY (MEV) 


Energy ods of _i-hyperons emitted from 815 K--meson captures at 
oat pds ee =*-hyperons ; (b) identified x--hyperons ; (c) 
- nesic decays ight 4 hi 
ae Pe : psie decays im flight, unidentifiable charge. Cross hatching 
lcates emission in association with a 7-meson. _ 


K~-Meson Interactions in Nuclear Emulsion 1201 


Figure 3 shows the distribution of the angle y between the initial directions 
of the hyperons and z-mesons from capture stars with no visible prongs 
other than these. In the histogram, we have distinguished between the 
contributions from positive and negative hyperons, but since the origin 
of the correlation is almost entirely kinematical they appear on one diagram. 
The curves in the figure show distributions of ~ which would be expected 
if the production of hyperons and z-mesons is isotropic in their centre 
of mass system, assuming two effective momentum distributions for the 
capturing proton and no interaction of either particle with the parent 
nucleus. The solid curve assumes a Gaussian momentum distribution for 


Fig. 3 


EVENTS 


, ia aS 
20 60 100 140 180 


INCLUDED ANGLE X (DEGREES) 


Distribution of the angle y between initial directions of =-hyperons and 
m-mesons emitted together from 815 K~ captures at rest. Solid cross 
hatching indicates =--hyperons, broken cross hatching, =*, no 
hatching, =+. The solid curve is the distribution predicted for a 
Gaussian momentum distribution for the capturing proton, the broken 
curve that for a Fermi gas distribution. 


the proton with average kinetic energy 19 Mev (Brueckner et al. 1955); the 
broken curve, a Fermi distribution with cut-off at 20 Mev. Ascan beseen, 
the correlation is remarkably insensitive to the assumed distribution, and 
either fits the data rather well (the validity of assuming such types of 
effective nucleon momentum-distribution is discussed in § 4). 
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= 


In fig. 4 is presented the energy spectrum of all protons greater than 
50 Mev in energy from the 815 K~-meson captures at rest. These protons 
are the result of +-meson absorption and of the two nucleon capture process. 
There do not appear to be any distinct groups, but the spectrum is given 
for general interest. 


Vig. 4 
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Energy distribution of all protons of energy greater than 50 Mev emitted from 
815 K--meson captures at rest. 


3.3. Characteristics of the X-hyperons 

The decay branching ratio (nz+)/(p7°) of the +-hyperon was estimated 
from our data in the preceding subsection to be 1.3. 

The lifetimes of the hyperons, and their mean free path for interaction 
and scattering are so sensitive to scanning biases and so dependent on 
details of interpretation that we do not consider nuclear emulsion a 
reliable tool for their estimation in quantitative terms. 

It is possible to search for non-invariance of the =-hyperon decay 
interaction under the operations of parity and of charge conjugation, as 
the effect sought should be independent of scanning biases. We must 
look for an asymmetry in the emission of z-mesons with respect to the 
direction of polarization of the hyperon. A strong asymmetry has been 
observed for A®-hyperons in bubble chambers (Feld 1957), but not for 
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charged &-hyperons. Gatto (1957) has suggested that the capture of 
K~-mesons by nuclei may produce polarized hyperons because the Fermi 
motion of the nucleons within the nucleus can provide a preferred direction 
to define a polarization. In 2-prong stars with only a hyperon and a 
m-meson visible, the direction of motion of the capturing proton is defined 
by the momenta of the hyperon and meson. Thus any polarization of 
the hyperon will be perpendicular to the plane defined by the initial 
directions of the hyperon and the meson, provided that parity is conserved 
in the strong hyperon-producing interaction. An asymmetry in the 
direction of the decay 7-meson owing to parity violation in the decay 
process will therefore occur about the plane determined by the L-z 
production. 

In presenting the results of these measurements, we use the convention 


(PsA Pz) ° Pip 


GUE i tap 
[PsA P, . P,p |sin x 


as suggested by Smith and Inman (private communication). P, and 
P_ are the initial momenta in the laboratory system of the associated 
=-hyperon and z-meson respectively, while P_,, is the initial momentum 
of the z-meson from the decay of the &-hyperon in the hyperon centre of 
mass system. Table 5 shows the distribution of this angle. 


Table 5. The Distribution of cos w’ is Shown for the Different Decay 
Modes at Rest (r) and in Flight (f) 


| Cos w’ > 0 cos w’ <0 


it+p+7° (r) 3 3 
Stepix® (£) 7 G 
Yt>n-+_at (r) #4 4 
Lt+n+7 (f) 4 2 


It is important to note that failure to observe an asymmetry in w’ 
may well be attributable to a lack of polarization of the &-hyperons with 
respect to the plane of production. Such polarization could not occur for 
an s-state capture of the K~-meson by the proton in the nucleus and the 
lack of asymmetry could be considered to be evidence for this possibility. 
In such a case, the suggestion of Treiman (1956), that the spin of the 
= -hyperon may be investigated by looking for an anisotropy in the direction 
of the decay 7-meson (in the hyperon centre of mass system) with respect 
to the initial hyperon direction, may be expected to be valid. This latter 
angle, 6’, may be also investigated for an asymmetry about 6’ = 90° which 
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would possibly occur if the parity doublet theory of Lee and Yang (1956) 
were correct. The distribution of cos 6’ found for all hyperons is shown 


in table 6. 


Table 6. The Distribution of cos 6’ is Shown. cos 6’ is Divided into 
Intervals of 0-500 


cos 6’ 1-0-0:5 0-5-0-0 0-0-(—0-5) | —0-5-(— 1-0) 
xut>p+7° (r) 4 4 0 4 
ut>+p+z7° (f) 6 2 9) 6 
x<t+n-+r+ (r) 3 5 3 Ss 
Meson + a= (f) i 6 S ps 

Fig. 5 
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Prong number distribution from 30 =~-hyperon captures at rest. No estimate 
of the number of zero prong stars is made. 


Finally, the =~-hyperon capture stars may be mentioned. For the 
reasons given in § 2 and in the discussion of § 4, we have not attempted to 
determine the number of zero prong capture stars; we merely record the 
characteristics of identified X~-hyperon stars. In fig. 5 the prong number 
distribution isshown. Figure 6 shows the energy distribution for particles 
of charge Z=1 (assumed to be protons) emitted from the capture stars. 
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Fig. 6 


25 


Ga 


No. OF PRONGS 


0 jtierag el eel lea 


20 40 60 80 00 ~=-:120 
PRONG ENERGY (MEV) 


No. OF PRONGS 
n 


4 8 V4 16 20 24 
PRONG ENERGY (MEV) 


Energy distribution for particles of Z=1 (assumed to be protons) emitted 
from 30 &~-hyperon capture stars. 


§ 4. Discussion 

A problem much discussed in the literature is the determination of the 
branching ratios for production of the various hyperons in the K~-meson 
absorption process. The difficulty lies in the identification of &~-hyperons 
by their endings, an unknown fraction of which produce zero prong stars. 
A further complication is the ordinary scanning loss for all types of events, 
especially those involving the emission of fast 7+-mesons. We can estimate 
the ratio, S, of the number of ~-hyperons produced by the reaction 
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K-4+n>2-4+7° sine ere, Regu em 


to the number produced by 

K-+p7>2-4+7t ere rd a) E3) 
with some rigour by making scanning loss (due both to individual observers 
and to zero prong events) cancel out, and also any correlation in direction 
of emission of the particles. 

Let N(+) be the number of &+-hyperons seen, N(+7) the number of 
x++a-, N(—) the number of X-, N(—7) the number of Dar For 
each particle, i, we can define a quantity x(7), being the probability of 
escape of particle i from the nucleus, and another, s(¢), the scanning 
efficiency for that particle. These last quantities may be found through 
observations of the reaction 

K-+p>it4+77 US. ey ee 


which, if we exclude those from the two nucleon interaction by their 
energy, is the only source of &+-hyperons. The number of &*-hyperons 
seen is then 
N(+)=P(3)x(2*)s(2*) 

and 

N( +7) =P(3)a(2*)s(Z*)a(n—)s(r-) F( + ) 
where P(3) is the probability of reaction (3) occurring and F(+) is a 
factor depending on the angular correlation of the hyperon and 7-meson 
in the laboratory system of coordinates. The last factor allows for 
correlations. which would be important if K~-meson absorption takes 
place in the nuclear surface. Thus 

a(n )s(7-).F( + )=N(+7)/N(-+ ). 

Similar considerations of reactions (1) and (2), which occur with prob- 
abilities P(1) and P(2), lead to 

N(—7) _ P(2)a(2-)s(27)a(art)s(at).F(—) 

N(-) (PQ) + PQ) (es) 
We now assume that x(7~)s(7-) F( + ) = a(a*)s(7*)F(— ), which is reasonable 
to the extent that 7-meson interactions are charge symmetric, scanning 
losses for 7+-mesons are independent of their charge, and that the correla- 
tion is kinematic in origin. Combining these results we have 


1+ P(1)/P(2) = N(+7)N(—)/N(+)N(—7)=14+S8. 


If the X-nucleon interaction is charge symmetric, it is easy to see that the 
contributions to (+) and N(—) by charge exchange processes approxi- 
mately cancel. The weighted mean value of S, compounded from the 
results of each individual observer, is 0-9+ 0-7, which may be compared 
with the value S=0-87 calculated, assuming charge independence, from 
the results of the isotopic spin analysis of Alvarez et al. (1957) on the 
absorption of K~-mesons in a hydrogen bubble chamber. This figure 
takes into account capture by both light and heavy nuclei in the ratio 3/7 


2 
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and the fact that the heavy nuclei have N#Z. As Alvarez et al. point 
out, however, in this case it is not known whether the absorption takes 
place from the ls or 2p Bohr orbit, and so the comparison could be 
meaningless. The value for S found by Dilworth (1957) for all nuclear 
emulsion data available at that time, and using similar methods is 0-6. 

We now propose to discuss the results of § 3 which have a bearing upon 
the absorption of K~-mesons by nuclei, and try to show that they favour 
the K~ surface absorption model. As was mentioned in the introduction, 
the calculations of Jones (1958) indicate that absorption of K--mesons 
captured in Bohr orbits by the Coulomb potential of heavy nuclei should 
take place for well defined large orbital angular momenta, if these orbits 
are circular (/=4 for Ag and Br; as low as /=2 for the light elements). 
Experiments designed to observe x-rays from the K--meson electro- 
magnetic transitions on this model should not show transitions involving 
Bohr orbits with angular momentum below about /=4 for medium—heavy 
nuclei. In addition, the very short mean free path of K--mesons in 
nuclear matter implied by the scattering and absorption of fast K~-mesons 
makes it a priort unlikely that any mechanism could introduce them far 
into the nuclear volume. 

The probability of absorption at a distance 7 from the centre of the 
nucleus is approximately P(r) = const. r?|%|%p(r), where & is the perturbed 
K— wave function and p(r) is the nucleon density distribution. P(r) has 
been calculated for two nuclei representative of light and heavy elements, 
for different values of the complex nuclear potential V, and for several 
different values of /+. The best estimate of V is Re(V)= — 20 Mev, and 
Im(V)= —50 Mev. The resulting distributions are plotted for the bromide 
(= 4) and for carbon (/= 2) in fig. 7. We take p(r) to be identical with the 
proton density distribution determined by high energy electron scattering 
experiments (Hofstadter 1956). V,, is the Coulomb potential. P(r) has 
been plotted for the most likely set of the possible parameters because the 
result is insensitive to reasonable variations of these. 

Absorption of the K~-meson in a region relatively remote from the 
nucleus immediately explains the observed correlation displayed in fig. 3 
between the directions of L-hyperons and z-mesons from capture stars: 
if the hyperons were produced inside the nucleus, the strong scattering 
to be expected of the hyperons would destroy the correlation. The 
relatively large probability of both hyperon and meson escaping the 
nucleus, suggested by the large proportion (0-68) of &+-hyperons seen 
with z--mesons is also hard to understand on the basis of the volume 
absorption ‘model. Finally, the absence of observed &~-hyperons in 
association with 7~-mesons indicates that the charge exchange reaction 


rrr 

+ The calculation of P(r) was made by Mr. J. M. Charap of Trinity College, 
Cambridge, with the kind permission of Dr. M. V. Wilkes of the Cambridge 
University Mathematical Laboratory, using the EDSAC 2 digital computer. 
Details of such calculations will be published soon by Mr. Charap. 
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Fig. 7 


ARBITRARY UNITS 
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Plot of the nuclear matter density distribution p(7), K~-meson capture prob- 
ability P(r) from the most probable Bohr orbit, and the Coulomb 
potential V¢, against the distance, 7, from the centre of the nucleus, 
for two nuclei representative of nuclei of photographic emulsion. See 
text for other parameters used in the calculations. 


X°+n—> 2X +p does not contribute appreciably to the number of observed 
2~-hyperons as we should expect if the absorption took place inside the 
nucleus. 

The only alternative explanation of these effects seems to be that the 
mean free path in nuclear matter of normal density for scattering and 
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absorption of &-hyperons is of the order of the nuclear radius. The 
estimate of Gilbert and White (1958) is that a nuclear mean free path for 
absorption of 5+% x 10-8 em (cross section 11+ ’ mb) would be required to 
reproduce the observed escape probabilities. Information on this point 
is rather scarce, but preliminary data have been reported (Venice—Padua 
Conference 1957) on scattering and absorption of fast E-hyperons which 
indicate at least geometrical cross sections for interaction. The existence 
of a X+-proton compound, recently reported by Baldo-Ceolin et al. (1957), 
with small binding energy, would imply a fairly large zero energy scattering 
length in either the triplet or the singlet spin state and consequently a large 
hyperon—nucleon scattering cross section at low energies. In nuclear 
matter, of course, the effective scattering cross section would be decreased 
by the effect of the exclusion principle. Results obtained with the 
Berkeley hydrogen bubble chambers (private communication from L. W. 
Alvarez) also indicate a large hyperon—nucleon cross section. A mean 
free path of greater than about 5x 10-3 cm for X-hyperons in nuclear 
matter, such as is probably required to make the results consistent with 
the hypothesis of volume absorption, therefore seems rather improbable. 

It is possible to reconcile the actual observed abundance of K~-meson 
absorption products with the surface absorption model. Using the 
m-meson interaction parameters given by Frank et al. (1956), and assuming 
isotropic emission of 7-mesons from points of origin distributed through 
the nuclear surface with probability P(r) shown in fig. 7, the escape 
probability of such z-mesons is found to be 0-66. If 85°% of K~-meson 
absorptions produce z-mesons (two-thirds of which are charged), we 
should observe 37% of absorption stars to have z-mesons. This is quite 
close to the generally accepted experimental figure, of a little over 30%, 
if a reasonable scanning loss is assumed. For &*-hyperons, the escape 
probability will be about 0-5 if the nucleus is very black to them, while 
for =~-hyperons, the Coulomb field will capture about 30% (the fraction 
of X+-hyperons with energy <2V,,) and the escape probability is reduced 
to about 0:3. Estimates of the X-hyperon production probabilities can 
be made using charge independence and the results of absorption of 
K--mesons in hydrogen (see Jones 1958); and using these our model 
predicts that 6° of stars will yield &+-hyperons and 12% &~-hyperons. 
These figures are consistent with experiment. Such very rough figures 
may well be misleading. We should like to point out, however, the 
simplicity of the assumption (strong hyperon—nucleus interaction) needed 
to reconcile all aspects of the experimental situation (y-distribution, 
a and = abundances) with the surface absorption model, and to contrast 
this with the special and probably implausible assumptions demanded by 
the volume absorption model. 

If further investigation confirms that the K~-meson absorption occurs 
in the surface of large nuclei}, it will be a reaction of rather unusual 


_ +The same considerations of course apply to the absorption of stopped 
anti-protons. 


P.M, 4N 
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properties. It will provide a tool for investigating the nature of dilute 
nuclear matter. It will also be complicated both by the fact that the 
K--meson has high angular momentum, and that at the low relative 
momentum involved, the strongly exothermic part of the total reaction 
cross section should be proportional to (k)~ for s-state interaction with a 
free nucleon, where k is the centre of mass momentum. Both these 
factors modify the effective momentum distribution of the capturing 
nucleon. It can easily be shown, using the independent particle model 
of the nucleus, that the probability of absorption of a K~-meson by a 
nucleon with ground state wave function ¢, depends on the final state 
momentum K through a factor 


Ke) fdr, exp (—1iK . rj)vb;,(ri) h(t 4) . 


d; is the wave function of the captured K~-meson, and the interaction 


operator assumed is 
A(r—r,;) =vd(r — rj). 


For the heavy nuclei, the K~-meson should, as we have remarked above, 
be absorbed from a state of high angular momentum. In this case, 
¢,(r;) <0 inside the nucleus, and there is no reason to expect that the 
effective momentum distribution will coincide with a nucleon momentum 
distribution found from the analysis of other high energy nuclear reactions 
occurring more uniformly throughout the nuclear volume. It may even 
be different for neutrons and protons considered separately. 

Finally, we may remark on the difficulties involved in the interpretation 
of the other experimental data, particularly the energy spectrum of fast 
protons and the number of =~-hyperons emitted from the K~ capture 
stars. In the case of the protons with energy greater than 100 Mev, 
the different processes mentioned in § 3.2 contribute and it is certainly 
not correct to assume that only the protons produced by the absorption 
of a K~ by a nucleon pair are important. The determination of the 
number of zero prong X~ capture stars is difficult. One possible method 
would be to compare the distribution of the angle y for those 2-prong stars 
which have no visible hyperon with the distribution for Xa stars shown in 
fig. 3. However, we regard this method as unreliable since at energies 
of more than about 10 Mev the direct interaction processes become 
important and it can no longer be assumed that the distribution in x for 
protons will be isotropic: they will to some extent carry the memory of 
the direction of the hyperon or meson from which they originated. 

In discussing the escape of the reaction products from the nucleus using 
an optical model, the surface interaction and the high angular momentum 
involved are both important. Because the interaction is in the nuclear 
surface, but close to the body of the nucleus, the escape probability is 
insensitive to the reaction cross section for geometrical reasons. Thus, 
exact calculation of the absorption branching ratio inside the nucleus by 


direct use of the results of Alvarez et al. (1957) would be questionable in 
the present case, 
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§ 5. CONCLUSION 


We feel that the present data on K--meson interactions in heavy nuclei 
indicate fairly clearly that absorption takes place in the nuclear surface. 
The simplicity of the assumptions needed to explain the data is a further 
point in favour of the model presented as compared to the complex and 
explicit assumptions needed to justify the volume absorption model. 
However, the model is so insensitive to the hyperon—nucleon interaction 
that little information on this point can be obtained from a study of 
K~-meson absorptions, and it must be investigated separately. 
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ABSTRAOT 


The yield and fracture stresses of unirradiated and irradiated En2 steel have 
been measured at 20°, —78°, and —196°c. From the change in yield stress, 
o, (shear), with grain size, 2d, the values of o; and k,, in a Petch type equation, 
o,=ai+k,d 4, have been determined. The term k, is unaffected by 
irradiation to within 5% and o; increases at arate approximately proportional 
to the cube root of the dose. The effect of irradiation on the grain size tran- 
sition from ductile to brittle fracture has been determined experimentally at 
—196°c and estimated from the change in o; and the value of k, for tests at 
—78° and 20°c. The results are in accord with the predictions of Cottrell’s 
relation (a; dt? +k) k,=Byw for the transition. The transition surface 
relating grain size, yield stress and temperature has been obtained. 


§ 1. InTRODUCTION 

COTTRELL (1958) has recently developed a theory of brittle fracture which 
overcomes some difficulties of previous theories and accounts for the 
brittle cleavage of single crystals. In this theory it is proposed that the 
crack nucleus is formed at the junction of two intersecting {110} slip 
planes by the coalescence of dislocations on these planes: 

ga{1 1 1] + 4a[i 1 1}+ [001]. . reyes (L) 
The new dislocation is a pure edge lying in the (001) cleavage plane and 
if further dislocations run into it a cleavage crack is formed with a length 
of the same order as the slip length. When the joint size of the crack and 
slip bands is greater than the Griffith value the crack will continue to 
grow and produce cleavage fracture, but if it is smaller than the Griffith 
value the crack will not grow beyond a limited size and will become 
blunted by plastic flow at its tip. A calculation of the critical size shows 
that in grains of diameter 2d the transition occurs when the relation 

(0,04? + k,,)k, = Byp Pret Foon beens eed 

is satisfied, where shear stress o,, grain diameter 2d, and k,, are related to 
the shear stress for yield o,, by a Petch type equation 
RRL eat es ree eek), 
B is a constant depending on the mode of deformation, being unity for 
uni-axial tension and one-third at the root of a plastic notch in hydrostatic 
tension, y is the effective surface energy of the crack, and yp is the shear 


modulus. 
GTS Sr as 1 


+ Communicated by the Authors. 
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The constant k,, is a measure of the stress required to operate dislocation 
sources and og, measures the stress opposing the motion of dislocations. 
Heslop and Petch (1956) have shown experimentally that o; can be 
divided into two terms 


0;=01+ 9% . ° . ° ° . . . . (4) 


here a, is the resistance of the lattice itself, and o,, the hardening caused, 
for example, by impurities, precipitates, and other dislocations in the 
lattice. 

For steels irradiated with neutrons, Churchman et al. (1957) have 
shown that hardening results from an increase in the frictional resistance 
o,. The results suggest that radiation hardening is caused by dispersed 
regions of radiation damage, such as clusters of point defects, that are 
sufficiently large to withstand thermally activated movements of 
dislocations. 

Equation (2) shows that the transition will depend on grain size, mode 
of deformation, and yield stress which is sensitive to temperature and 
strain rate. In the present work o, and k, have been measured for En2 
steel at three testing temperatures after irradiation doses up to 
9x 10!¥neutronem-?. The combined effect of grain size, temperature, 
dose and yield stress has been determined. Providing k, is unaffected by 
irradiation the results can be used to check the validity of Cottrell’s 
transition equation which predicts that if k, and Byy are constant the 
product of d? and o; should be constant at a given transition temperature. 


§ 2. EXPERIMENTAL 


Specimens were machined from En2 steel with a composition given in 
table 1. Double shouldered Hounsfield type test pieces with a lin. 
parallel gauge length and 0-126in. diameter were used. The yield and 
fracture stresses were measured at 20°, — 78°, and —196°c ina hard beam 


Table 1. Composition of En2 Steel 


C Si Ny 12 Mn Ni Cr Mo Co Al N 
YEN % % % % % % % % % 
0-15 | 0-12 | 0-051 | 0-026 | 0-59 | <0-15|] <0-10| <0-15 | <0-06 <0-02 | <0-011 


tensile machine with a strain rate of 10-2 per minute. An optical lever 
system magnified the beam deflection. 

A range of ferrite grain sizes (300 to 2850 grains/mm2) was obtained by 
heat treatment and all the specimens were tested in the annealed condition. 
The specimens were irradiated in the BEPO and Windscale graphite 
reactors. The positions in the pile were chosen such that the maximum 
flux in a temperature range 60-100°C was utilized. The specimens were 
irradiated to doses between 1018 and 9-1 x 10!8neutronem-2, The thermal 
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and epithermal fluxes were approximately equal and an estimated 10°, 
of the epithermal neutrons have energies greater than 1 Mev. 


§ 3. VARIATION OF k, AND o, WITH TEMPERATURE AND IRRADIATION 
Doss 

In figs. 1-4 the lower yield stress and fracture stress results for three 
temperatures and after doses up to 9-1 x 10! neutron cm~? are summarized. 
For unirradiated specimens k, and o; increase as the temperature of 
testing is reduced. The values, k, = 2-89 x 10’c.g.s. units, o,=4-48 x 108 
dynecm~ at room temperature agree reasonably well with those of 
Cracknell and Petch (1955), i.e. k, = 2-52 x 107¢.g.s. units, o,= 3-54 x 108 
dyne cm~? for En2 steel of similar composition. There is an approximately 
linear variation of k, with temperature, dk,,/07'= — 8-0 x 10*c.g.s. units. 


Fig. 1 
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Yield and fracture stresses of irradiated specimens at — 196°C. 
(a) For doses 1 x 1018 and 1-5 x 10!® neutron em-. 
(6) For doses 5-7 x 1019 and 9-0 x 10!° neutron em-, 


After irradiation it was found that k,, is unchanged, in agreement with 
Churchman et al. (1957) who reported that in tests at 20°C k,, remained 
constant to within 15%. Greater accuracy was not possible because of 
the scatter in the results which was mainly due to difficulties in measuring 
the grain size of irradiated specimens. This also applies to the results 
shown in figs, 1-4 except for the 1 x 1018, 3 x 10!8, and 8 x 108 neutron em? 
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dose irradiations tested at —78°0, in which the grain size was controlled 

more accurately. From the latter results it was concluded that k 

remains constant to within 5° with irradiation. ; 
Fig. 3 
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Irradiation increases oc, at a rate approximately proportional to the 
cube root of the dose, as shown in fig. 5. This agrees with work by 
Blewitt (1955) on copper. It is interesting to compare fig. 5 with the 
results of Heslop and Petch (1956) for the variation of o; with carbon and 


+ A somewhat better fit is obtained by taking the increase in oj proportional 
to fo4, 
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nitrogen content. They obtained a linear increase of o; which was 
independent of temperature. The irradiation results show a slight 
temperature dependence, e.g. 00;/0p1 = 2-93, 3°28 and 3-79 dyne cem~?/108 
neutron cm-2, at 20°, —78°, and —196°c respectively. 


Fig. 4 
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The fracture stresses at 20° and —78°c are not significantly affected by 
antes For tests at — 196°C irradiation rapidly embrittles specimens 
of all grain sizes and in many cases fracture occurs without any detectable 


yield. After high doses some of the specim : 
low stress value. pecimens fractured at a relatively 
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§4. Ducrmte-Brirrie Transition ox UNIRRADIATED SPECIMENS 


For the range of grain sizes and temperatures used all the unirradiated 
‘specimens showed some ductility and only specimens tested at — 196° 


Fig. 5 
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became brittle after irradiation. For tests at —196°c the reduction in 
area at fracture varied from 4 to 47% with the range of grain sizes used, 
but at —78° and 20°c the reduction in area was greater than 60% for all 
grain sizes. Since the coarse-grained specimens of this steel are ductile at 
—196°c the transition temperature cannot be determined without using 
lower temperatures. Low (1953) has given a method of determining the 
transition grain size which can be used for the present results. The lower 
yield stress o, and fracture stress o, are plotted against d—?. The 
o,/d—? plot follows the Petch relation (eqn. (3)) and o,/d~1? the 
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Griffith-—Orowan relation : 


ope const x da48 Vix eh 8 nD Me (5) 


which is a straight line passing through the origin. The intersection of 

o,/d-1? and o,/d~*? gives the transition grain size. In fig. 1 the method 
has been used to obtain the transition grain size at — 196°o, 
d-¥2— 17-3 em", o,/d-¥? = 0-2 x 10° c.g.8. units. 

The equation of the fracture stress line according to Cottrell’s theory is 


oe (Be )an Meee. i: 


y 

Taking »=8x104dynecm~, k,=4-6x 107c.g.s. units, and the results 
from fig. 1 the effective surface energy, y= 1-16 x 10tergscm™*. This value 
is less than half that obtained from results by Low, y= 2-8 x 10*ergs cm~, 
where k,=9-98 x 107¢.g.s. units and d-/?=11-2cm™, for a low carbon 
rimming steel containing 0-07% carbon, 0:03°% manganese, 0-02% phos- 
phorus, 0:033% sulphur, 0-020% oxygen, and 0-005% nitrogen. 

For tests at — 78° and 20°c the line plotted through the fracture stress 
results does not extrapolate to the origin. This is due probably to a 
change in the mode of fracture, from ductile cleavage to ductile fibrous 
(Crussard et al. 1956), in specimens tested well above the transition 
temperature. In fractures with a reduction in area greater than 50% the 
fracture line has been observed to deviate from linearity (Cottrell 1958). 
Thus the transition grain size cannot be determined directly by the 
method used for tests at —196°c. The transition grain size can be 
estimated, however, using eqn. (6). Assuming By to be independent of 
temperature in the range used, the slope o,/d~¥? can be calculated at 
— 78° and 20°c from the measured values of k, at these temperatures, and 
the value of By at —196°c. At —78°c, o,/d-V?=2-51 x 10%c.g.s. units 
and at 20°0, o,/d-?=3-21 x 108c.g.s. units. These values have been 
plotted in figs. 3 and 4 to obtain the transition grain sizes. 


§ 5. Errrect or [rrapIATION DosE ON TRANSITION GRAIN S1ZE 


It is now possible to determine the effect of irradiation on the transition 
grain size directly from figs. 1-5. Since £, is constant with irradiation the 
fracture line is unaffected and the transition grain size is the intersection 
of this line and o,/d~4? at various doses. In table 2 the information from 
these plots is summarized. 

As mentioned in §1, Cottrell’s equation for the transition predicts that 
if k, is constant, then at any one temperature, the product o,d'? should 
be constant. It has been shown that k, is constant with irradiation to 
within 5% and o; increases, thus d¥2 should decrease at a corresponding 
rate. Only the results obtained at —196°c can be used to check this 
prediction since the transition grain sizes at — 78° and 20°c were calculated 
assuming Cottrell’s relation to be true. From fig. 2 it is possible to 
determine the transition grain size for the high doses fairly accurately, 
but after 1x10'§ and 1-5x10!neutroncm- the results are rather 
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Table 2. Values of k,, 0; and Transition Grain Size at Three Testing 
Temperatures for Specimens Irradiated 
Between 0 and 9-1 x 10! neutron em~* 


Temp. of | Dose, ¢ | Dose, ¢1/3 hy oj Transition | o, d'/2 
testing | n.cm~?| n.cm-? ¢.g.s. |dynecm-?| grain size C.2.8. 
°c 10 x 108 x 107 x10? a, cust cl xe) OF 
—196 0 0 4-64. 2-66 17-5 1-52 
1 1 4-64 2-98 20-0 1-44 
15 2-47 4-64 3°61 26-0 1-39 
57 3°85 4-64 4-15 30-0 1-38 
90 4-48 4-64 4-40 32-0 1-38 
—78 0 0 3°70 0-72 3°5 = 
i 1 3°70 0-92 4-4 =S 
3 1-44 3°70 1-07 5-1 = 
8-6 2-05 3°70 1-27 6-0 = 
18 2-62 3°70 1-62 7 = 
91 4-50 3-70 2°22 10-5 a 
20 0 0 2-89 0-45 1-6 — 
1 1 2-89 0-62 2-2 = 
15 2-47 2-89 1-0 3°5 = 
57 3°85 2°89 1-45 5:0 == 
90 4-48 2-89 1-87 6:5 — 


» @ DR CHeals (ao 


QO; dynes em>*x 10? 


GRAIN SIZE d 2 cm 2 


Transition surface relating yield stress, grain size and temperature, 
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scattered and the transition grain size values obtained are only approxi- 
mate. The o; results are shown in fig. 5, for doses up to 1-5 x 101 neutron 
em-2, The high dose values have been obtained by extrapolation. The 
product o,d which is given in the last column of table 2, remains 
reasonably constant as o; increases. At —196°c there is a linear variation 
of transition grain size with the cube root of the dose (fig. 6). 

In conclusion, the combined effect of grain size, o;, and temperature on 
the transition can be obtained from the above results and the transition 
surface separating brittle and ductile fracture is plotted in fig. 7. 
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ABSTRACT 


A study based on electrical resistivity measurements at liquid nitrogen 
temperature has been carried out on the isochronal annealing out of quenched- 
in vacancies in 99-995% Al. The process of annealing takes place in two 
stages; the first stage is centred at about room temperature and takes place 
at a lower temperature the higher the quenching temperature; also the 
observed migration energies (in the range 0-37 to 0-58 ev) and the mean 
number of jumps for annihilation (in the range 104 to 10°) are very sensitive 
to the initial concentration of defects, that is to the quenching temperature. 
The second stage takes place at about 160°c and its amount increases 
largely with quenching temperature; its activation energy results in 
130-2 ev. The results have been interpreted as due principally to a 
clustering of vacancies (first stage) and subsequent dispersion of clusters by 
self-diffusion (second stage). 


§ 1. INTRODUCTION 


Ir is now a well-established fact that in many close-packed metals it is 
possible, by rapid quenching from high temperature, to freeze-in lattice 
defects which are in thermal equilibrium at the high temperature. 
Theoretical considerations on the energy of formation support the 
hypothesis that such defects should be vacancies (single or possibly pairs) 
rather than interstitial atoms. Cottrell (1957) has given a short review 
of some recent results for various metals. 

After the vacancies are quenched-in the question arises about the 
way they are annealed out to reach the equilibrium concentration. The 
present investigation was undertaken to throw some light on the annealing 
process in aluminium; the concentration of defects has been evaluated 
by means of the increase of resistivity they produce. In this paper 
results obtained specially by isochronal annealing are described; results 
connected with isothermal annealing will be reported in a forthcoming 
paper. 

The possibility of freezing-in thermally generated defects in aluminium 
was shown by Maddin and Cottrell (1955) by the increase of yield stress 
of aluminium single crystals by quenching and by ageing after quenching ; 
Levy and Metzger (1955) confirmed it by the decrease of the internal 


+ Communicated by the Authors. A short preliminary report on this work 
was presented in Italian by T. Federighi and F. Gatto, at the Meeting of the 
Societa Italiana di Fisica, at Padua, September 22, 1957, 
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friction of aluminium by ageing after quenching. In each case the observed 
effects were referred to a pinning of dislocations by vacancies. The 
assumption that the quenched-in defects were vacancy-like defects was 
then supported by the observation of Takamura (1956) who reported 
that during the annealing process a contraction of the sample was observed. 
Subsequently, the possibility of freezing-in vacancies in aluminium was 
shown by an increase of electrical resistivity by quenching (Wintenberger 
1956 and Panseri et al. 1957). The former showed that the rate of annealing 
of vacancies was largely increased by 10% plastic deformation ; the latter 
reported that probably some clustering of vacancies takes place in the 
isothermal annealing at room temperature as the process was quicker 
the higher the quenching temperature. 

Since the present investigation was started, other results obtained by 
electrical resistivity measurements have been published. Bradshaw and 
Pearson (1957) reported values of 0-76 ev for the formation energy of 
the quenched-in defects and 0-44 ev for their migration energy ; De Sorbo 
and Turnbull (1957) reported a value of 0-37 ev for the migration energy 
and quoted a value higher than 0-6 ev for the formation energy. Finally 
Wintenberger (1957 a) reported a value of about 0-57 + 0-01 for the migra- 
tion energy after air quenching. Till now however, as far as we know, the 
details of the annealing process after quenching had not been investigated. 


§ 2. EXPERIMENTAL PROCEDURE 


Aluminium of nominal purity 99-995°% was used. Principal impurities 
were Cu=0-0003%, Fe=0-0004%, Si=0-0015%. 

Owing to the very small resistivity variations these were always 
computed from resistance variations of wire samples of ¢ 1 mm and length 
50 cm with four soldered contacts. To simplify operations, samples were 
wound as a lamp coil filament and supported by an insulating rod. After 
soldering, samples were always stabilized by annealing some hours at 
630°c. Usually resistance measurements were carried out in liquid 
nitrogen using a dummy sample to account for the variations of temper- 
ature; only some additional measurements, as specified hereafter, were 
carried out at 20°c in a thermostatic oil bath. The usual potentiometric 
system, using a Wenner type potentiometer, was employed in each case. 
The sensitivity of the equipment was about 1 x 10-5, 

The resistivity variations have been computed by the formula Ap=pAR/E 
where AR and F are respectively the variation of resistance and the 
resistance of the fully annealed specimen, and p the resistivity of aluminium 
(p= 0-228 pohm cm at —195°c and p=2-65 pohm em at 20°C). Results 
have been expressed usually in muohm em = 10-° ohm em, a unit which 
is very useful for quenching experiments. 

Quenching operations were carried out by manual extraction of samples 
from the furnace and by rapid immersion in the cooling medium (brine 
at about 2°c); immediately after that, samples were dipped in liquid 
nitrogen. This procedure, which is not, of course, the most convenient 
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for accurate reproducibility of quenching operations, was selected in this 
research for its great simplicity. By this technique, wires thinner than 
1 mm could not be used as they would cool too much during the extraction 
from the furnace, before the water quenching. 

The annealing treatments, after quenching, were carried out at various 
temperatures in the range — 100 to 240° by using proper liquid baths. 


§ 3. REsuLts 
3.1. Formation Energy 


Although it was not the principal object of this work, it was necessary 
at the beginning to determine the increase of resistivity which was possible 
to obtain by the adopted technique, with the quenching temperature. 

Preliminary experiments had shown that by annealing a quenched 
sample for a few minutes at 240°c the same value of electrical resistance that 
existed before quenching was within experimental error obtained (see, 
however § 3.5) ; namely, all the excess of resistivity was practically annealed 
out at a such temperature. Therefore the quenched-in component of 
resistivity was usually computed from the changes of resistance, at 
liquid nitrogen temperature, after quenching and after the subsequent 
annealing at 240°c. 

Results of many observations are shown in fig. 1, plotted as usual, 
on a semilogarithmic scale against 1/7’; in the same figure the results of 
Bradshaw and Pearson are indicated. 

The following points may be noted: (1) notwithstanding the relatively 
large wire diameter, the quenched-in component of resistivity Apo 
increases with quenching temperature in the range 360 to 620°C; (2) at 
low temperatures, until about 420°c, a straight line can be drawn which 
gives H,=0-76+0-03ev for the formation energy in agreement with 
Bradshaw and Pearson’s results; (3) at higher temperatures a deviation 
below the straight line is observed; (4) the whole results are lower by a 
factor 0-65 than the Bradshaw and Pearson ones; the straight line of the 
present results can be described by the relation: 


Apo = 800 exp (—0-76/k7') wohm cm Steel er enol) 


and the pre-exponential factor is to confront with the value 1240 pohm cm 
reported by Bradshaw and Pearson. 

It is interesting to observe that although a fraction of defects is lost 
by using a large diameter wire as in the present case, the fraction of 
quenched-in vacancies is still sufficient to leave a role to vacancies in 
increasing the rate of diffusion in pre-precipitation in aluminium alloys 
(Federighi 1958). 

3.2. Isochronal Annealing 

In order to know in which range of temperature the defects were 
annealed out, a first set of isochronal annealing was carried out after 
quenching from several temperatures; the selected procedure for these 
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Semilogarithmic plot of quenched-in resistivity with the reciprocal of the 
absolute quenching temperature. 
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Isochronal annealing of 2 minutes at each temperature, after quenching from 
various temperatures. 
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tests was an annealing of 2 minutes at increasing temperature with steps 
of 20°C, in the range between —100°c to 240°; hence after each 
treatment the variation of the electrical resistance in liquid nitrogen was 
determined. 

This first set of results is shown in fig. 2; they have been normalized, 
namely Ap/Apy is given instead of Ap. Ap is the residual increment of 
resistivity after the treatment at each temperature and App is the experi- 
mentally observed increase of resistivity (see table 1), which of course 


for higher quenching temperature is lower than the value Ap, obtained 
from eqn. (1). 


Table 1. Data for Curves of Figs. 2 and 3 (Ap, and Ap,” are Mean Values, 
except for 360°C) 


uenchin - 
‘Saeed Apo Apo Apa’ Em Jumps 

°6 muohm cm | muohm cm Apo ev n 
600 19-5 3°5 0-18 0:58 + 0:03 1x 104 
550 Lig 1-55 0-14 0:42 + 0:02 6:6 x 10" 

470 4-7 0:5 O-l 0:37 + 0-07 6:5 x 10° 

400 1-6 ~0 ~0 0-45 + 0:03 1-4x 109 
360 0-7 if) 0 — 


The first general observation is the existence of a very strong influence 
of quenching temperature. 

In particular the following facts can be deduced: (1) for high (600°c) 
and mean values (470°c-550°c) of the quenching temperature, the 
-annealing of defects takes place in two stages; in consequence of these 
results the quenched-in resistivity Ap, can be subdivided into two parts: 


Apo=Apo + Apo’ 


where Ap,’ and Ap,” represent respectively the first and the second stage ; 
(2) for low quenching temperatures (400°C, 360°c) the recovery takes 
place in one stage; that is Apy)’~0; (3) the fraction of recovery which 
takes place in the second stage increases with Ap», namely with quenching 
temperature (see § 4.3.1); (4) the first stage which is centred at about 
room temperature takes place at lower temperatures the higher the 
quenching temperatures; for instance, after quenching from 600°C a 
fraction Ap/Ap,=0-6 is observed at about —30°c; after quenching from 
400°c the same fraction of recovery is observed at about 55°c; (5) the 
range of temperature of the second stage (if existing) does not seem to be 
practically dependent on the quenching temperature, and it is centred 
at about 160°c. 
The presence of a second stage in the annealing out of quenched -in. 
defects in aluminium had been noted earlier by Panseri et al. (1958) in 
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a study of Al-Mg alloys. Also Wintenberger (1957 b) reported that after 
a rapid quenching, only a fraction of quenched -in resistivity was annealed 
out at room temperature. 

Following these results other similar isochronal annealings, after quench- 
ing from four temperatures, were carried out, with the same procedure, to 
test reproducibility and the influence of time of annealing (2 or 15 min). 
Results are reported in fig. 3 (results of fig. 2 also are included). In table 1, 
the mean values of Ap, and Ap,” and their ratio for curves of figs. 2 and 3 
are reported (for other data of table 1, see below). 


Fig. 3 


Al 99,995 % 


L 


=O0 


0 50 100 150 MG 250 -100 250 0 50 100 150 


Isochronal annealing of 2 or 15 minutes at each temperature, after quenching 
from various temperatures. 


Since at the beginning of the anneal from —100°c the resistivity 
remained constant until a sufficiently higher temperature, experimental 
points in fig. 3 (and fig. 2) are reported only from that toupee at which 
the variation began to be observed. It is necessary to note however 
that after quenching from 400°c (and 360°C) an abnormal scattering of 
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results, sometimes with a small increase of resistivity, was noted before 
the true variation shown in figs. 2 and 3. At present it is not sure if this 
increase is real, and if it is, the reason is not clear. 

The reproducibility of curves of fig. 3 appears to be good enough for the 
high temperature quenching and not so good for the lower temperatures ; it 
should be noted however that in figs. 2 and 3 the scattering of experimental 
points is higher, the lower the quenching temperature, due to the large 
decrease of Ap, (see table 1). 

For the interpretation of the results it is useful to introduce the number 
of jumps 7 for annihilation of defects in the first stage, given by (Cottrell 
1957): 

m~12.¢.10% exp (—#_ /kT). eames femed 2) 


where H,,, is the migration energy and ¢ the lifetime of the defects (the 
entropic factor has been put equal to unity). From the data of fig. 3 it 
is possible to deduce the values of nm and of #,, They are evaluated 
for a fraction of recovery Ap/Ap,=0-6, and included in table 1. The 
values of #,,, are practically constant for any selected fraction Ap/Ap, in the 
first stage of recovery for each couple of curves drawn in fig. 3; the errors 
reported in table 1 have been computed from the scattering of experi- 
mental points. 

A noticeable dependence of #,, and n.upon quenching temperature 
appears to exist; this dependence, notwithstanding the indetermination 
of the results, can be concluded to be real and will be discussed in § 4. 

It is interesting to note that the rapid annealing observed after quenching 
from 600°c is related to the low value of and not to a low migration 
energy (which, on the contrary, is the highest observed). 


Mm 


3.3. Some Isothermal Results 


The effect of quenching temperature observed in isochronal annealing 
should be evident also in isothermal annealing curves as reported by 
Panseri et al. (1957) by means of electrical resistivity measurements 
at room temperature. To show this effect better, a set of isothermal 
annealing curves at 0°c was therefore obtained, by means of more 
convenient liquid nitrogen measurements. 

The quenching temperatures were selected in the range 460 to 600°C. 
After a time of isothermal recovery sufficiently long to get a large fraction 
of recovery, samples were left to rest 16 hours at room temperature (20°C) 
to deduce the final value of the first stage ; then the samples were annealed 
4 hour at 240°c to anneal out the second stage. 

Results, again in normalized form, are reported in fig. 4; the effect 
of quenching temperature and the presence of the second stage (shown 
by the limit value of each curve) are evident. 


3.4. The Second Stage 


3.4.1. Lsochronal annealing 
Some characteristics of the second stage will now be reported in this 
and in the following sub-sections. 
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In fig. 5 (left) isochronal annealing curves of the second stage are 
shown ; they have been deduced from the data of fig. 3, after quenching 
from 600°C; the results are normalized, namely Ap”/Ap," has been 
plotted, where Ap” is the residual increase of resistivity after each treat- 
ment, and Ap,” has been evaluated at 80°o (table 1). Similar isochronal 
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Isothermal annealing at 0°o after quenching from various temperatures. 
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annealing curves of the second stage, again after quenching from 600°c, 
are reported in fig. 5 (right); in these measurements the first stage has 
taken place in 16 hours at 20°c followed by 3 minutes at 80°0. 

The examination of results shows that the way in which the first stage 
takes place has a slight influence on the range of temperature of the 
second stage; namely, the second stage is not independent of the first. 
For example, a fraction Ap”/Ap,” = 0-6 in 2 minutes is observed at 170°C 
in fig. 5 (left) and at 185°o (right). 

From the shifting of the isochronal annealing curves with the time of 
annealing, an activation energy @ can be evaluated assuming that the 
time to observe a fraction of recovery is proportional to exp (Q/k7T). The 


result is @=1-:3+0-2ev, independently, for the data reported in fig. 5, 
left or right. 


3.4.2. Dependence of Ap” upon Apo 


The results obtained in many tests about the dependence of Ap,” on 
Ap, are shown in fig. 6 on logarithmic scales. Experimental points bave 
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Logarithmic plot of 4p,” with Apo (for explanation of experimental 
points see the text). 
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been subdivided into three classes; black circle points have been deduced 
from the data of fig. 3 (or table 1), white triangle points have been deduced 
from the isothermal curves of fig. 4, while white circle points have been 
deduced from proper experiments in which the first stage has taken place 
in 3 minutes at 80°c. The examination of results shows that despite a 
noticeable dispersion of experimental points the second stage Ap,” is 
related to the total increase of quenched-in resistivity Ap) by the law 
Apo” =a. Ap,” where m has a value of about 3 and a=3-8x 107". Also 
the second stage is somewhat higher if the first stage has taken place at 
low temperatures; again there seems to exist a slight dependence of the 
second stage from the way the first stage has taken place. 

The ratio Ap,"/Ap, is obviously given by 3-8 x 10~*4/Ap, and increases 
with App. 


3.4.3. Dependence of Ap,” upon 1/T 


The data of Ap,” of fig. 6 can be plotted against 1/7’ to deduce a 
fictitious formation energy for the defects of the second stage, as shown 
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The fictitious formation energy so deduced is about 1:2 +0-lev; about 
the same value could obviously be evaluated from the reported 
dependence of Ap)” vs Apy and from the dependence of Apy upon 7’. 

It is to be pointed out that the dependence of Apy” on quenching 
temperature is relatively strong. 


3.4.4. Effect of quenching speed 


The values of the second stage Ap)”, and of the total increase of 
resistivity Apo, are very sensitive to quenching speed, as shown in table 2. 
In these experiments the second stage has been evaluated after annealing 
for 3 minutes at 80°c. The time to reach a temperature lower than 80°c 
was less than 1 minute with air quenching and less than 4 minute for a 
stream of air; after this time samples were dipped in liquid nitrogen. 


Table 2. Influence of Cooling Medium on Ap, and Ap,” 
(Results are Mean Values of Several Observations) 


Cooling Apo Aas 
medium mpyohm em | myuohm cm 


Brine at 2°o 19-5 
Water at 30°o 

Oil 

Blow of air 

Still air 


3.4.5. Dependence of the second stage on the temperature of measurement 

As the second stage can be evaluated on a sample by measurements in 
liquid nitrogen, or at 20°c, control is possible if it is independent of 
temperature. 

A series of 6 tests has given the following result : A*p9"/Apo” = 1-17 + 0-03 
where the star indicates the measurements at 20°c. A slight dependence 
from the temperature of measurement seems to be present. 


3.5. Plastic Deformation during Quenching 


It is possible to obtain a rough estimate of the amount of plastic 
deformation during quenching by electrical resistance measurements. 
For this purpose it is necessary to evaluate the variation of resistance of 
the sample between the fully annealed state before and after the quenching 
(when all the introduced defects are annealed out). 

The resistance R of the sample can be expressed by R=p.f where pis 
resistivity and f a form factor (dimensionally f is a length Bt) during small 
plastic deformation f is altered to f+ Af, so that it should be possible to 
detect a permanent increase of resistance AR, given by AR,=p. Af. 
The plastic deformation ¢ should then be of the same order as Af/f=AR,/R. 
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As stated before, by carrying out measurements at liquid nitrogen 
temperature all quenched-in resistivity appears to be annealed out, 
between experimental errors, after the annealing at 240°c; therefore Af/f, 
if not zero, is very small. However, by carrying out measurements at 
20°C (before quenching and after a full annealing after quenching) a 
permanent increase of R can be detected more easily, due to the fact that 
AR, increases with p, that is, with the temperature of measurement. 

Experience has shown that AR,/# is scattered in the range 1-5 x 10-4 
and that there is not an appreciable dependence upon quenching tempera- 
ture in the range between 400 to 600°C. This very small plastic deformation 
cannot account, of course, for the increase of resistivity by quenching or 
justify the dependence of the annealing process upon quenching 
temperature. 


§ 4. DIscussION 


4.1. The First Stage 


4.1.1. High temperature quenching 


The presence of a double step in annealing and the variation of the 
migration energy with quenching temperature show clearly that more 
than one kind of defect should be present. Probably several tentative 
interpretations can be advanced; here we limit ourselves to show that it 
is possible to give a simple picture of the annealing process in terms of 
clustering of vacancies. At present this picture is only of a qualitative 
nature ; really it is not easy to give an exact quantitative picture of the 
annealing process of quenched-in defects (for instance, see Koehler e¢ al. 
1957 b) as too many factors are involved, especially if there is a tendency 
to clustering; for instance, one should know the migration energies for 
each kind of defect (vacancies, divacancies, trivacancies, and more complex 
aggregates), their binding energies and their contribution to resistivity. 

The present picture is based principally on three independent hypotheses : 
(i) the quenched-in defects are vacancies; (ii) an appreciable binding 
energy between vacancies exists; (iii) the contribution to resistivity of a 
cluster of vacancies is smaller than the contribution of the same dispersed 
single vacancies. 


Then the annealing process, in the case of the high temperature 
quenching (600°C), can be synthetized as follows. 

As soon as the annealing temperature increases enough to allow 
sufficient mobility of point defects, the concentration of which should be 
of the order of 104, vacancies begin to couple together, forming divacancies, 
with a number of jumps of the order of 104, which is much smaller than the 
number of jumps to reach dislocations (~ 10°); if the binding energy of 
a divacancy is sufficiently high it can live long enough to meet another 
vacancy before dissociation, and form a trivacancy more stable than a 
divacancy; the creation of trivacancies will be slowed up when their 
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concentration has increased to reach about the residual concentration of 
vacancies ; from this moment trivacancies begin to be reached directly by 
_ other vacancies and grow up to larger clusters ; therefore if the concentra - 
tion of trivacancies is large enough, only a few vacancies will find a way to 
dislocation. In conclusion, if the initial concentration of defects is 
relatively high and the interaction energy between vacancies sufficiently 
strong, one should expect at first, in isochronal annealing, principally a 
coagulation of vacancies in clusters of some unspecified shape and not a 
true annihilation of vacancies at dislocations or other sinks ; these clusters 
should be then annealed out subsequently at higher temperature by 
normal self-diffusion. 

This picture is in agreement with the observed results as long as we 
suppose that the first stage, observed in the range from — 80 to — 20°60, 
is caused to a large extent by a smaller effect of clusters on resistivity 
and not to a true decrease of the concentration of vacancies. There is 
really an expectation for a smaller effect of clusters on resistivity 
(Dexter 1956). The second stage should then be due to the dispersion 
of clusters. 

The experimental migration energy (0-58 ev) observed in the first stage 
should probably be connected with the motion of single vacancies, in 
agreement with Wintenberger’s (1957 a) result; really, if divacancies move 
much quicker than vacancies (Bartlett and Dienes 1953) they can reach 
another vacancy, to form a trivacancy, as soon as they form, and therefore 
the rate of clustering should be controlled by the migration of single 
vacancies. 

The suggested interpretation can account easily for the presence of the 
two stages, for the low value of jumps n (~ 10+) observed in the first stage, 
and the value of Q of the second stage (1-3 + 0-2ev) which is of the order 
of magnitude expected for self-diffusion (1-4+0-lev) in aluminium 
(Nowick 1951). 

It is also interesting to note that H,,+ H,=0-58ev +0-76ev=1-34ev is 
again of the expected order of magnitude for the activation energy for 
self-diffusion in aluminium. This fact may be considered as good 
evidence that in aluminium diffusion takes place wa vacancies, or 
vice-versa that the quenched-in defects are vacancies. 


4.1.2. Lower quenching temperatures 

A very important point to be considered is what happens when the 
initial concentration of vacancies is lowered, as by quenching from 550°c 
and 470°c. 

There is no doubt that the isochronal annealing curves must change. 
Actually at a given temperature the time for clustering is higher the lower 
the initial concentration of vacancies; therefore to observe a given 
fraction of recovery, in a fixed time, it is necessary to increase the annealing 
temperature; but at higher temperature the divacancies and the 
smallest clusters are not so stable and dissociation processes begin to be 
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competitive. The result is that the process of coagulation is relatively 
slowed down and many vacancies succeed in reaching dislocations giving 
a straight contribution to the decrease of resistivity ; one should observe 
therefore an increase in the mean number of jumps n, a decrease of the 
fraction of the clustered vacancies and a decrease of the mean migration 
energy due to the presence of a relatively higher fraction of divacancies 
(whose number of jumps to reach another vacancy is now higher). These 
facts are really observed. 

If the concentration of defects is again lowered, as by quenching from 
400°c, the same picture as above takes place, but the temperature at 
which one could observe (in a fixed time) the coagulation is now 
sufficiently high so that, due to the more rapid dissociation, no cluster 
can form, and also the concentration of divacancies is relatively smaller ; 
therefore one should observe the presence of no cluster, an increase in the 
mean migration energy (as the relative concentration of single-vacancies 
is now preponderant) and lastly, a high value of n, as now the whole 
decrease of resistivity should be due to a true annihilation of vacancies at 
dislocations. 

To conclude, this picture is qualitatively in agreement with actual 
experimental results. The migration energy of divacancies should be, of 
course, not higher than 0-37 ev, to justify the result obtained by quenching 
from 470°c; this value is to be compared with the value 0-58 ev assumed 
for single vacancies. 


4.2. The Second Stage 


As said before, the second stage should be due to the dispersion of 
clusters of vacancies which have originated during the first stage. 
This interpretation is in agreement with the presence of the two stages 
and with the activation energy @ for the second stage considered to be 
equal to the activation energy for the self-diffusion in aluminium. 

Of course, by resistivity measurements only, it is not possible to decide 
the shape of these clusters, for example, whether they are spherical or 
whether vacancies prefer to condense to form stacking faults (see, for 
instance, Kuhlmann-Wilsdorf 1958). 

Coulomb and Friedel (1957) have recently given a theory for the forma- 
tion of clusters of vacancies along dislocations. Although dislocations are 
surely preferential places where clusters are formed, the present results 
seem to show that clustering takes place anywhere in the bulk of the 
sample ; actually with the mechanism of Friedel one should expect always 
n~10° in the first stage and Q lower than the activation energy for 
self-diffusion of aluminium in the second stage. 

However, at present one cannot establish clearly the role of dislocations 
as nucleation agents for clusters; this point will be probably better 
understood after a study of the influence on the first stage of small plastic 
deformation after quenching and also by the influence of plastic 
deformation during quenching using wires of different diameter. In the 
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same way it is not possible at present to establish if the clusters are 
nucleated near some kind of impurity ; the use of zone-refined aluminium 
could probably resolve this. 

It will be shown now that other interpretations of the second stage are 
improbable. For example, it might be due: (i) to a precipitation of 
impurities set in solution before quenching ; (ii) to a fraction of vacancies 
trapped near some impurity or dislocation; (iii) to the annealing out of 
dislocations introduced by plastic deformation during quenching 
operations. 

With the first hypothesis a saturation for Ap,” should be observed at 
about 400—500°c, when all impurities are set in solution; besides a scarce 
dependence of Ap,” from quenching speed should result. There is no 
agreement with experimental results. 

With the second hypothesis a lower value for Q of about 0-6-0-9ev 
should be expected (namely, the sum of migration energy 0-58 ev and the 
binding energy with impurity atoms which would not be higher than 
0:2-0-3ev); also a quite different dependence should be expected of 
Ap,” upon 7’; if the concentration of sinks were high, Ap,” should be 
nearly equal to Apo, at least until a sufficiently high quenching temperature 
is reached; if the concentration of sinks were small, Ap,” should be 
practically constant with quenching temperature. 

As for the last hypothesis there is no doubt that during the quenching 
operation there is a plastic deformation, but this is too insignificant to 
account for the variation of resistivity of the second stage; for instance 
applying a formula of Wintenberger (1957a) a plastic deformation of 
5 x 10-4 should give an increase of resistivity of about 1-8 x 10-3? muohm cm, 
which is much smaller than the values of about 2 muohm cm experimentally 
observed (see table 1). Besides, the second stage should not depend very 
strongly upon quenching temperature, as in fact it does. Only the 
hypothesis of clustering therefore appears to account for the experimental 
results. 

4.3. Additional Remarks 
4.3.1. On the annealing after irradiation 

McReynolds et al. (1955) have reported that in aluminium irradiated 
with neutrons at —150°o a recovery of resistivity takes place in the 
range from — 80 to — 30 with a migration energy of 0-55ev. Several facts 
suggest that this recovery might be due to the clustering of vacancies in 
the same way as after quenching: (1) the migration energy is about the 
same; (2) the number of jumps computed from the data of McReynolds 
et al. (fig. 3 of their paper) is 1-4 x 10%, which is lower than the value 
observed after quenching from 600°C; this fact was expected as it is 
possible to deduce, from the same data, that the excess of resistivity 
after irradiation was ~76mpohmcm, that is about four times higher 
than the Ap, produced by quenching from 600°c; (3) the isothermal 
annealing curves are of the second order (as after quenching from 600°c, 


as will be shown in a forthcoming paper). 
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However, there are other facts which leave some doubts : (1) McReynolds 
et al. do not report the existence of a second stage (it might be that 
sufficient attention was not given to this possibility) ; (2) also the critical 
shear stress is recovered with resistivity and there is no clear mechanism 
for this in terms of clustering. 

Also Blewitt et al. (1957) have reported (fig. 20 of their paper) the 
existence of a stage in the range from — 80 to — 20, after irradiation with 
fast neutrons at 10°K, which they attribute to the annealing out of 
vacancies. 


4.3.2. On vacancies lost during quenching 


As said before, a fraction of vacancies is lost during quenching from 
high temperature, that is the Ap) experimentally observed is lower than 
the Ap, given by eqn. (1). These lost vacancies could have, of course, 
some bearing on the subsequent annealing process; for instance they 
could act as nuclei for the clustering process. Therefore, in view of the 
impossibility of establishing from the actual data the precise role of 
these lost vacancies (if there is one), it is necessary to consider the actual 
results as peculiar of the adopted quenching technique; namely they 
cannot be, a priori, extrapolated to other quenching conditions without 
an experimental confirmation. 


4.3.3. On the annealing of vacancies in gold 


Koehler et al. (1957a) have published recently their results on the 
recovery of vacancies in gold, after quenching from several temperatures. 
It is useful to mention here some essential points of similarity or dis- 
similarity between the results for gold and for aluminium. The principal 
results of Koehler et al. can be synthetized as follows: (i) the migration 
energy decreases from 0-8ev for a quench from 700°c to 0-6ev for a 
quench from 1000°c; (ii) the annealing ‘ half-life ’ time decreases largely 
with increasing quenching temperature; (iii) a decrease of the volume of 
the specimen, proportional to the decrease of resistivity, is found; (iv) 
with quenching temperature an increasing fraction of the quenched-in 
resistivity is not annealed at low temperature (0-10 for a quench from 
1000°c). 

Koehler et al. (1957b) interpret their results in terms of single and 
multiple vacancies, especially to explain the effect of quenching tempera- 
ture on the migration energy and on half-life time. The presence of larger 
clusters seems to be precluded especially by the results of point (iii). 

On the contrary, in the case of aluminium the clustering is needed 
to explain, inter alia, the reduction in the number of jumps and the 
increase in H,, as the quenching temperature is increased from 470 to 
600°C. This last fact is not observed in gold, as shown by point (i). In 
the same way, if we estimate the number of jumps n in gold, by using 
eqn. (2) and the half-life times and migration energies reported by 
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Koehler et al., we find for the highest quenching temperatures n ~ 108, and 
this is too large a value to be connected with clustering. . 

In conclusion, although multiple vacancies appear to exist in each 
metal, the tendency to clustering appears to be much stronger in 
aluminium than in gold. It is to be noted, however, that probably some 
clustering occurs in gold also, as shown by results of point (iv) 


§ 5. CONCLUSIONS 


The principal experimental results can be synthetized as follows: 
(1) notwithstanding the relatively large wire diameter, a large fraction 
of vacancies can be quenched-in ; a formation energy of 0-76 ev is obtained 
from the lower quenching temperatures; (2) isochronal annealing after 
quenching from mean and high temperatures show the presence of two 
stages; after quenching from lower temperatures (400°C) only one stage 
is observed ; (3) the first stage is centred at about room temperature and 
takes place at lower temperatures the higher the quenching temperatures ; 
(4) the second stage is centred at about 160°C; there is some evidence of a 
slight dependence upon the way the first stage has taken place; (5) the 
amount of the second stage is related to the total increase of resistivity 
App, by the expression Apy’=a.Tp,™ where m~3 and a~3-8 x 10-2; 
(6) the migration energy for the defects of the first stage is 0-58 + 0-03 ev 
for quenching from high temperature and somewhat lower (~0-4ev) for 
lower quenching temperatures ; (7) the number of jumps n for annihilation 
in the first stage increases from ~ 104 for quenching from 600°c to ~ 10’ 
for quenching from 550°c, and to ~ 10° for lower quenching temperatures ; 
(8) the activation energy for the second stage is about 1:3+0-2ev; 
(9) the fictitious formation energy of the second stage is ~ 1-2 ev. 

The experimental results have been interpreted, in a qualitative way, 
in terms of clustering of vacancies; namely the first stage observed should 
be related principally to a coagulation of vacancies in clusters of some 
unspecified shape and the second stage to the subsequent dispersion by 
self-diffusion of these clusters. 
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ABSTRACT 
The following conclusions are reached: 


(i) In a medium comprising two identical interpenetrating ion streams. 
four electric space-charge phenomena may occur. The ‘ fast’ and ‘slow ’ 
travelling waves do not concern us; the two effects of interest are non- 
travelling in a system in which the medium as a whole is at rest. 


(u) The first is most probably a set of evanescent waves. These have 
erroneously been interpreted as spatially growing waves and as the basis of 
the electron-wave and similar amplifiers. 

(iii) The second effect is not an oscillation at all but an instability in the 
medium. A single space-charge cloud, however low its density, will, if its 
linear dimension is large enough, continuously increase in density with time. 
The mechanism disclosed is quite simple and does not require, as previously 
suggested, a series of ion clouds (that is, a spatial wave). 

(iv) The effect may occur widely in nature, for example in the solar atmo- 
sphere (leading to patchiness and increased radio emission) in interstellar space 
and in the ionosphere (sporadic —#). It also provides a physical picture of 
the action of the electron-wave tube. 

(v) An earlier criticism of Bailey’s theory of radio waves growing in a 
medium which drifts along a magnetic field is consistent with special relativity 
theory. The consistency depends on a suitable definition of ‘ growth’, 
being increase in amplitude away from the source of the wave. The original 
definition—increase in amplitude in the direction of travel—indicates growth 
of virtually all waves to suitably moving observers. 

(vi) The earlier criticism is also confirmed in discussing a previous Laplace 
transform analysis of the same situation plus a reflecting boundary. A new 
type of wave (the ‘ transverse space-charge wave ’) described in this analysis 
is shown to be spurious and a corresponding theory of solar radio noise un- 
tenable. 


§ 1. INTRODUCTION 


ELxEcrric space-charge effects in interpenetrating ion streams have been 
the subject of numerous investigations. However, the different phenomena 
which may be present do not seem to have been adequately isolated, 
identified and described physically. : 

important in nature, in interstellar space, in the sun’s atmosphere and in 
the ionosphere, as well as in the laboratory. For this reason a better 
understanding is desirable and an attempt is made here to provide this 


for the simple case of two equal density streams. 


+ Communicated by the Author. 


The effects concerned may be 
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The different effects are isolated and it is found that, apart from 
ordinary travelling waves which do not grow or decay (collisions being 
neglected), there are two phenomena of interest. The first is an 
instability which may be described as a “ growing ion cloud’. It has been 
disclosed by previous mathematical analyses and its physical mechanism 
is discussed below. It is thought to be the basis of the two-stream 
laboratory amplifier and to have possible astrophysical significance, 
particularly in the field of radio astronomy. 

The second effect has sometimes been interpreted as a spatially growing 
wave but there seems no evidence for this belief and the phenomenon is 
probably an ordinary evanescent wave as previously suggested 
(Piddington 1956 a). 

It seems appropriate to include here a discussion of * growing ’ radio 
waves in a plasma drifting along a magnetic field. The reason is that 
the phenomenon involved is similar to one of those in the double ion 
stream and the discussion almost identical. A previous conclusion 
(Piddington 1956 b), that the growth was spurious, was based on a steady- 
state analysis. It is confirmed by reference to a transient analysis and 
also shown to be consistent with the principle of (special) relativity. 


§ 2. THE DispeRsION EQUATION OF THE DOUBLE-ION STREAM 


Consider two interpenetrating streams, each of infinite extent and 
each consisting of N protons and JN electrons per cm?. Initially we 
consider an observer at rest relative to the whole plasma, the streams 
having velocities +v. The dispersion equation for disturbances having 
the form exp 7(wt— kx) then has the form (Piddington 1956 a, eqn. (10)) 

Wo oe 
jah eae ee a 
where w.’?=47Ne?/m and ¢« and m are the electron charge (e.s.u.) and 
mass respectively. 

Equation (1) is a biquadratic and so indicates that two pairs of waves 
may propagate in the medium. Each pair comprises two similar waves 
propagating in opposite directions. Some characteristics of these two. 
types of wave may be displayed in a diagram where k is plotted against 
w as in fig. 1 (a) where only positive values of w and k are used and where 
the scales are simplified by putting v=1 and wWy= Vs. 

For values of w > 1 and k>1, that is, for frequencies and wave numbers 
above certain critical values, the waves are ordinary travelling waves of 
the form expi(w,t—k,«) with both w, and k, real. The phase velocities 
of these waves correspond to the two values of w,/k, so that those 
corresponding to curve Al are slow travelling (velocity <v) and those 
corresponding to curve B1 are fast travelling waves. Neither of these 
meres shows either growth or decay and so they do not concern us directly 
nere. 


Nearer the origin O we have two other phenomena, which are the ones 
of interest. 
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(i) The A2 curve, obtained by assuming k real in eqn. (1) 

The values of w for k<1 are imaginary} and it might fe presumed 
that the corresponding waves remained stationary in the medium but 
grew or decayed with time. To settle this point a Laplace analysis of a 


Fig. 1 


Wave Number (k) for curves At Bi A2 
Space Growth Factor (k) for Ba 


Frequency (s) for curves Ai Bi Ba 
Time Growth Factor (lw) for A2 


3T k @ curve) 
tk (Bz curve) 


S | 2 3 
(b) 
Graphical representations of the dispersion equations of two types of 
electron gas. 
(a) The double-stream gas. 
(b) Thermally agitated gas. 
+ The critical wave number k=1 applies when v=1 and w)=/%. In general 
the values of w are imaginary for 0<k< ¥/(2),/v. 
AP 
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transient situation is required. No boundaries of the medium need be 
assumed ; the transient need only comprise an initial weak spatial distortion 
of the medium whose subsequent development is studied. A previous 
analysis of this type is discussed in the following section. The discussion 
is necessary because the results, as related to the dispersion equation, 
appear to have been misinterpreted. 

(ii) The B2 curve, obtained by assuming w real. 

The values of k for w <1 are then imaginary. The B curves of fig. 1 (a) 
may be compared with those of fig. 1 (b) for a gas whose ions have thermal 
motions only. The latter correspond to travelling waves (Bl) and 
evanescent waves (B2). The B2 waves of the double ion stream may 
similarly be reflected waves shown schematically in fig. 2 (a) or, alterna- 
tively, exponentially growing waves as shown in fig. 2(b). This ambiguity 
is discussed in § 4 below; it cannot be resolved by the simple substitution - 
analysis methods described above and has been the cause of some 
confusion. 

Fig. 2 


Amplitude 
Amplitude 
ie 
ey 


(a) (b) 


Spatially exponentially varying standing waves. 


(a) Attenuated (evanescent) waves. 
(b) Growing waves. 


§ 3. THE A2 Curve—Grow1ne Ion Ctoups 


In an unbounded two-stream medium Twiss (1951 a) has assumed an 
arbitrary initial spatial distribution of charge at time t= 0 within a narrow 
region 0<a<d. He found that this ‘ion cloud’ drifted with the plasma 
as a whole and at the same time steadily increased in amplitude. This 
simple and important result was later complicated by the introduction of 
a time modulation which is discussed in §6 below but need not concern 
us as yet. 
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To relate Twiss’ analysis to the steady-state phenomena described 
above consider a single spatial (8-function) distribution of charge carried 
along by the drifting plasma. This distribution may be identified with a 
Fourier integral or, what is more easily visualized, a Fourier series having 
a large number of terms. Each term, or each wave, is defined by a 
particular value of & or a particular wavelength 27/k. Each value of k 
is related by eqn. (1) or graphically by fig. 1 (a) to two values of w. Thus 
each term of the Fourier series corresponds to one or two points on the 
curves of fig. 1 (a). Every term of the Fourier series must correspond to a 
wave uniform throughout space and described by a real wavelength and 
real value of k. They correspond therefore to the curves Al, Bl and A2 
(but not B2) of fig. 1(@). The travelling waves (Al and B1) cause the 
disturbance to spread, as found by Twiss, The A2 effect, with w imaginary, 
causes it to increase in amplitude with time. Hence as far as growth is 
concerned, the analysis is concerned with the A2 effect, having imaginary 
values of w. 

The same phenomenon has also been discussed by Twiss (195la) in a 
steady-state analysis. He determined w for real values of k and stream 
velocities w,, Up. 


. Pie p2(oy gg \2) 121 ae 
om (S54) beef 1+ U;) ae aa ] - (2) 
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Two of these roots are conjugate complex when 
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Twiss referred to the growth corresponding to the imaginary part of w 
as ‘instability ’ in the medium. On transforming to a coordinate system 
moving with the gas (w,+w.=0) and putting w.—w,=2v, eqn. (2) reduces 
to eqn. (1) and Twiss’ ‘instability ’ is identified with the A2 curve of 
fig. 1 (a). 

In relating his Laplace transform and steady-state analyses Twiss 
seems to have made an error. He concluded that the former, together 
with its indication of growth, was a different phenomenon to his 
‘instability ’ which he considered unimportant. He believed that the 
Laplace transform analysis revealed a second phenomenon involving 
growth which he called ‘ amplification ’ and which, in the steady-state 
analysis corresponded to imaginary values of k. This phenomenon he 
considered important being responsible for the operation of two-stream 
amplifiers. 

These conclusions are inconsistent with the present discussion. As 
shown above the growth indicated by the Laplace analysis can only be 
associated with the A2 effect (k real, w imaginary). ‘This effect 1s 
important; a physical explanation is given in § 5 and in § 6 it is shown 
that operation of the two-stream amplifier may be explained by this 
effect. On the other hand there seems no evidence that Twiss' 
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‘amplification ’ effect is important. It corresponds to k imaginary, 
w real, that is to the B2 curve. In the following section it is suggested 
that this is probably an evanescent wave. 

Using steady-state or substitution analyses, the A2 phenomenon (é real, 
w imaginary) has been investigated by Bohm and Gross (1949a) and 
Kahn (1957) who obtain, within certain wavelength bands, a dispersion 
equation of the form w?=— Ak, where A is real and positive. This 
equation suggests the possibility of growth but may not be acceptable as 
full proof (Piddington 1956a) in the absence of a description of the 
mechanism of growth. 


§ 4, Tur B2 Curve—ProBaBLy REFLECTING WAVES 


The roots of the dispersion eqn. (1) corresponding to the B2 curve of 
fig. 1(a) (w real, k imaginary) may be represented either by fig. 2 (a) or 
fig. 2(b). The amplitude varies exponentially in space and, at a given 
point, it varies sinusoidally with time. The steady-state analysis cannot 
distinguish between the two possibilities of decaying and growing waves 
and since no appropriate Laplace transform analysis has been made the 
ambiguity remains. 

However, it has been shown (Piddington 1956a) that previous 
demonstrations of growth of these waves were not valid. Ifthe explanation 
given below of the growth mechanism in double-stream amplifiers is 
acceptable, then there remains neither theoretical nor experimental 
evidence that the B2 waves are growing waves. 

Furthermore, when considered in connection with the double-stream 
amplifier, the waves of fig. 2(b) seem improbable, if not physically 
unrealizable. Between the input electrode and infinity (on the simple 
theory there is no limitation in length of the amplifier), at a given instant 
the perturbation space-charge is all of one sign and increases (spatially) 
exponentially. A half-cycle later it is all of the opposite sign, the change 
being caused by an applied signal which may be vanishingly small. 

On the other hand, waves of the form shown in fig. 2 (a) are well known 
and understood—they are evanescent waves. It seems very probable 
that the double-stream B2 waves are of this type. 


§ 5. THe MECHANISM OF THE GROWING Lon CLOUD 


There does not seem to be a satisfactory physical explanation of the 
instability or growing ion cloud phenomenon represented mathematically 
by certain roots of eqn. (1) and by curve A2 of fig. l(a). Perhaps the 
omission is due to the fact that the effect has generally been regarded 
essentially as a wave phenomenon. Growth was thought to be due to a 
weak initial disturbance, say at «=0, being carried along by one stream 
and growing due to velocity modulation. By interaction, the signal was 
then transferred to the second stream and carried back to x=0 in an 
amplified form (Bohm and Gross 1949 a). 

There seems little justification for assuming that signals are carried 
along by the separate streams; rather do disturbances travel with or 
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propagate through the medium asa whole. Furthermore, it will be shown 
that the instability or growth applies for a single cloud and is simpler than 
_ the above ‘ feedback ’ mechanism. It may be described briefly as follows : 

consider a neutral stream passing through a cloud of negative space- 
charge. On entering the cloud each electron will be slowed down by the 
space-charge field and on leaving, each electron will be accelerated to its 
original velocity. The corresponding effects on the heavy ions are 
neglected because of their much greater mass. The electron velocity in 
the cloud is lower than that outside the cloud and to maintain continuity 
of flow the electron stream density must be greater inside the cloudy. 
Thus the retardation of the electron stream contributes to the initially 
assumed excess electron density. This contribution may be greater than 
the whole of the initially assumed cloud, in which case the cloud density 
increases spontaneously. Two ion streams (or one stream plus a fixed 
grid system) are necessary because a space-charge cloud which tended to 


form in a single stream would be swept along with that stream and could 
not grow. 


Electron Density 


2k ae ey 


The electron density distribution in a hypothetical double electron beam. 


Consider a neutral medium through which are flowing two electron 
streams of uniform density N and velocities along the x-axis of + vp. 
Assume for the moment that a space-charge perturbation of the form 
shown in fig. 3 may exist ; in the region — 7/2k <a <7/2k the total electron 
density is increased from its normal value of 2N to 2N+n where 
n=n,coska. The perturbation is caused by inserting a cloud of electrons 
of density n(x) attached to a rigid but permeable grid. ie function Ns 
must be determined as well as those m) and ,, which are perturbations 
in the stream densities caused by the charged grid. We have, of course 


M=NytNg tg. + - + + 2 + + | (4) 


The presence of the space-charge cloud changes the stream velocities 
from +, to Vp +2, and —v+%. To make the picture clearer physically, 
OE 


+ This situation corresponds to the * travelling pulse ’ solution of Bohm and 
Gross (1949 b). The correspondence does not extend, however, to the mech- 
anism of growth developed here; the mechanism actually proposed by Bohm 
and Gross (1949 a) is mentioned above. 
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shielding grids may be assumed at the cloud boundaries so that the 
external field is everywhere zero and the stream velocities + Up. 
Poisson’s equation gives (e.8.u.) 


— = —4men, cos kx, (=ha <i <t0)is = ee 


which on integrating and noting that by symmetry EH=0 when x=0, 
yields 


Teme sin ke. , 2 see 
Let us assume, subject to verification, that the situation is stable so that 


a/at=0. The equations of force and continuity for the first stream and 
for small perturbations are 


Or et a 
°0 Bar tn i 
CG] 
ag WN Ouch teat) = . . . . . x . (8) 
Combining with (6) and integrating again 
i wg Mo cos kx (—in<kx<4n) (9) 
iss ork? 5) 2 2 . : : . 


where w,?=47Ne?/m. The value of n, is the same, so that 


‘ 2 
Ns Ny COR ha ice ais): shy She OL 
This equation points directly to the physical reason for growth of ion 
clouds. 

When 4/(2)wo/upk<1 so that nz is a positive quantity we have the 
stable situation of a cloud of electrons of density n, attached to a rigid 
frame inserted in the double stream. The slowing down of the streams, 
given by 


ai 2 
a ATTN p€ 


V4= cos kx, 


mug k 
causes an increase in density in each stream so that the resultant cloud 
has density n>n;. The stability results from the mechanical forces 
exerted by the frame in holding the original cloud together; otherwise 
this would presumably dissipate as a set of travelling waves. 

Now increase the dimensions of the cloud until 4/(2)w,/vk=1 and 
ng3=0. ‘The situation is simpler, as the original cloud held by the frame 
is dispensed with. A perturbation of this critical size (the wavelength 
A= /(2)mv9/w9) is self-sustaining and again our assumption that d/ot=0 
is justified. 

On further increasing A a new situation arises; an initial assumed cloud 
has such a powerful effect in slowing and concentrating the streams that 
the resulting cloud is more dense than that assumed (nz negative), the 
effect, being independent of the assumed density of the original cloud. 
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The more dense cloud in turn causes an even more dense cloud and so on 
to continual growth—for disturbances beyond a critical size the medium 
is unstable. The assumption 0/@#=0 no longer holds so that the full 
equation of force and continuity become necessary to determine rate of 
growth; however the simpler equations served their purpose in showing 
when and how growth occurred. The growth criterion k< /(2)a9/V9 is 
the same as that defining curve A2 of fig. 1 (a). 

The above arguments may be modified to show that a cloud of positive 
space-charge (an electron deficiency) also grows, the electron streams 
being attenuated in passage. Growth might also be anticipated when a 
single electron stream moved past stationary positive ions (an electric 
current). 

The effect is essentially an instability rather than an oscillation (where 
an initial distortion is resisted and replaced by one of opposite sign). 
However, as seen below, to a suitably moving observer a travelling spatially 
growing wave may be simulated. 


§ 6. GRow1ne Ion CLoups IN NATURE AND IN THE LABORATORY 


In solar atmosphere clouds of ionized gas are seen to move with velocities 
up to 10? km sec~! or more. Near the borders of these clouds there must 
be double-stream regions where space-charge clouds may grow. The 
minimum extent of an electron cloud which grows spontaneously is given 
by A= 4/(2)mv,/H . In the corona where N=108 cm-? the values of A 
for cloud velocities of 10°, 10®cmsec—! (v,=5x 105, 5x 108 cm sec-*) 
are 0-004, 0-4em. This may be compared with the collision distance of 
electrons, about 108 cm for a gas temperature of 10° °K. It is evident 
that the extent of the double-stream region, limited to the electron collision 
distance, is ample to provide growing electron clouds so that every slight 
irregularity tends to grow into a dense space-charge cloud. 

Possible secondary effects of small ion clouds in the solar atmosphere 
are of interest. They create localized electric fields which accelerate 
electrons passing through or near the clouds. ‘This in turn causes a type 
of thermal radio emission as discussed by Hoyle (1953) and Piddington 
(1954). The latter showed that under ideal conditions, when all the 
electrons exist in separate clouds, the emission may be increased over 
the free—free value by a factor 0-002 (v,/f)2N. v, is the electron (thermal) 
velocity, f the radio frequency and N the electron density. For a gas 
temperature of 106 °K, v, ~7 x 108 cm sec! and for an electron density 
of 108 cm-? the factor at 100 Mc/s is about 10%. It seems unlikely that 
the ideal conditions required would ever be remotely approximated but 
even so, the factor is so large that the effect may be important in contri- 
buting to solar radio emission. 

Solar radio emission also contains components of such high intensity 
that their origin would seem to require the ordered motion of large numbers 
of electrons rather than random (thermal) motions of individual electrons. 
The space-charge clouds seem to provide a possible basis for such an 
effect, provided there are localized irregularities in the gas density and 
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velocity. ‘The emission would result from the ‘collision’ of two clouds 
rather than of a single electron with a cloud. 

The electrostatic energy of an electron cloud is provided at the expense 
of the kinetic energy of the streams which are consequently slowed. 
Khan (1957) has considered this effect in connection with the collision of 
interstellar gas clouds whose relative motion, he concluded, might be 
arrested because of the instability of the system. The theory raises 
interesting possibilities but does not, however, consider the motions of 
the heavy ions which carry most of the gas kinetic energy. 

Yet another effect is revealed when account is taken of collisions between 
particles of the two streams. This provides ions which do not share the 
streaming velocities and so may become permanent parts of the ion clouds. 
A cloud of electrons will attract positive ions and so create a co-existing 
neutral gas cloud. This effect may be important in the solar atmosphere 
where there is evidence of small-scale irregularities in the gas density, 
It may also play a part in creating ionospheric irregularities such as 
sporadic-# for which no satisfactory explanation seems available. 

The growing ion cloud or A2 effect also appears to be the basic mechanism 
of the two-stream amplifier (electron-wave tube and similar devices). 
However, whereas the A2 effect shows growth in time (k real, w imaginary), 
the two-stream amplifier shows spatial growth (k complex, w real). The 
apparent discrepancy lies in the assumed methods of introducing the 
waves into the medium. In the former case a wave of uniform amplitude 
is assumed already present throughout the medium at a given instant 
so that k=k, and w= —iw;, where k, and w; are real and positive. In the 
second case a sinusoidally varying input signal creates a series of identical, 
equally spaced ion clouds in the medium moving with velocity U (the 
two streams have velocities U+v, U—v) past the emitter. These drift 
along, each growing in time at the rate given by exp w. In a time 
interval t, a cloud moves a distance Ut and so shows an apparent spatial 
growth given by exp w,x/U. Thus the appropriate value of k is k, + ia,/U ; 
as before, for an observer moving with the stream w= —iw,. In the two- 
stream amplifier the observer is assumed fixed as the emitter so that a 
Galilean transformation, involving a velocity change —U, is necessary. 
The value of & is unchanged, in this transformation the value of w is 
increased by an amount Uk so that 


heck tw; = j lw; 
k,.+ U 5 WwW tw;+ O(ky+ a) 
Digi. 


These expressions correspond to the type of wave found in the two-stream 
amplifier. 


eee EEE Ee 
} Twiss extended the simple situation described in § 3 above by assuming a 
time modulation, at a new arbitrary frequency Q, of the spatial distribution of 
ee in the narrow region 0<a#<d. Since this region drifts through the gas 
‘he result is a series of equally spaced ion clouds released with increasing time 


delays. The physical growth mechani is 
\ phys ‘chanism is unchanged but t si i 
corresponds with the two-stream amplifier. : 2 a 
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The early theory of double electron-stream amplifiers (see, for example, 
Pierce 1948) has been criticized (Piddington 1956 a) on the grounds, 
briefly, that the same theory would show ‘growth’ of any reflected wave 
in any medium Pierce and Walker (1956) stated in reply that “no 
definitive conclusions can be drawn from it (the steady-state analysis) 
alone’. They added that a Laplace transform analysis is capable of 
providing such conclusions. 

This difference may now be clarified in terms of the phenomena referred 
to here as the A2 and B2 effects. The former (the growing ion cloud) is 
an instability in the medium by which kinetic energy of relative drift 
of the two streams is converted to electric potential energy. It is 
essentially a continuous, non-oscillatory process but as seen above may 
be made to appear as a spatially growing oscillation. The B2 effect was 
the basis of the earlier criticism, being accepted by Twiss (1951 a) and 
others as an ‘amplification’ effect or spatially growing oscillation of the 
type illustrated in fig. 2(b). There appears to be no evidence for this 
conclusion ; the effect is probably an ordinary evanescent wave as previously 
suggested. 


§7. GRowine Rapio WAVES 


As seen above, an evanescent wave may appear to be a spatially growing 
wave when both observer and emitter drift through the propagating 
medium. This applies equally for radio as for electric space-charge 
waves and is the basis of a previous criticism (Piddington 1956 b) of a 
theory of growing radio waves. For completeness the earlier discussion 
must be extended as follows: 


(a) To include the boundary problem, solved by the method of Laplace 
transforms by Twiss (1951 b), who concluded that real wave growth 
could occur. This is done in § 8 below. 

(6) To show that, contrary to a recent assertion (Bailey 1956)}, the 
previous criticism is consistent with Einstein’s principle of (special) 
relativity. This is done in the following paragraphs. 


The dispersion equation for radio waves of the form exp U(wt — ka) in 
a plasma with mass drift along a magnetic field has been given in a simple 


form by Twiss (1951 b). 
k= +ia+U6 rae ee LL) 


where « is a function of w and is real within the wavebands concerned, 
U is the velocity of the gas relative to the observer and f is a (real) constant. 
When U is small, even vanishingly small, one wave has a positive phase 
velocity (in the +x direction) and increases in amplitude in that direction. 


+ Repeated references in this note to ‘views about the transfer of energy ’ 
together with a simple illustration that the transfer of kinetic energy of a par- 
ticle is not invariant, imply that the earlier criticism of Bailey’s theory was 
based on energy considerations. This was not the case: as shown below a 
discussion of relative wave amplitudes was sufficient to show the theory un- 
tenable. Energy did not and need not enter into the arguments. 
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Bailey (1950) interpreted this as a spatially growing wave but an alternative 
explanation is that the wave is the left-hand wave of fig. 2 (a) moving 
to the right and growing to the right but decaying spatially away from its 
source at «=0. The wave moves because the gas drifts relative to 
observer and emitting boundary (Piddington 1956 b). 

The basic equations used in the original theory (Bailey 1950), the theory 
itself and its criticism are all independent of considerations of energy. 
We are concerned with the amplitude of the wave as a function of time and 
space. Further, since only small velocities need be invoked, transforma- 
tions between relatively moving systems are Galilean. Hence wave 
amplitudes and the direction in which the amplitude increases are inde- 
pendent of the observer. 

It is the physical significance of this (invariant) direction of spatial 
growth that is the point at issue. Bailey relates it to the direction of 
positive phase velocity which, in the wavebands concerned, is merely 
the direction of drift of the gas. Using this criterion any evanescent 
wave may show growth. A more satisfactory criterion of growth is that 
the wave amplitude should increase in the direction away from the 
boundary at which the wave is introduced into the medium concerned. 

The reality of wave growth defined in this way is invariant and so 
renders the original criticism consistent with the principle of special 
relativity. 

§ 8. THE TRANSIENT ANALYSIS 


Twiss (1951 b) has described and analysed a laboratory experiment 
corresponding to Bailey’s theory of radio waves growing in a plasma with 
a mass drift along a magnetic field. The experiment is illustrated 
schematically in fig. 4. Media I and III are vacua and Medium II a 
plasma drifting along a magnetic field. A transverse electromagnetic 
(radio) transient is assumed incident on the boundary «=0 at time t=0. 
The boundary x=d is a metallic reflecting surface and the field in the 
region 0<x<d is investigated at times t> 0. 

Before undertaking the Laplace transform analysis Twiss derived the 
appropriate (steady-state) dispersion equation. He modified the corres- 
ponding equation for the two well-known circularly polarized radio waves 
(ordinary ” and ‘extraordinary’) propagated along a magnetic field in 
an ionized gas 


Ch? = w* — w?(1 + Q/w)-1 er es ak) 
where w, and © are the plasma and gyro-frequencies and ¢ the velocity 
of light. The modification was a Galilean transformation to an observer 
moving with the emitter (w assumed real) relative to the gast. The 


{ The dispersion equation is then of the third degree and Twiss has claimed 
that a third type of wave, a “ transverse space-charge wave ’’, is present. 
Before considering wave growth it is well to be sure which wave is under 
discussion. To this end it is shown in Appendix I that the ‘ transver 
charge wave ° has no reality and that the only waves 
‘ ordinary ’ and ‘ extraordinary ’ waves viewed by a ( 


se space - 
present are the original 
slowly) moving observer. 
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appropriate dispersion equation then has the familiar form (neglecting 
electron collisions) of eqn. (11) where Bis real and, within certain wavebands, 
« is also real. One of these waves evidently grows spatially in the +a 
direction and one decays. Twiss used the transient analysis to find if 
the growing wave appeared in Medium II (fig. 4). 

Part of the transient solution revealed waves travelling in the + a 
direction and decaying. Twiss concluded that if another part revealed 
waves travelling in the —x direction, then these must be growing waves. 
Part of the solution (Twiss 1951 b, § 5, line 4) did correspond to waves 
travelling in the —x direction. This was shown by the fact that the 
contribution made by this term to the electric field at any plane « was 
identically zero until a time ¢Sc-(2d—a) where x=d is the reflecting 
boundary of Medium II. The waves concerned must have been reflected 
and so must be travelling in the — 2 direction. Twiss concluded that this 
result indicated the presence of growing waves. 


Fig. 4 
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Boundary conditions for * growing ’ electromagnetic waves. 


The error in this argument results from the fact that the wave shown 
travelling in the —z direction had an infinitely high frequency. The 
term in Twiss’ transient solution with which we are concerned comprises 
a set of radio waves of all frequencies (0<w< _«). Above the critical 
frequency, that is in the frequency range where « of eqn. (11) is imaginary, 
the waves are ordinary travelling waves, not growing or decaying. In 
particular the wave of infinite frequency is not affected by the medium 
at all and so travels with the velocity of light and returns to the plane « in 
time t=c71(2d—x). It was this wave which Twiss used in testing for 
reflected waves and the answer was just as anticipated from simple theory. 
However, this wave is not a growing wave so that the presence of a growing 
wave in Medium II has not been demonstrated. 

According to eqn. (11) and to Bailey’s theory, growth may occur only 
in certain wavebands where « is real so that the waves have the form 
exp tax exp i(wt—UBx). Viewed in a frame moving with the gas 
(U =0) these are standing waves and so cannot be reflected at the boundary 
at x=d and cannot travel against the direction of gas flow. Hence there 
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are no exponentially spatially varying waves in Medium II travelling 
in the —x direction; no growing waves can be present. 
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APPENDIX 1 


Waves in a Drifting Gas 


When ionized gas drifts along a magnetic field the dispersion equation 
for electromagnetic waves is modified from eqn. (12) above by the trans- 
formation w,=—Uk, k,=k. Each of the two equations is now of the 
third degree suggesting that instead of two ‘ordinary’ and two ‘extra- 
ordinary’ waves there are now three in each case. We will consider the 
physical significance of the third wave. 

It might be suggested as obvious that the mere motion of an observer 
along an existing wave (or waves) could not create another wave. Twiss 
explains the third wave as due to the observer’s drift past successive 
layers of gas in which the electrons have different (transverse) motions, so 
creating a wave of velocity ~U. However, all the electron motions are 
integral parts of the radio waves which existed prior to the drift and are 
only modified by the drift. In a linear medium they cannot create a 
fresh wave. It remains to explain the fact that three different values of 
k satisfy the new dispersion equation. 

This is best done by plotting values of k against w as shown in fig. 5 
for either radio wave (the extraordinary wave is chosen) in the absence 
of drift. For any particular value of w, say w,, there are two values of 
k= +k, found by drawing a vertical line to cut the curves. Similarly for 
w, we find +k,. Now introduce a drift and for the two unchanged values 
of & (in a Galilean transformation k,=k) we find two new values of 
frequency given by 


Wig = Wa Uk,, W1p = Wy Uk,. 


In effect the diagram is now drawn in skew axes, the new points corres- 
ponding to @,, ,, having moved to the left of the old as shown. As 
the rate of drift increases slowly the points move further to the left but 
that corresponding to w,, moves faster and eventually overtakes w,, at 
the point marked @,. The skew axis is now parallel with the dashed line 
which passes through three points marked &,, k, and k,, corresponding 
with three waves having these propagation constants. 

This does not mean that a new type of wave has been introduced but 
that three extraordinary radio waves, having different wavelengths, all 
appear, to a moving observer, to have the same frequency w,. The 
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reason is that, because of their different wavelengths, they have different 
Doppler shifts, thus compensating for the differences in frequency noted 
by a stationary observer. Figure 5 shows only real values of w and k 
but the same arguments hold when one or other is imaginary. 


A graphical representation of the dispersion equation of transverse (radio) waves 
propagated along a steady magnetic field. 


If this conclusion is correct it removes an inconsistency with a previous 
investigation (Piddington 1955) in which it was concluded that in a magneto- 
ionic medium only four different waves are possible. Apparent additional 
waves evident to moving observers appear to be spurious. The present 
conclusion would also invalidate Twiss’ (1952) theory of sunspot radio 


noise. 
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ABSTRACT 


The energy spectra of the charged particles emitted from aluminium on 
bombardment with 14 Mev neutrons have been measured for the angular 
intervals 16—48° and 132-164°. Groups of protons corresponding to known 
levels in 2’?Mg were found. The ratio of the differential cross section for the 
forward angular interval to that for the corresponding backward interval 
was found to be 2:6. 


§ 1. INTRODUCTION 


Ly previous studies of the reaction ?’Al(n, p) at 14mev (Allan 1957, Colli 
and Facchini 1957, Brown et al. 1957, and Haling et al. 1957) peaks in the 
energy spectra of the protons emitted have been observed and a 
preference for proton emission in the forward direction has been noted. 
However in none of these studies were the statistics and resolution 
sufficient to make a detailed comparison of the spectrum at forward 
angles with that at backward angles. 

This paper is a report of such a comparison in which photographic 
emulsions were used to detect the protons. 


§ 2. EXPERIMENTAL PROCEDURE 
2.1. Apparatus 

A schematic diagram of the apparatus is given in fig. 1. The 150kev 
deuteron beam of the Glasgow H.T. set was incident on a tritium— 
zirconium target producing neutrons by the T(d,n)4He reaction. At 90° 
to the deuteron beam and in the horizontal plane, a wedge-shaped iron 
collimator, 10in. in length and 1-5in. in height with a tapered rectangular 
hole through its centre, produced an approximately uniform beam of 
neutrons with an energy of 14:1+0-1mev. At the centre of the brass 
exposure box the neutron beam was 1-5cm in height and 3cm in width. 
The box, 7in. x 3in. x 1-5in., was lined with gold, 400 in thickness, to 
shield the photographic emulsions from external sources of protons (gold 
has a very small (n,p) cross section). An aluminium foil 14mg/cm? in 
thickness was stretched across the centre of the box with its plane 
perpendicular to the neutron beam. Ilford C2 nuclear research emulsions, 


+ Communicated by Professor P. I. Dee, F.R.S. 


On the Charged Particles Emitted from Aluminium 1257 


Fig. 1 
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400 pin thickness, were placed below the beam on either side of the target 
foil with the emulsion surface horizontal. The neutron beam on emerging 
from the exposure box passed through a nuclear emulsion inclined at 30° 
to the beam and then to a plastic scintillator 2 metres away which was 
used to monitor the exposure. 


2.2 Examination of the Emulsions and Subtraction of Background 


The emulsions were scanned using x50 oil immersion objectives and. 
x 10 eyepieces for tracks of charged particles which entered the surface 
of the emulsion and appeared to originate at the aluminium foil. 

The areas scanned in the forward and backward directions were 
symmetrical about the target foil and subtended angles at the foil of 
16-48° and 132-164° respectively to the direction of the incident neutrons. 

The range, angle of dip and the projected angle of each track accepted 
was measured. Assuming that the tracks were produced by protons, the 
‘proton ’ energies were obtained using the range energy tables of Rotblat 
(1951). Some of these tracks will be due to deuterons and «-particles, the 
latter will have little effect on the spectra because an a-particle with the 
maximum energy possible would correspond to a 3 Mev ‘ proton ’ in fig. 2. 


Table 1 
Angle 16°-48° | 132°-164° 
Area scanned 0-513 cm? | 0-441 em? 
Tracks accepted 2123 741 
reese |. te 
pai from Al 1816 623 


The background for the forward direction was determined by scanning 
the corresponding backward area for tracks whose projected angles and 
dips fell within the ranges of acceptance used in the forward direction. 
The background for the backward direction was determined using a 
similar procedure. This method is probably valid since the areas scanned 
are both close to and symmetrical about the centre of the exposure box 
and well away from the ends where most of the local background is 
produced. The numbers of tracks found and the background subtraction 
are given in table 1. 


2.3. Correction for Energy Loss 
Protons originating within the aluminium foil lose energy before 
emerging from the foil. Allowance was made for this loss by adding to 


the observed energy of each ‘ proton’ the energy it would have lost in 
traversing half the thickness of the foil, taking account of the angle of 
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traversal. Thus protons having ranges in the emulsion corresponding 
to energies of 2 Mev, 6 Mev and 10 Mey were considered to have had initial 
energies of 2-75 + 0-75 Mev, 6-4+0-4mev and 10-2 + 0-2 Mev respectively. 


2.4. Determination of the Neutron Flux 


The extent of the beam and the flux of neutrons in it were determined 
from the emulsion exposed directly to the beam. At various positions 
across this emulsion, a known volume of emulsion was scanned for 
protons resulting from the collision of the incident neutrons with the 
hydrogen present in the emulsion. These recoil protons produce the 
majority of the tracks observed within the emulsion and are easily 
identified by applying the relationship, 


Ep = (14-1 + 0-1) cos?, 6 
where Hp is the energy of the recoil proton and @ is the spatial angle 
between the direction of the recoil proton and the incident neutron. 


§ 3. RESULTS 


The energy spectra of the ‘ protons ’, in the centre of mass system, for 
the two angular intervals are given in fig. 2. There are peaks in the 
spectra for both forward and backward angles. The solid line arrows 
drawn on fig. 2 indicate the position of known energy levels in the nucleus 
"Me, the residual nucleus of the (n,p) reaction, measured using the 
reaction **Mo(d,p)?’Mg, e.g. Hinds et al. (1958). Reading from left to 
right, they are the ground state, 1-0, 1-66, 3-50, 3-56, 3-75, 4:13 and 
4:75 Mev levels. The ground state (n,p) transition is not observed but 
the high energy peaks in both the energy spectra are in good agreement 
with the positions of the higher levels. 

Some of the particles classed as protons will in fact be deuterons, the 
maximum possible energy for a deuteron would correspond to a ‘ proton ’ 
energy of 5:-3mey. Thus the peak at 5-6mev consists mainly of protons 
suggesting energy levels of ?’Mg in this region. At lower energies a peak 
occurs at about 4Mev, an energy which might be expected to lie in the 
continuum of levels for 77Mg. 

The group of lowest energy ‘ protons ’ of the spectra occur at different 
‘proton energies ’ for the two angular intervals. This effect also appears 
in the results of Haling et al. The broken line arrows shown on the energy 
spectra indicate the ‘ proton energies ’ corresponding to deuterons from 
(n,d) transitions to the ground, 1-8 and 3-0 Mey levels of *6Me (Endt and 
Kluyver 1954). It should be noted that the positions indicated do not 
allow for the assumption that the tracks were produced by protons made 
in correcting for energy loss in the target foil and in determining the 
centre of mass ‘ proton energy’ from the observed value. if the low 
energy ‘ protons’ were in fact deuterons these effects each displace the 
low energy end of the spectrum by about 0-3 Mev. Thus at the forward 
angles the spectrum would be displaced about 0-6 Mev towards lower 
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energies while at the backward angles the displacements would cancel out. 
The two low energy peaks in the spectra are thus consistent with the 
reaction 27Al(n,d)2°Mg, indicating that a large proportion of these 
particles are deuterons. 

The ratio of the cross section averaged over the angular interval 16-48° 
to that averaged over the interval 132-164° for the different groups 


observed in the energy spectra is given in table 2. 


Table 2 


ee 


‘ Proton ’ energy 9.0-3:8 | 3-8-5-0 | 5:0-6-4 | 64-86 | >86 All 


in Mev energies 
016°-48° 2.8 2.8 2.1 2.1 3-5 2-63 
9 132°-164° 


In obtaining the values given for the two low energy groups, the 
majority of the particles in the groups were assumed to be deuterons and 
accordingly allowance was made for the displacement of the spectrum for 
deuterons at forward angles as described above. 

All the energy groups show a marked preference for emission at the 
forward angles. The ratio of the differential cross section for all particles 
emitted in the forward angular interval to that for the backward interval 
is 2-63 in good agreement with the value of 2-69 obtained by Haling e¢ al. 
for the ratio of the differential cross section at 30° + 15° to that at 150° + 15°. 
However the absolute value obtained for the average differential cross 
section over the interval 16° to 48° is 19-5 +3 mb/steradian while that of 
Haling et al. is 9-8 mb/steradian at 30° + 15°. 

The asymmetry of the angular distributions of the high energy proton 
groups is probably due to the occurrence of direct interaction effects. 

The preference for emission at forward angles of the low energy groups 
is another indication that they are deuterons. If they were protons from 
the (n,np) process as observed in heavier elements (March and Morton 
1958), their angular distribution would be expected to be symmetrical 
about 90°, whereas a pick-up reaction of the type (n, d) would be expected 
to lead to the emission of deuterons at forward angles. 
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ABSTRACT 


On the faces of a type I diamond octahedron of particularly good quality 
there have been observed clusters of trigons which appear to be associated 
with a centre of strain in the interior of the stone. The trigon patterns show 
features which can only have been the result of surface dissolution. These 
observations are discussed, and it is concluded that most trigons are etch 
pits; there is evidence that these particular trigons are centred on dislocations. 

Possible explanations are considered for the fact that the usual laboratory 
etchants develop pits oppositely oriented to trigons. 


§ 1. LyrrRopUCcTION 


THE origin of the triangular depressions, commonly referred to as 
*‘ trigons ’’, which are seen on all or practically all natural octahedral faces 
of diamonds, has long been a subject for dispute. The observation of 
these features has been notably advanced by Tolansky, who has given a 
fairly comprehensive description of them in published papers and a 
recent book. The two competing theories for their origin are (1) that they 
are etch pits or (2) that they are regions which have been surrounded by 
growth layers, but not filled in. Tolansky has argued strongly, but not, 
as he claims, conclusively, for the latter view. He has two lines of 
evidence : one that he found a case in which the base of a flat-bottomed 
trigon was at the same level as a neighbouring flat area, with an intervening 
plateau ; the other that laboratory etching experiments do not reproduce 
the naturally occurring patterns. Chance cannot be excluded as the 
cause of the first of these : or indeed, it could be the result of a growth 
process without implying that this is responsible for all ‘ trigon’ 
markings. They vary considerably among themselves, and the general 
family resemblance may to a large extent be simply attributed to the 
trigonal symmetry natural on a (111) face of diamond. The laboratory 
etching experiments provide negative evidence which can never be made 
complete. The case from these experiments is considerably weakened by 
the results of Omar and Kenawi (1957) who have etched diamonds (by 
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heating in low pressure air on a molybdenum strip) to produce patterns 
which have all the appearance of natural trigon markings, with one 
exception. The exception is that whereas the sloping sides of the natural 
pits on (111) are of the form (1,1, 1—8), with its permutations, making a 
bounding triangle in opposite orientation to the triangular face of the 
octahedron, those produced by Omar and Kenawi, in common with 
those of most other laboratory etching experiments, have the form 
(1,1,1+6) making a triangle in the same orientation as the octahedral 
face. We may refer to these respectively as triangles in negative and 
positive orientation, and note that no other orientations of triangles are to 
be expected from the crystal symmetry : pits with any other bounding 
faces of a single form would have hexagonal outlines. 

The observations to be described here were made on a good quality 
* blue-white * octahedral diamond of type I. On each face of this stone 
there is a fairly well-defined group of trigons, the approximate centre of 
which is in a line, perpendicular to the face, with the corresponding centre 
on the opposite octahedral face. The four lines intersect at a point within 
the stone (not its geometrical centre) which was the centre of a birefringent 
pattern observed in polarized light. Two such opposed faces, with the 
accompanying birefringent pattern, are shown in fig. 1 (a)—(c), Pl. 82. 

A more detailed study of these trigon groups was made. For this 
purpose some method of contrast enhancement was necessary, and a very 
simple adaptation of a Vickers projection microscope was found to give 
good results. The 45° reflecting plate in this instrument (as used for 
metallography) is a half silvered aperture in a fully silvered plate, 
mounted on a movable slide. This may be partially withdrawn so that 
the fully silvered edge of the plate reflects the illuminating beam and at 
the same time obstructs part of the light reflected from the specimen. It 
has been pointed out by Zernike (1946) and demonstrated by Pratt and 
Lomer (1952) that such a stop in or near the back focal plane of the 
microscope gives rise to a ‘schlieren’ effect which greatly accentuates 
apparent differences in surface slope. The best results were obtained in 
our case with the 16mm objective, probably because the reflecting plate 
coincides most nearly with its back focal plane, while the loss of resolving 
power due to diminution of the objective aperture is not (as with lenses of 
higher power) particularly serious. A further gain in contrast can be 
obtained by inserting into the illuminating beam, slightly off centre, the 
circular stop provided for dark ground illumination. 

The advantage gained with this arrangement can be seen by comparing 
figs. 2 and 3, Pl. 83, of the same field on the unsilvered diamond surface. 
Figure 2 is taken in normally incident light with the conventional metal- 
lurgical microscope; fig. 3 is a photograph by the schlieren method. 
Interferometric measurements have shown that the sides of the fainter 
pits to be seen in fig. 3 makes angles of the order of 10’ of are with the main 
surface. Since this photograph was not taken with the stop in the position 
of maximum contrast, and these pits are not at the limit of visibility, it 
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does not seem unreasonable to suggest that with a silvered surface the 
angular sensitivity of this system could be increased to a minute or so of 
arc. 

The trigons seen may be divided into two classes, ‘ strong ’ and ‘ weak ’, 
according to the inclination of their surfaces. All in the strong class, and 
almost all in the weak class, are of complete, i.e. non-truncated, inverse 
pyramidal form. In both classes a feature repeatedly observed is the 
precise contact of a corner of one triangular outline of a pit with a side of 
another. We consider this a very significant diagnostic feature for the 
origin of these pits. Wherever such contacts occur there is an associated 
depressed area, outside the individual pits, on either side of the contact. 
Such depressed areas have gently sloping floors, with orientations which 
may be approximately described as permutations of (1+ ¢,1—¢«,1—e’), 
surrounded by steeper walls with the same (1,1,1—6) slopes as the 
individual pits : each such outer wall is an extension of a face of an 
individual pit within it. This description applies to the greater part of 
the boundary of a depressed area, but there are portions of the boundary 
of a large depressed area where the step is relatively shallow and curves 
away from the (110) type of direction, and rarer places where it becomes 
straight again with a (1,1,1+6) slope. Depressed areas may contain 
within them further depressed areas, whose steeper boundaries are again 
extensions of the faces of contained pits. The overall depth of all pits of 
the strong class on the one hand, or of all the weak class on the other, is 
about the same when reckoned from the outer plateau, that is, including 
the depth of the depressed areas in which they lie. 


§ 2. THe ForMATION OF THE PATTERNS 


The process producing the pits, whether it be dissolution within, or 
growth without, must by its continuation enlarge them, and so, in time, 
cause neighbouring independent pits to impinge. The growth theory 
must, however, be subdivided into two versions : (i) that the trigon 
patterns are formed only in a final stage of growth; and (ii) that they are 
present throughout the growth history of a diamond. Tolansky’s 
discussion implies version (ii). For this version it is necessary that 
processes by which the pits are filled in, and are therefore conceivably 
capable of shrinking, also operate. On the assumption of etching, the 
outline of a pit may shrink by dissolution of the surrounding area. 
Consider therefore two adjacent pits as in fig. 4. They may approach 
the condition of contact either by enlargement from non-contact as in 
fig. 4 (a) or by shrinkage from overlapping as in 4 (c). The maintenance 
of contact indicated by our photographs, figs. 1 and 2, implies that as 
soon as this condition is passed, a new process operates so as to restore it. 
Passage from contact to fig. 4 (a) introduces no novel crystallographic 
feature into the situation: but passage from contact to fig. 4 (b) brings 
Into existence two salient edges of form (1, 1, 2+ 5) where the two pits 
meet, no salient edges of this kind being present before. We conclude 
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then that the development was in the direction contact to 4 (6). On the 
etching theory the interpretation is straightforward. It is narra to 
suppose that the presence of these salients, being of other form than the 
plainly resistant {110 edges, leads to accelerated dissolution, and produces 
the observed pattern in the manner detailed below. There does not 
however, appear to be any means of interpreting the patterns from a 
growth theory. There is no basis in theory or analogous observation to 
suggest that the salient edges would accelerate growth. If, nevertheless, 


() 


When two pits overlap, salient edges at A and A’ provide points of easy dis- 
solution. Dissolution of the neighbouring surface proceeds until the 
position of exact contact is restored. A depressed area is formed en- 
closed by steps which are extensions of the original trigon faces. 


they did that, so restoring contact by partially filling in the pits, we should 
be left without an interpretation of the adjacent depressed area which 
accompanies contacting pits. To explain the whole figure on the growth 
theory we require to understand how a change of situation at the point of 
contact causes growth layers to hold back at the remote boundaries of 
the depressed area. There does not seem to be any way in which this 
can be explained. Since on the other hand etching will explain the 
figures we conclude that etching produced them. 
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Along the (110) directions of diamond there are zig-zag chains of 
carbon atoms like those in the long-chain paraffins. Such chains parallel 
to any chosen (110) direction contain all the atoms of the crystal. The 
macroscopic habit of the crystal (always involving faces of the (110) 
zones, i.e. predominantly {111}, but also {001} and approximate {110} 
faces which are in fact more nearly (110)-zone cylinders) together with 
the trigonal surface markings indicate that both in growth and in etching 
the extension or contraction of these (110) chains goes easier than the 
initiation of new ones or the breach of existing ones. We may reach a 
closer understanding of this by considering figs. 5 and 6. Figure 5 (a) 
shows in (111) projection the surface of a diamond which departs from (111) 


Fig. 5 


(a) Negative [110] step in (111): the (>) Kinked negative step as in a 
monomolecular step in a (1+e, l—e, 1—8) vicinal face 
(1, 1, 1—8) vicinal face. 


by having one monomolecular step along [110] of the kind present in 
(1, 1, 1—68) vicinal faces. We shall call this a negative [110] step in the 
(111) face. All atoms of the peripheral zig-zag at this step are triply 
bonded to the crystal. An additional atom added anywhere is only 
singly bonded to the crystal. Figure 4 (b) shows this step with one 
kink in it, such as would be present in a (1+.e, l—e, 1—6) vicinal face 
The terminal atom at the kink is doubly bonded to the crystal : its removal 
leaves another doubly bonded atom, the removal of which leaves a kink 
of the same form as before. The same applies to the addition of atoms 
The diamond crystal grows when atoms attach more frequently to such 
doubly bonded positions than they leave them, and dissolve when the 

leave more frequently than they attach. Growth or dissolution at aed 


a step may be analysed into the production of i 
; D such kinks 4 
tion along the step. sand their prop age 
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Figures 6 (a) and 6 (6) show alternate configurations for a step along 
[110] of the kind present at a (1, 1, 1+6) vicinal face. We call this a 
positive [110] step in (111). In either of these configurations the outer- 
most row of atoms is doubly bonded to the crystal, and this bonding is 


Fig. 6 


(c) (d) 
(a) Positive [110] step in (111): the monomolecular step characteristic of 
(1, 1, 1+) vicinal faces. _ 
(6) Alternate configuration of the positive [110] step in (111). eee 
(c) A kink in the positive [110] step, as in (I+, 1—<; 1438) ee aces. 
(d) Alternate configuration of the kinked positive [110] step in (EEL 


not changed when neighbouring atoms of the same row are removed. 
Removal of the outermost atom row from either of these configurations 
produces the other. Such steps can therefore be attacked piecemeal. 
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Their dissolution does not depend on kink propagation following a relatively 
difficult breach of the chain. The same applies in reverse for the addition 
of atoms in the growth process. The kink in this step, figs. 6 (c) and 6 (d), 
in fact provides one ‘protected’ atom, triply bound into the crystal, but 
a single ‘protected’ atom is of no great importance for the kinetics of 
dissolution or growth. The kink in a positive [110] step may be regarded 
as a short length of negative [101] or [011] step. Steps in any other 
orientation in the (111) face may be regarded as kinked positive or 
negative (110) steps. For example, there are two alternate configura- 
tions of a [112] step, in each of which there is an outermost row of atoms 
doubly bonded to the crystal, so that this type of step is also suitable for 
piecemeal attack, or growth. 

This discussion leaves us in need of an explanation for the appearance 
of positive vicinals, (1, 1,1+6), in diamonds etched with potassium 
nitrate, or with oxygen, in the laboratory, and to a minor extent on nautrally 
etched diamonds as in fig. 3. It is suggested as a tentative hypothesis 
that in the presence of oxygen the positive (110) step becomes stabilized 
as —C—O-C-—O- zig-zag chain, such as is present in the polyoxymethylene 
polymers of aldehydes, as shown in fig. 7. This suggests that to imitate 


Positive step as in fig. 6 (b), with an additional line of oxygen atoms. 


the natural patterns one should attempt to etch diamond with reagents 
furnishing only, or predominantly, monovalent atoms and radicles, such 
as halogens or hydrogen. Neither of these seems likely to have been the 
actual natural etchant, and etching by the water-gas reaction, which 


would be expected to attach predominantly hydrogen atoms and hydroxyl] 
radicles to the diamond surface, is a plausible suggestion. 
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§ 3. TRIGoNS as INDICATORS OF DISLOCATIONS 


The pyramidal form of most of these pits suggests that they may be 
centred on dislocations. This suggestion is supported by their location, 
in clusters on each face, around the foot of the perpendicular from the 
interior strain-nodulus: it is reasonable to suppose that the crystal grew 
around this centre, and that dislocations emanating from it would extend 
approximately along the normals to each face. The fact that there are 
essentially two populations of pits, deep and shallow, as may be seen in 
fig. 3 lends further support to this suggestion. Elementary dislocations in 
diamond would have Burgers vectors of the form ${110). 110) directions, 
considered in relationship to a (111) face, fall into two classes, half of them 
being parallel to the face, and the other half making an angle of 54° 44’ 
with it. The latter, having partially screw character with respect to the 
face, would be expected to provide points of particularly easy attack, 
and so give rise to ‘deep’ pits. The former would make points of easier 
attack than the face in general, because of concentrated strain, but would 
probably not facilitate attack so much as the screw dislocations: they are 
more likely to be responsible for the shallow pits. This particular stone 
may be an unusually fortunate example. The relative rates of attack 
at the two kinds of dislocation end would be expected to vary widely 
with conditions in the medium: so also would the ratio of rate of attack 
at the dislocation itself to rate of lateral growth of the pit, which determines 
the steepness of its sides. In some circumstances there might be no 
attack except on the screw dislocations, and it is quite likely that only in 
fortunate cases do we have conditions of etching which produce pits on 
both classes of dislocations, while giving a good discrimination between 
them. 


APPENDIX 


This Appendix gives the results of a sample survey made by one of us 
(E. M. W.) of the slopes of the faces of pyramidal trigons on three diamonds. 
The first (classed as ‘ blue-white’) is the stone with which the main text 
of the paper is concerned. The second is classed as ‘fine white’ and the 
third as ‘glassy’. The figures for slope have been derived from measure- 
ments of lateral dimensions from photomicrographs, and of depth obtained 
by using either multiple beam Fizeau fringes or fringes of equal chromatic 
order (Tolansky 1948). 

The reference surface with respect to which depth is measured is always 
the immediately surrounding flat, or nearly flat, area, which may itself 
be a depressed area with respect to the main surface. On this account 
slope has been taken as the preferable basis of comparison rather than depth. 
In general, however, ‘steep’ and ‘deep’ correspond ; that is to say the 
lateral dimensions of pits on the diamonds studied do not vary very 
widely, nor does this variation appear to be correlated with depth: shallow 
trigons are on the average as wide as deep ones. 
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Not all diamonds are suitable for study, and some selection of specimens 
is necessary. Because of optical requirements the octahedron surfaces 
must be smooth, so only good quality stones can be used. The density 
of trigons on a face must not be too great, otherwise it becomes difficult 
to identify individual trigons in the focal plane of the spectrograph. On 
the other hand there must be adequate numbers of deep and shallow 
trigons to permit a numerical survey, and this condition eliminated some 
of the best quality stones which had no deep trigons and only a few shallow 
ones. The surface density of trigons on the stones finally chosen was of 
the order of 10? to 104 per cm?, which was found to be a satisfactory 
compromise between the requirements. 

In the course of the investigation some trigons with flat bases were 
observed. Although the number of flat-based trigons varies from stone 
to stone, they are in general rare. They have been disregarded in this 
survey. 

The histograms, figs. Al, A2, A3 show the distribution of slopes 
for the three diamonds specified at the beginning of this Appendix. Each 
shows a bimodal distribution. Mean depths and slopes for the two groups 
in each of the three stones are presented in table Al. 


Table Al. Mean Depths and Slopes of Pyramidal Trigons 


Deep trigons Shallow trigons 
Depth | Slope | Number | Depth | Slope | Number 
Stone ; : : 1 : : 
A minutes | measured A minutes | measured 
Blue white 9200 90 18 950 13 27 
Fine white 5570 140 6 500 14 16 
Glassy 3315 49 27 650 Hi 12 


The variation in steepness found (within either the deep or the shallow 
group) would fairly certainly have been less if a survey had been limited 
to a restricted area on a stone, instead of pooling the results for a variety 
of places on several faces of any one stone, as was in fact done. That 
this is the case is shown by the results (table A2) for the five measured pits 


on the ‘blue-white’ stone (all of them deep ones) which fall in the area 
of Plate. 


Table A2. Blue-White Diamond 


Trigon Depth Slope Reference 
number (A) (min) surface 
1 12170 Te main surface 
2 8670 92 depressed area 
3 7850 93 depressed area 
4 11180 109 main surface 
5 7690 107 depressed area 


“4 


Number 


Number 


Number 
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Fig. Al 


100 
Slope (minutes) 


120 140 160 


Histogram of trigon slopes: ‘ blue-white ’ diamond. 


Fig. A2 
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Histogram of trigon slopes: ‘ fine white ’ diamond. 


Fig. A3 
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: - oO ’ om 1's 
Histogram of trigon slopes: glassy ’ diamond. 
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All of these have somewhat greater slope than the average of the deep 
pits on the stone, and the variance among them is only about half as 
great as that for all of the deep pits. 

It will be noticed that the larger variation in depth is mainly due to the 
different reference surfaces with respect to which measurement is made. 
This likewise affects the lateral dimensions of the pits. Variation in 
lateral dimensions is to be expected on this account, and also because of 
variation of conditions over the surface of the stone. Thus, non-uni- 
formity of lateral dimensions does not necessarily indicate non-uniformity 
of environment during the establishment of the surface topography, but 
uniformity of lateral dimensions is an indication of uniformity of environ- 
ment. Consistently with this, the bimodal resolution in slopes is least 
clear in the histogram for the third stone, on which trigons had a relatively 
wide variation in lateral dimensions (about 3:1), is of intermediate 
sharpness for the first, with variation of about 1-5: 1 in the lateral dimen- 
sions, and is very sharp for the second, for which the lateral dimensions of 
trigons was remarkably uniform (about 1-1:1 for the deep pits, 2:1 for 
the shallow ones). 

No statistical significance attaches to the relative numbers of deep 
and shallow trigons in this survey, which are influenced by the varying 
requirements of measuring technique. 

To sum up, it appears that at least for these three stones it is meaningful 
to classify trigons into two groups, deep and shallow, differing in the slope 
of their pyramidal faces by a factor of 5 to 10, despite considerable variation 
within groups, and wider variation from stone to stone. Other stones 
occur, with a relatively wide variation in trigon sizes, for which no such 
clear resolution into two groups is possible. 
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ABSTRACT 


Some experiments on artificial etching of diamonds are described. 
Etching in steam and chlorine at temperatures of around 800°c produces etch 
pits in ‘ positive ’ orientation (i.e. opposed to that of natural trigons). Wet 
hydrogen in the same temperature range and above gives inconclusive results. 
Etching in this medium may be greatly accelerated by the presence of some 
catalyst. 

Diamonds heated in fused kimberlite at about 1450°c develop sharp-edged 
triangular pits in ‘ negative ’ trigon orientation. The bearing of this obser- 
vation on the natural history of diamond is briefly discussed. 


§ 1. INTRODUCTION 


In another paper Frank et al. (1958) have described surface features 
on diamond which strongly suggest that the shallow triangular pits 
traditionally called ‘trigons’ can be formed by etching. It is natural to 
ask whether pits can be artificially etched which resemble those produced 
naturally; Omar and Kenawi (1957) have shown that pits developed by 
heating diamonds in low pressure air reproduce all the characteristics of 
trigons except that they are oppositely oriented. This note describes 
some experiments on diamond etching, one of which resulted in triangular 
pits in the same orientation as trigons. 


§ 2. ETcHING IN GASES 


Frank e¢ al. suggested that a non-oxidizing atmosphere of monovalent 
radicals such as the water-gas reaction might be responsible for natural 
etching. As a first attempt to imitate the process, therefore, a small 
octahedral diamond (whose natural surface features had been recorded 
by photography and replicating in Perspex) was heated in a quartz tube 
through which steam was continuously passed, to a temperature of 750°C 
for one hour. This resulted in a surface covered with flat-bottomed 
triangular pits, mainly with slightly rounded corners, about 2 in size 
and in ‘positive’ orientation, i.e. parallel to the edges of the octahedral 
face and opposite to the trigon orientation. It was thought possible that 
some oxygen had been dissolved in the water from which the steam was 
generated, and a deoxidizer was added before a second heating to 820°C, 
again for one hour. This treatment etched the diamond more deeply, 
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but again with positively oriented pits, and also delineated many ring 
cracks and similar evidences of mechanical damage. 

Other possible sources of monovalent ions are hydrogen or the halogens. 
An attempt was made to re-etch the diamond on which the first experiment 
had been made, with wet hydrogen at 820°c for an hour. This did not 
alter the appearance of etch pits previously formed, but it was noticeable 
that some very deep cavities were scored in the surface. A further attempt 
was made at 900°c, but the furnace burnt out after a few minutes. After 
all this treatment the stone appeared opaque and dark grey, suggesting 
the presence of a surface layer of graphite. 

Another diamond was etched in wet hydrogen at 850°c for fifteen 
minutes (again terminated by the burning-out of the furnace). Most of 
the stones used in these experiments were cleaned in organic solvents and 
nitric acid before heating; this one, however, was only washed in xylol. 
Comparatively few pits were formed, most of which were so rounded that 
no crystallographic orientation was apparent. Through a low power 
stereoscopic microscope another phenomenon could be seen; a number of 
winding tunnels penetrating deeply into the interior of the stone (fig. 1, 
Pl. 84). This remarkably localized attack, originating mainly at points 
on the edge of the crystal, suggests some kind of catalytic action. It is 
known that the formation of methane by the direct union of carbon and 
hydrogen is catalyzed by hydrated iron oxide, particles of which may have 
been transferred from the steel tweezers used for manipulating the 
specimen. 

Further attempts to etch a diamond in wet hydrogen at 1100°c and 
1300°c again resulted in the development on the surface of rounded pits 
bearing no obvious relation to the crystal symmetry. 

There remained the possibility that one of the halogens might be a 
suitable etching reagent for diamond. As a supply of pure liquid chlorine 
was available, this element was tried first. A tube of combustion glass 
containing a diamond was evacuated to backing pump pressure (about 
10 mm of mercury) filled with chlorine to about half atmospheric 
pressure, and sealed. The tube was then heated for 17 hours at 800°c in 
a muffle furnace. 

After this treatment the surface of the stone showed sharp-edged 
triangular pits which in many ways resemble natural trigons (fig. 2 (a), 
Pl. 85). Both flat-bottomed and pyramidal pits can be seen; but they 
are smaller than most trigons, ranging from the limit of resolution of the 
optical microscope up to 1-2. As in other experiments, they are in 
positive orientation with respect to the crystal face. The etching also. 
revealed a number of curious features which seem to have no crystallo- 
graphic habit; the roughly rectangular markings in fig. 2(b) and (c), 
Pl. 85. It is worth noticing a similar feature seen on the surface of the 
unetched blue white diamond examined by Frank et al., an electron 
micrograph of which is shown (fig. 3, Pl. 86). It seems to be a series of 
small islands surrounded by a ‘moat’. The replica surface has been 
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tilted to increase contrast in the photograph, so that the true angular 
relationships are distorted, but it also has no obvious relation to the 
symmetry. Such markings may represent etching of small inclusions 
in the diamond surface. 

The foregoing experiments on etching of diamonds in gases are therefore 
generally inconclusive. The production of pits in positive orientation 
in water vapour and chlorine might in each case be ascribed to the undoubted 
presence of a trace of oxygen, while etching in a reducing atmosphere of 
wet hydrogen seems to have no well defined habit. It was therefore 
decided to try to dissolve a diamond in conditions resembling its natural 
environment. 


§ 3. Ercuine In BLuE Grounp 


It was observed by Luzi (1892) that a diamond dissolves at a temperature 
of 1770°} in fused kimberlite or blue ground, the altered igneous rock in 
which diamonds are found. After heating for half an hour, very deep 
corrosion figures of an oval or half-round shape were found, but as the 
diamonds were blackened, no doubt by graphite, it seems uncertain how 
far the effect was the result of graphitization followed by attack of the 
graphite. Luzi also noticed that his specimens were frequently covered 
with a red coating, which is stated to have been an oxide of iron. 

We have carried out some experiments with this material, a sample of 
which from the Bultfontein mine was kindly given to us by Mr. White 
of the Geology Department, University of Bristol. The specimen was a 
small octahedron, with faces about 1 mm?; careful microscopic examina- 
tion revealed only one or two small, shallow trigons on each face. This 
stone was surrounded by pulverized blue ground in an alumina crucible, 
which was heated in graphite container by induction. A transparent 
quartz jacket enclosed the graphite, over which a slow stream of nitrogen 
was passed to minimize burning, and also any oxidation of the kimberlite. 

The graphite was heated to 1450°c in about five minutes. The temper- 
ature was measured by sighting an optical pyrometer on to the inside 
surface of the container; the temperature of the specimen itself may 
have been considerably higher than this figure, which should be regarded 
as a lower limit. The point at which melting of the blue ground began 
was also difficult to determine, because it started at the bottom of the 
crucible, but it appeared to be about 1200°c. Melting was accompanied 
by considerable effervescence, probably as a result of the evolution of 
water vapour from the hydrated components of the rock. In the first 
experiment this effervescence caused the formation of a frothy cap of 
material over the top of the open crucible. Heating was continued for 
about ten minutes, and the sample was then cooled fairly rapidly by 
switching off the heater. 


+ The temperature scale does not seem to have been. specified by Luzi, but 
Gardner F. Williams (The Diamond Mines of 8. Africa) states that it was 
Reamur. It seems remarkable, however, that any diamond should survive for 
half an hour at 2200°c. 
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The porous material formed by the solidified froth was chipped away, 
and the diamond (which seemed to have been encased in a bubble of the 
melt) was recovered. It was blackened by a surface layer of graphite, 
and one side—presumably that actually in contact with the fused mineral 
throughout the heating, was considerably corroded with no recognizably 
crystallographic etching. This face had also a reddish tinge. The other 
faces exhibited a considerable number of sharp-edge triangular pits in 
‘negative’ orientation, some of them several microns deep (fig. 3). 

Since the diamond had been graphitized to an unknown depth, the 
possibility that the markings had been etched in the graphite could not 
be excluded. A second crystal, similar in size and quality to the first 
was heated in blue ground. The temperature rise was even more rapid 
than in the first experiment: 1450°c was attained in about two minutes, 
and the fused mineral did not form a cap obscuring the interior of the 
crucible. The heating was continued for seven minutes, at the end of 
which time it was noticed that the diamond had floated up to the surface 
melt. The heater was then switched off and the crucible cooled as rapidly 
as possible by increasing the flow of nitrogen. 

This stone was found to be deeply etched, but with no graphite coating. 
(It is possible that surface graphitization had occurred and that the 
graphite was burnt away by a trace of oxygen in the nitrogen atmosphere 
after reaching the top of the melt.) The microstructure of the diamond 
near an apex of the octahedral face is shown in fig. 4, Pl. 86. The surface 
(previously almost featureless) is covered with shallow pits in the ‘negative’ 
orientation, characteristic of the natural trigons on diamond, fairly 
uniform in size and outlined by what seem to be grooves, giving a curious 
“wormeaten’ appearance to the crystal. There are also a number of pits 
which are extremely deep; the particular cluster in fig. 5, Pl. 87 are about 
15 » in depth, and others of 20, and more have been seen. Figure 5 
is a photograph taken with the microscope focused on the bottom of 
another group. The shape of the pits seems to be roughly that of a steep- 
sided inverted pyramid, with a curious bump (indicated by a bright spot 
in the centre of the triangle) at the apex. It will be noticed that in spite 
of the very rapid rate of etching at these pits the corners are quite sharp, 
indicating that the ‘negative’ etch habit is very strongly marked. 

The face which had been resting on the surface of the molten kimberlite, 
and had thus been in contact with it during the entire period of heating 
was, as before, so heavily attacked that any crystallographic etching 
had been obscured. This face, like the corresponding face in the previous 
experiment, was a reddish colour. 


§ 4. Ercn anp Dissotution Hasrr or KIMBERLITE 
Since pits can be etched in negative orientation, and Frank et al. have 
shown that features characteristic of dissolution are associated with 
trigons, it seems reasonable to assume that trigons are etch figures. It 
follows that the various curved shapes of diamond are also the result of 
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dissolution of the octahedral form, since any reaction which etches the 
negative (110) steps forming trigons will also develop the same type of 
step at the edge of the (111) face, as shown in fig. 6. The edges are thus 
rounded off and the octahedral face is gradually replaced by a number of 
curved surfaces. Fersmann and Goldschmidt (1911) showed how this 
dissolution habit leads to the so-called dodecahedral and hexakisocta- 
hedral forms of crystal. 

Since all the observed surface features and morphology of diamond 
can be accounted for by dissolution processes, it is of interest to know at 
what stage in the history of the crystal these could have occurred. 


Cc B 


Formation of trigon (abc) and dissolution of octahedral face from edges (A BC}. 


§ 5. Buur Grounp anp THE Natura History or Diamonp 


Blue ground, or kimberlite, is found in South Africa in pipes and dykes 
formed by eruptions which probably occurred towards the end of the 
Cretaceous period; the upper layers have since been extensively eroded, 
possibly by some thousands of feet. The existing surface layers have 
been weathered to ‘yellow ground’. Kimberlite is a complex mineral 
whose composition varies with the time and place of eruption, but in all 
cases it clearly derives from material having a high proportion of olivine, 
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(Mg, Fe),SiO,. This has been largely altered by hydration to serpentine, 
H,Mg,Si,Qy. 

The rock has a porphyritic structure : that is, it consists of a fine grained 
matrix, in this case of mica and olivine, in which are embedded much 
larger crystals. In kimberlite these may be of, for example, olivine, 
phlogopite (Magnesia—mica), enstatite and diopside (pyroxenes), garnet, 
and occasionally diamond. Since any of these frequently occurs as 
inclusions in most of the others, it appears that they were all formed under 
the same conditions of great pressure and temperature at depths of several 
hundred kilometres, and were transported by the rise and eruption of 
fluid magma to their present position. After the eruption the liquid 
phase rapidly solidified to form the matrix of small crystals. 

Two facts are relevant to the present inquiry. First, it was remarked 
by Williams (1932) that diamond inclusions in diamond are often rounded 
octahedra. Obviously a fluctuation in carbon saturation in the neigh- 
bourhood of a growing crystal caused partial resorption, after which 
conditions for further crystallization were restored. 

A second observation worth noting is that most of the large transported 
crystals show clear signs of partial dissolution, arguing some universal 
change in conditions in the parent magma. This is generally assumed to 
have been a gradual fall in pressure during the rise of magma; it seems 
reasonable, therefore, to suppose that it was responsible also for the 
natural etching of diamond. It has been shown that at atmospheric 
pressure diamonds both dissolve and graphitize fairly rapidly in molten 
kimberlite, and that in an environment free of oxygen graphite forms at a 
temperature below 900°C in a few hours. Since stones are commonly 
found free of graphite (and from the traces which the formation of graphite 
leaves on a surface) and relatively lightly etched, the inference is that 
immediately after the eruption the kimberlite pipes were cooled by at 
least 400-500°c, and possibly much more, in a period of the order of 
minutes. This was no doubt the combined result of the release of pressure 
and of the consequent boiling off of water vapour and carbon dioxide. 
This cooling must in fact be an essential condition for the preservation of 
large diamonds; since the rate of cooling was presumably greatest at the 
surface, the upper layers may have been much richer in diamonds than 
the existing pipes. These upper layers were no doubt the source of 
alluvial diamonds. It is significant that the yield of all mines decreases 
with depth. 

The particular reaction responsible for etching trigons cannot yet be 
certainly identified. However, the available evidence indicates that 
dissolution took place in the fluid magma, and in such an environment 
anes must readily reduce the Fe,O3 present. It seems likely, therefore, 

his reduction is the principal agent of natural dissolution of diamond. 
The reddish substance deposited or adsorbed on the surface of diamonds 
artificially etched in kimberlite, also suggesting that this oxide is involved, 
lends some support to this hypothesis. 
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ABSTRACT 


Measurements have been made of the transmission of neutrons of wave- 
lengths between 5 and 10 A through 10 em of unirradiated and pile irradiated 
Brazilian quartz crystals. The total transmission cross section of the 
unirradiated crystal varied as o,=0-55-+1-63A in the measured wavelength 
range. After irradiation with 6-4x 1018 fast n° cm~*, as indicated by the 
[S(n°p)P] reaction, additional scattering was found. This amounted to 
0-18 barns per molecule between 12 and about 9A rising continuously 
below 9 A. The measurements are discussed in relation to optical absorption 
measurements and radiation damage theory. 


§ 1. IyTRODUCTION 


Antat et al. (1955) have shown how the concentration of defects in pile 
irradiated graphite can be determined from the induced long wavelength 
neutron scattering. The concentration is calculated in terms of experi- 
mentally determined parameters but a relatively large irradiation dose 
(ca. 101 fast n° em-*) is required to observe the effect. 

Because of the possibilities of collision induced recombination occurring 
during such heavy irradiations, for example, a displaced atom filling a 
neighbouring vacancy, direct comparison of the experimentally determined 
concentration with low dose damage theory is probably not valid. The 
object of the present experiment was to use the method to calibrate some 
other property—optical absorption in quartz—which could then be 
compared with radiation damage theory at lower doses. 

In addition an attempt has been made to measure the wavelength 
dependence of the induced scattering. This is important because it was 
assumed by Antal et al. (1955) that the scattering is by isolated point 
defects. For wavelengths of the order of angstréms we should expect 
the scattering by isolated nuclei of dimensions of the order 10-44 to be 
independent of wavelength. 


§ 2. EXPERIMENTAL 
The experiment became possible in this country following the develop- 
ment by Dr. P. Egelstaff and his colleagues (1957) of a small volume of 


+ Communicated by Professor R. W. Ditchburn. 
t Present address: Diamond Research Laboratory, South Africa. 
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liquid hydrogen moderator in the Harwell pile BEPO. These authors have 
already described the details of the neutron source, liquid air cooled beryl - 
lium filter, chopper and counting equipment. 

In our experiment the neutron transmission has been measured of two 
erystals (A and B) of Brazilian quartz having dimensions 2x 2x 10cm. 
The crystals were supplied by Hilger and Watts Ltd. and were described 
as being of the highest optical quality. The uniformity of the visible 
absorption induced by irradiation is known to be an indication of quality, 
although the band is associated with aluminium impurity. In crystal B 
the visible absorption after irradiation was completely uniform throughout 
the crystal. 

The neutron beam was defined by two cadmium sheet stops theapertures 
of which were just less than the specimen dimensions. A specimen holder 
was built which allowed each of the crystals to be moved reproducibly 
into the beam and which provided one straight through position. 
Measurements of the transmitted intensities J, and J, and of the incident 
intensity J) were made, counts being taken usually for 10 minutes at a 
time and repeated in a sequence (I, I) I, 1, I) Ip...). 


§ 3. RESULTS 


Both crystals were measured in the wnirradiated condition. At this 
stage although some counts were made for \ > 104 the greatest numbers of 
counts were obtained betwen 4 and 104. Approximately 10* counts 
were made on the beam transmitted by each specimen at each of six wave- 
lengths between 4-4 and 6-04, and 2x 103 at six wavelengths between 
6-2 and 104. The measured transmissions decreased with increasing 
wavelength from 73% to 55%, the measurements on the twe crystals 
agreeing within + 1°% and frequently within +4%. These results indicated 
that the crystals could be assumed to be identical for our purpose. 

Crystal B was then irradiated in a hollow uranium slug in BEPpo. This 
facility was kindly provided by the Metallurgy Division and the dose 
was estimated from the *2S (fast n°, p)?2P reaction. The threshold for 
this reaction is about 2mev and the dose received by the quartz was 
6-4 x 1018 n°em-2 of energy between 2 and 4mev. The temperature 
during the irradiation was approximately 40°c. . 

Subsequently the transmissions of both A and B were pea atte 
At this stage many more counts were made, varying from 10° for each 
crystal at wavelengths near 8A to 1-8 x 104 near 114. For the whole 
wavelength range we take the measurements on crystal A to be those for 
the unirradiated state and compare with those for B after irradiation. 
The transmission of the irradiated crystal was between 5% less, near 
6A, and 2°% less near 11 A, than the unirradiated crystal. ae 

No attempt was made to record accurately the change in length o 
crystal B with irradiation, but any such change must have been less than 
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0-5%. Inany case we have used the length w, cm measured after irradia- 
tion. Then we have: 

I,=I,exp(—o,%,N), 

Ig =1, exp (— op%pN’) =Lo exp ( —op tpl) 
where N and N’ are the numbers of SiO, molecules per cm* before and 
after irradiation, and the o’s are total cross sections of all processes leading 
to removal of neutrons from the beam. 

The densities of the two crystals were measured by weighing in air and 
water and the change was found to be: 

Ap|p= — (0-784 + 0-002) %. 
A similar value of — (0-81 + 0-10)% was found by flotation in tetrabromo- 
ethane/ether for a 0-2 mm disc which was irradiated at the same time and 
in the same can as the large crystal. 

The results for o, are shown in fig. 1 (a) for wavelengths between 5 and 
124. Some measurements were made below 4A using a crystal spectro- 
meter and the transmission decreased rapidly with decreasing wavelength, 
corresponding to the occurrence of increasing numbers of Bragg reflections. 

In fig. 1 (0) the quantity (c,’—«o,) is plotted as a function of wavelength. 
For all measurements the resolution of the chopper was such that the 
wavelength spread in the beam at each wavelength was about 0-8 A. 


§ 4. DiscussIon oF RESULTS 
4.1. Unirradiated Quartz 

If the specimen had been polycrystalline the limiting wavelength (2d) 
for Bragg reflections would have been 8-54. However, the orientation of 
the crystal was such that an a-axis (e.g. the normal to (120)) was 6° from a 
line parallel to the beam. The orientation was determined from a back 
reflection Laue photograph. In this orientation we expect Bragg scattering 
to occur at 7-76 and 6-88 A corresponding to (110) and (100) type reflections. 
These were not observed. Bragg reflection might be expected to be spread 
in a good quartz crystal over a wavelength range of about 0-01 A, but the 
wavelength spread of the beam was 0-8 A so that it is not surprising that 
other peaks were not resolved. Measurements at lower wavelengths 
using the spectrometer showed a strong peak near 3-5 A, corresponding to 
a (200) reflection. 

From fig. 1 (a) it will be seen that o, increases linearly with wavelength, 
but the extrapolation to zero does not go through the origin. Between 
5 and 124 we find 

o, (barns per molecule) = 0-55 + 0-163. 
This cross section will include contributions from absorption (Si), inelastic 
scattering by lattice vibrations, and scattering by impurity or other defects. 
Isotope and spin scattering are negligible for SiO,. 

According to Bacon (1955) the absorption cross section of Si is 0-06 b. 
per atom at 1-084 and since absorption is directly proportional to A we 
have a contribution to o4 of 0-056, b. per SiO, molecule. The theory of 


}otai Cross section 
(barns per molecule) 


Gp ay: 
(barns per molecule) 


Scattering of Long Wavelength Neutrons by Quartz 1283 


inelastic scattering has not been worked out exactly for quartz. However 
Cassels (1950) calculated for a body-centred cubic crystal whose Debye 
temperature (#) was 453°K and whose atoms were of mass (M) 56, that 
H(S)/S increased linearly with wavelength above about 54. (8S) is the 
inelastic scattering cross section and S the coherent: scattering cross section. 
He found that H(S)/S at a given wavelength was approximately propor- 
tional to 1/176. Using these results reasonable agreement is obtained with 


Fig. 1 


Neutron wavelength (A) 


(a) The variation of the transmission cross section of quartz with neutron 
wavelength: ....... before irradiation ; x x x x x after irradiation. 
(b) The variation of the radiation induced scattering cross section 
with wavelength. 


the measurements on aluminium (Cassells 1951), magnesium (Squires 1952) 
and carbon (Kothari and Singwi 1957). For quartz we have Ss; =2-0b. 
per atom, S,=4-2 b. per atom (Bacon 1955) and @=760°K. Then regarding 
the crystal as being composed of atoms of mass 20 units and S=3:5b. 
we find, using Cassel’s result, that in the region beyond about 
54 H(S)=0-11A, b. per molecule. 

These two contributions, therefore, would be expected to give a cross 
section for \>ca. 5 4 of o oc 0:166A, b. per molecule, which is in fortuitously 
good agreement with the experimental value of 0-163A. 
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We have now to consider whether the wavelength independent part of 
the total cross section can be ascribed to impurities or other defects. 
Absorption is wavelength dependent and need not be considered. 
Hydrogen seems to be the only likely impurity having a sufficiently high 
scattering cross section, 81 b. per atom, to give the required contribution. 
Approximately 0-2 atm %, of H would have to be present. This is a higher 
concentration than one would have thought probable, but there is little 
quantitative data available and the number is feasible. 


4.2. Irradiated Quartz 


The additional scattering induced by the irradiation may be represented 
by an effective cross section per molecule (aq) given by: 


1 14 1 is) 
pi, 29 eee il IES, SO ate ye Lg 
Bae Bites WieAmne oe) va (72 ae 
or if AN/N <1 by: 


Pn oc Le ee EE n 
Oqg—Op, -0.= 5 {2-Im(7*) xm (7)} . . ( ) 


It will be seen from fig. 1 (6) that, for wavelengths greater than 9 A, (ap’— a4) 
is constant within the experimental spread and is equal to 0-16 + 0-03 b. per 
molecule. As shown by the density measurements AN/N =0-8°% and the 
true o,—given by eqn. (1a)—should be greater than (c,’—o,) by 0-02 
for wavelengths between 9 and 124. Thus we have co, in this region 
equal to 0-18 + 0-03 b. per molecule. 

For wavelengths less than 9A the induced scattering cross section 
increases with decreasing wavelength. This indicates that the scattering 
in the damaged crystal cannot be ascribed only to isolated point defects. 
Four factors may contribute to the wavelength dependence : 

(i) Scattering by atoms in the amorphous regions which may be 
produced near the ends of tracks of knocked on atoms (Mitchell and Paige 
1956). 

(ii) Scattering by the atoms surrounding point defects, which atoms 
will not be expected to be in their normal lattice positions. 

(ili) Scattering by aggregates of point defects. The constituents of such 
aggregates would have to be separated by about 54 or less. 

(iv) A change in the inelastic scattering cross section arising from a 
radiation induced change in the vibrational spectrum of the crystal. 

The last factor could be eliminated by measuring the transmissions at 
temperatures low enough for inelastic scattering to be negligible. It is 
hoped that this may be done in a future experiment. For the present we 
note that a change of 30° in the Debye temperature would not introduce 
a serious error and that changes of this magnitude have not been reported. 


Following Antal e¢ al. (1955) it is assumed, therefore, that any change in 
the inelastic scattering may be neglected. 
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The net effect of (i) and (ii) is that a number of atoms, while remaining 
part of the SiO, network}, are displaced from their normal positions. 
It is likely that the increased scattering from these factors will be in the 
normally diffracting region, less than 9 A. 

Finally we should expect that aggregates, such as close pairs of like 
interstitial atoms having separations of about 5A, would scatter less at 
wavelengths beyond 9A than would two isolated interstitial atoms. 
If the effect is important the defect concentration calculated below would 
be an underestimate of the number of displaced atoms. 

With this qualification we interpret o, for wavelengths greater than 9 4 
as arising from scattering by point defects alone. Assuming that the dis- 
placed atoms are in the ratio of 2 oxygen to 1 silicon then: 


oq (b. per molecule) = 29(c5; + 200), tomes ae (2) 


where g is the fraction of atoms displaced, and o,; and oy are the scattering 
cross sections for silicon and oxygen nuclei. The essential feature of the 
method of Antal ef al. (1955) is that the atomic scattering cross sections are 
known from other nuclear scattering experiments. Bacon (1955) gives 
oy, = 2-0 and o9=4:2b. peratom. Thus for o,=0-18 + 0-03 b. per molecule 
we find g= (0-86 + 0-10)%. 

Using the result found by Kinchin and Pease (1955) that a primary 
knock on of energy H# produces #/2H, displacements we find that 
6-4 x 1048n° cm? having energies between 2 and 4Mev would produce a 
total fraction of displaced atoms of ca.10%. At such a concentration 
replacement collisions have to be considered so that a direct comparison 
with the experimental value is not possible. 

We make the comparison with low dose theory as follows. There is an 
absorption band in the ultra-violet (C-band, 5-7ev) induced by pile 
irradiation. For low doses the absorption increases approximately 
linearly as pmax(em)!=Dx 10-"’, where D is the equivalent thermal 
dose for the irradiation conditions described by Mitchell and Paige (1956). 
Crawford (private communication) have shown that for doses corres- 
ponding to jzmax>100cm~! the rate curve departs increasingly rapidly 
from the linear rate and eventually goes through a flat maximum. This 
could also be associated with the occurrence of replacement collisions. 

Theabsorption in the thin specimen, which was irradiated in the same can 
as the larger crystal, was 194 cm~ at the C-band peak and had an estimated 
width at 4-absorption of 1-3ev. The width agrees as well as can be 
expected with the value of 1-4ev found at lower doses. Assuming that 
the fraction (zx), of defects in the state of ionization which gives the 5-7 ev 
absorption, is the same for the two conditions of irradiation and that the 
f-value is independent of the concentration of centres, we have: 


fey =9nn Gun eet ae te ete (3) 
where g,,, is the concentration of displaced atoms determined above— 
the subscript HD referring to ‘high dose’—and gp is the fraction which must 


+ Point defects are not parts of the network. 
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be present to produce each 1 em of absorption in the low dose experiments. 
This required an equivalent thermal dose of 1017 n° cm. Thus g,, 18 
4-4 x 10-8% per D=10!7 n° cm~. 

According to the calculations of Mitchell and Paige (1956), 9,5 per 
1017 equivalent thermal dose is between 30x 1017/N and 45x 10!"/N 
where N, the total number of Si and O atoms per cm3, is 7°86 x 10. By 
calculation, therefore, we find that 3-8 x 10-°% <g,p < 5:7 x 10°08 

Finally, the measurements of the heavily irradiated crystals lead to a 
value of (af) equal to 3-7 x 10-? for the C-band absorption, where f is the 
oscillator strength for the transition. 


§ 5. CONCLUSIONS 


The main conclusions of this work are: 

(i) The transmission of long wavelength neutrons through crystalline 
quartz can be reasonably accounted for by inelastic scattering and Si 
absorption. 

(ii) The scattering between 5 and 12A induced by pile irradiation 
cannot be explained in terms of isolated point defects alone. 

(iii) The scattering for A> 9 A taken in conjunction with optical absorp- 
tion data leads to a reasonable numerical agreement with radiation damage 
theory. Because of the various uncertainties, however, the agreement is 
not as good as the numbers indicate—there possibly being an uncertainty 
of a factor of about 5. 

ACKNOWLEDGMENTS 


We should like to thank Dr. P. Eglestaff, A.E.R.E., Harwell, for making 
it possible for us to do this experiment. We are also grateful to him and to 
various members of his group for their advice—particularly to Mr. G. 
McCallum and Dr. F. J. Webb. 

For his interest and help in obtaining the crystals we wish to thank 
Professor R. W. Ditchburn, and P. T. W. is grateful to Dr. J. F. H. Custers, 
Director of the Diamond Research Laboratory, Johannesburg, for per- 
mission to work at A.E.R.E. during the Summer of 1957. 

Finally, we should like to thank Mr. J. D. Rigden for assistance with the 
density and optical measurements, and Dr. R. J. Elliott for helpful 
discussions. 

REFERENCES 
AntaL, J. J.. Diengs, G. J., and Weiss, R. J., 1955, Phys. Rev., 99, 1081. 
Bacon, G. E., 1955, Neutron Diffraction (Oxford : University Press). 


Butrerworts, I., Kcerstarr, P. A., Lonpon, H., and Wess., F. J., 1957, 
Phil. Mag., 2, 917. 


Senay ee M., nee Progress in Nuclear Physics: I, ed. O. R. Frisch 
ondon : Butterworth; Berlin: Springer): 1951 

Neoey ia pringer) ; , Proc. roy. Soc. 
Kinoutn, G. H., and Pass, R. 8., 1955, Rep. Prog. Phys 18, 1 

Kormarr, L. §., and Stvew, K. §., 1957, Phys. Rev., 106, 230. 

MitcHeLL, E. W. J., and Patcs, E. G. S., 1956, Phil. Mag., 1, 1085. 
Squirus, G. L., 1952, Proc. roy. Soc. A, 212, 194, ; 


[1287 3 


Activation Energies for High Temperature Creep of 
Polycrystalline Zinc} 


By W. J. M. Trcgarr 
Department of Metallurgy, The University, Sheffield 
and OLEG D. SHERBY 


Department of Metallurgical Engineering, Stanford University, California; 
formerly Department of Metallurgy, The University, Sheffield 


[Received July 1, 1958] 


ABSTRAOT 


The creep of polycrystalline zine above 0:5 Tin (Ty,= absolute melting 
temperature) appears to be controlled by two thermally activated processes 
acting additively. One process, predominating between 0-5 and 0-6 eS 
believed to be associated with dislocation climb since the experimentally 
determined activation energy for creep is about equal to the energy for self- 
diffusion, namely, about 21-0 kcal per mole. The other process, pre- 
dominating between 0-8 and 1:0 7,,, is believed to be associated with 
prismatic slip since the activation energies for these two processes are both 
equal to about 38-0 keal per mole. 


§ 1. IyTRODUCTION 


EXTENSIVE experiments on the high temperature creep of polycrystalline 
pure metals reveal that the experimental activation energy for creep, Q,, 
is about equal to the respective activation energy for volume self- 
diffusion, Q, (Sherby etal. 1954, Frenkel etal. 1955, Weertmanand Shahinian 
1956, Sherby and Lytton 1956, and Sherby 1958). This equality appears 
to be well established for twelve pure metals, namely, the face-centred 
cubic metals, Al, Ni, y-Fe, Au, Cu, and Pb; the body-centred cubic 
metals, «-Fe and B-Tl; the hexagonal close packed metals, Mg, «-TI, and 
Cd; and face-centred tetragonal indium. A recent publication on self- 
diffusion in platinum suggests that @,~@, for this metal also (Y,=56 kcal 
per mole (Sherby e¢ al. 1954) and Q,= 67 keal per mole (Kidson and Ross 
1957)). These correlations suggest that creep of most polycrystalline 
metals is controlled by dislocation climb. The activation energies for creep 
of single crystals under single slip conditions do not conform with the 
activation energies for volume self-diffusion (Gilman 1956, Lytton et al. 
1958). Such evidence suggests that dislocation climb might not be the 
important process in these cases. Nag 
Although the creep of zinc has been extensively studied, the activation 
energy for creep does not appear well defined. All known experimental 
data on activation energies for creep and self-diffusion of zine are given in 
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tables | and 2. Sherby et al. (1954) suggested that Q, for polycrystalline 
zine was about equal to Q, (about 23 kcal per mole) from analyses of 
limited data on creep of zinc obtained in the temperature range 40-70°C 
by Pomp and Lange (1936). However, Cottrell and Aytekin (1950) 
obtained a value of Q, for polycrystalline zine considerably lower than 
@,, in fact about 12 kcal per mole in the temperature range 95-120°c. 
Further, Gilman (1956) suggested that creep of polycrystalline zinc, at 
least above 250°c, was controlled by prismatic flow for which Q, was 
equal to 38 keal per mole... In view of the discrepancy in the results for 
low temperatures and the possibility that @, for polycrystalline zinc may 
be temperature dependent above 0-5 7,,, in contrast to the other metals 
so far examined, the present investigation was undertaken to determine 
Q, over the temperature range from 70°—400°c., 


Table 1. Activation Energies for High Temperature Creep of Zinc 


Temp. Creep 
Aealincole} Reference Purity (%) en ene 
A. Polycrystalline zine 
Pomp and . 
26 ; 7 as 
26-0 Lange (1936) 99-0 40-70 | 450-750 
Cottrell and 
». Y i= ‘ Q 5 
12-0 Aytekin (1950) Spec. pure 95-120 | 80-120 
38-0 ‘ ss 
Gaited) Gilman (1956) 99-999 298 15-40 
B. Single crystals of zine 
F Cottrell and : 
29-0 (Basal) eae (1950) Spec. pure 85-110 2-5 
5-0 (Basal) Gilman (1956) 99-999 180-407 | 0-8-2 
38-0 (Pris- : 
Bhan Gilman (1956) 99-999 | 255-407 | 40-100 


§ 2. MATERIAL AND EXPERIMENTAL TECHNIQUES 


The material used was ‘Crown Special’ zine (> 99-99 wt% zinc) with 
traces of silver (0-0001%), cadmium (0-002%), copper (0:0002%), iron 
(0-003°%) and lead (0-0005%), and was received in the form of sananste 
bar, §in. diameter. These bars were warm rolled at approximately 
150° to j@ iu. diameter and then machined into solid cylindrical Specimens 
+ in. in diameter and 3in. tall. The initial batch was annealed in an 
argon atmosphere at 395°c for 1} hours. The resulting average grain 
size was determined to be about 0-06 in., although some fine grains were 
still present in the final annealed material. Later, a further batch was 
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a 


annealed at 403°o for 1 hour in a salt bath (KNO, and NaNO) yielding 
a smaller average grain diameter of 0-015 in. 

Creep tests were performed under compression since suitable equipment 
was available. This equipment is described in another paper on the 
creep properties of alpha and beta thallium (Sherby 1958). It was 
possible to measure creep strains to a sensitivity of 0-00025, and all strains 
recorded in this investigationare reported as truestrain ¢, where « =In (1+), 
e being the engineering strain. The stress was maintained constant by 
periodically adding lead shot to the direct load weight pan as creep 
proceeded. It was possible to maintain the stress constant to +4°% of 
the predesired stress throughout any one test. Tests were usually 
terminated by «=0-15 and virtually no barrelling of the specimen was 
observed to take place up to such strains. 


Table 2. Activation Energies and D, for Self-Diffusion in Zine 


Single crystal Polyerystal 
Investigators and 
reference ee eee 
Il QL Dll Dol @ D 
kcal/mole cm?/sec kcal/mole cm?/sec 
Jaumot and Smitht 
(1956) 22-0 | 24-9 | 0-076 | 0-39 22-7 O19 
Jaumot and Smitht 
(1956) 22-0 | 23:5 | 0-068 | 0-12 23°5 0-11 
Shirn et al. (1953) 21-8 | 24-3 | 0:13 0-58 — — 
Miller and Banks (1942)} 20-4 | 31-0 | 0-05 0-09 —- —- 
Liu and Drickamer 
(1954) 19:6 | 25:9 | 0-021 | 1-60 _ — 
Gertsriken and 
Pryanishnikoy (1957) | — 26-4 4:3 
+ Sectioning technique. t Absorption technique. 


Temperatures were maintained constant by using either (a) silicone oil 
baths for the lower temperature creep tests (60°C to 250°C) or (6) salt baths, 
containing equal mixtures of KNO, and NaNOs, for tests above 0 
Specimen temperatures were maintained constant to within + 5°. 
Cyclic temperature creep test techniques (Tietz and Dorn 1956, Sherby 
1958) were used to determine creep rates under a given constant stress 
and from these data activation energies were calculated. A typical 
example of such a creep test is shown in fig. 1 where the creep rates as a 
function of strain for a given stress, under cyclic temperature conditions, 
are recorded. The observed experimental activation energies, (,, were 
calculated from the relationship 

eens h Wea oi) Ve Sen ae ON OF) 
é, exp(—@,/RT,) 


Ss 
P.M, 4 
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where ¢,=creep rate at given creep stress and strain at temperature Te 
in degrees absolute, and €,=creep rate at same stress and strain at 
temperature 7’, and &=gas constant = 1-987 calories per mole. 
Metallographic observations were carried out on flat surfaces prepared 
on specimens by mechanical abrasion followed by chemical polishing in 
Vinaver’s solution (Vinaver and Druelle 1955) till all the surface 


deformation was removed. 
Fig. 1 


22:3 kcal/mole. 


3 = 2000 PSII. 
10 o 


Ol 02 ‘03 ‘04 -OS 06 -O7 “08 Oo 10 ‘ll 


Creep rate (e)—true creep strain (¢) curves for use in evaluation of the activation 
energies for creep of pure polycrystalline zine (99-99°%) under constant 
Stress. 


§ 3. RESULTS AND Discussion 


All of the values for the experimental activation energies for creep of 
polycrystalline zinc determined in this investigation are recorded in fig. 2. 
The temperature range covered was from about 350°K up to about 680°K. 
No clear cut variation of the activation energy for creep was noted with 
the creep strain, ef. the typical example in fig. 1, and furthermore there 
appeared to be no systematic variation of @, with the creep stress (fig. 2). 
However, fig. 2 reveals a very definite dependence of Q. on temperature. 
The results suggest the presence of at least two processes with differing 
activation energies. A process with a low activation energy appears to 
be operating below 0-6 7',, while a process with a high activation energy 
operates above 0-8 7',,. These results are thus quite different from those 
obtained for other pure polycrystalline metals where only one thermally 
activated process appears to operate above about one-half of the absolute 
melting temperature (Sherby et al. 1954). 

Considering first the low temperature results, it is seen that Q. for Zn 
is about 21 keal per mole which is approximately the value reported for 
self-diffusion (table 2). This equality between Q, and Q, suggests that 
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creep of polycrystalline Zn in this temperature range may be controlled 
by some type of dislocation climb process. Previous workers (Cottrell 
and Aytekin 1950, Cahn et al. 1953-54, Ramsey 1951-52) have examined 
polycrystalline zinc specimens after slow straining in this temperature 
range and have observed basal slip and subgrain formation, the subgrain 
boundaries being at right angles to the slip plane. Similar structures were 
observed in the present investigation. These observations suggest that 
the rate controlling process in this temperature range is the climb of edge 
dislocations away from pile-ups in the basal plane to form subgrain 
boundaries. The mechanism of such a climb process in h.c.p. metals has 
not yet been established. One of the difficulties associated with the 
development of a mechanism is the nature of the barriers in the basal 
plane. If only basal slip is operative, the blocking of edge dislocations by 
sessile dislocations is impossible. However, in polycrystalline aggregates 
of h.c.p. metals, the deformation is inhomogeneous and it seems likely 
that the barriers here are sessiles formed by interaction of basal and 


non-basal slip. 
Fig. 2 


as computed from two thermally — 
activated processes acting additively. 


40 


os) 
O 


coarse fine 
grained g, (PS!) grained 


Q. (kcal/mole) 
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Effect of temperature on the observed experimental activation energy for creep 
of polyerystalline zinc. 


Cottrell and Aytekin (1950) have also determined activation energies 
for creep of pure zinc polycrystals and single crystals in the low temperature 
range, but their results (table 1) do not agree with the values obtained in 
this investigation. This discrepancy is not clearly understood, but a 


482 
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possible explanation is that Cottrell and Aytekin only determined their 
activation energies over a small range of temperature, 25°C, and a small 
error in determination of creep rates would have resulted in a large error 
in the value of Q,. 

Considering now the high temperature process, it is seen that above 
0:87», Q, for zine is about 38-0kcal per mole. This value is identical to 
the value obtained by Gilman (1956) for creep of pure single crystals 
oriented for flow on the prismatic plane. Attempts were made to observe 
crystallographic slip after creep at 375°C. However, in harmony with the 
observations of Cahn et al. (1953-54), no visible slip could be associated 
with creep at such temperatures. This observation suggests that very 
fine slip may be operating at these temperatures. It is also possible that 
surface diffusion may obliterate the slip lines which would normally be 
visible. The fact that identical values of Q, were obtained for creep of 
polycrystalline zinc and prismatic creep of zinc single crystals, however, 
constitutes strong indirect proof that prismatic slip controls the creep of 
zine polycrystals above 0-8 7’,,,. 

The data in fig. 2 suggest the presence of two thermally activated 
processes. These two processes can contribute to creep by acting either 
additively or sequentially. Each process can be described by an equation 
of the form: 

€é,=w,;exp (—Q,/RT), o= Const. . .. aaee eee oe 


where €,=creep rate for ith process, W,=the weight factor for the ith 
process, and Q@;=thermal activation energy of ith process. If creep is 
controlled by two thermal processes acting sequentially then the observed 
rate of creep is primarily determined by the slowest process. If the 
thermal activation energies and W’s are different, an inversion in the rate 
controlling process can occur as a function of temperature. An exact 
analysis of the creep rate under such circumstances is difficult since it 
requires a knowledge of the dependence of creep rate of each process on 
the creep stress. However, this approach predicts a decrease in Q, as 
the temperature is increased and therefore cannot explain the present 
results on the creep of zinc. 

On the other hand, if two thermal processes act additively, the creep 
rate, at a given stress, can be described by the equation 


e=W, exp (—@,/RT)+w,exp(—Q,/RT). . . . . (3) 


This equation predicts a low value of activation energy at low temperatures 
and a high value at high temperatures with intermediate values of @ in 
the intermediate range of temperatures. The results on zine (fig. 2) agree 
with such an additive process. If we assume that eqn. (3) can describe 
the creep of zinc above 0-5 7',,, a value for the observed activation energy 

Q.,,, can be calculated, as . 


a, -— Bdlne _ wiQs exp (~ Qi RT) +0,Q, exp (— Qu RT) 
*  dQ/P) ~~ w, exp (—Q,/RT) +0, exp(—Q,/kT) 
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If one arbitrarily assigns a value of unity for W,, W, can be calculated 
from a knowledge of Q,,,, in the intermediate temperature range. Sub- 
stituting Q,,,=30-0keal per mole at 7'=500°K into eqn. (4), W, was 
found equal to 2-5x 107. It is now possible to calculate the variation of 
the observed activation energy with temperature giving the dotted curve 
shown in fig. 2. It is interesting to note that eqn. (4), with the above 
values, predicts a rather broad temperature range (150°c) of intermediate 
values of @,,,, in agreement with the experimental observations and this 
suggests that the proposed hypothesis might be substantially correct. 

The major difficulty of such a model is the large ratio of W,/W, needed 
to explain the experimental results. Recent work by Weertman 
(1955, 1957 b) makes it possible to calculate a theoretical value for W,/W,. 
Under conditions where dislocation climb is the controlling process, he 
has developed an equation for steady state creep which can be expressed 
as follows (Weertman 1955): 

SO Orpl = Gal jo «5 ee lb) 
where @,=activation energy for self-diffusion, c=creep stress and S is 
a structure constant to 

1207 oD, 
pwekT 
where ¢ is the critical stress to operate a Frank—Read source, L” is related 
to subgrain size, L’=grain size, M=number of Frank—Read sources per 
unit volume, =shear modulus, k= Boltzmann’s constant, 7'= absolute 
temperature and D)= diffusivity constant. Assuming values appropriate 
to zinc}, we obtain 
€=2°8 x 10° exp (—@,/RT7')o*?, hr. ear eee €(6) 
Furthermore, Weertman (19574) has also attempted to explain the high 
activation energy for prismatic creep of zinc single crystals by using 
Gilman’s suggestion (Gilman 1956) that the controlling process might be 
that of overcoming a large Peierls stress on the prismatic plane. Weert- 
man’s equation for creep under such conditions is similar to that for creep 
controlled by climb and can be approximately described for the case of 
zinc by ; 
é€=5-7 x 108 exp (—Q,,/RT')o*®, hr ao Pyaar aap 
where @ , = 38000 calories per mole and o is the creep stress in kg/cem?, 
Thus with Weertman’s models for creep, the ratio W,/W, is equal to 
about 200, which differs considerably from 107 as determined in this 
work. If the interpretation of the results obtained by means of two 
thermally activated processes is correct, then it is apparent that the 


"—(.02em, M=109/cm3, L'=0-lem, o=1 kg/cm? from Gilman’s 
oe ae creep Gas 1956), #=3-8 x 10° kg/em? (Gilman 1956), 
k=1-4x 10-22 kg em/°c, 7’=400°K, D,=360 cm?/hr (table 2). 

+t Weertman’s theoretical equation is considerably more complex a ee 
(7), but his calculated curve of € exp (Q)/RT) vo can be accurately describec 
by eqn. (7) for creep stresses ranging from 4 to 200 kg/cem?. 
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theoretical analyses will have to be modified before a complete description 
of high temperature creep in zinc can be given. Furthermore, Weertman’s 
equations predict that the creep rate should vary with (creep stress)*, but 
this is clearly not the case for the present work as shown in fig. 3. A 
recent modification of the original dislocation climb equation by 
Weertman (1957 a) has yielded a stress exponent equal to 4-5. Calculations 
of W,/W, from eqn. (7) and Weertman’s modified climb equation are now 
dependent upon the creep stress. At o=70kg/em?, W,/W,=7 x 10° 
which is still considerably below 10’. 


coarse 
grained 


: a ae. ee ee eee 
10 eee 7 a Misi. 
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Kuffect of stress on the steady state creep rate of polycrystalline zine (99-99%) 
at high temperatures. 


The results of this investigation raise the question of whether other 
hexagonal close-packed metals have similar changes in their activation 
energies for creep. One might expect that the creep behaviour of those 
metals whose slip systems change with temperature would reveal several 
activation energies. Unfortunately little data are available to check this 
idea. In magnesium, a change from basal to pyramidal slip occurs at 
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about 220°c (0-537) (Schmid 1931) but analyses of creep tests in the 
range 0:4 to 0-657',, (Sherby and Frenkel 1954) suggest that the 
activation energy for creep is constant and equal to that for self-diffusion. 
Tt seems desirable to determine activation energies for creep of magnesium 
from 0-5 to 1:07), using the cyclic temperature method, to check and 
extend these results. The creep of polycrystalline cadmium in the range 
0-5 to 0:87), (Frenkel et al. 1955) is associated with a constant activation 
energy about equal to that for self-diffusion. However, no observations 
of slip mechanisms appear to have been reported for cadmium at high 
temperatures. For polycrystalline thallium, creep in the range 0-5 to 
0-87), (Sherby 1958), at which temperature it transforms to a body- - 
centred cubic structure, is associated with an activation energy equal to 
that for self-diffusion. However, no investigation of slip systems in 
thallium has yet been reported. Clearly, further work is necessary to 
establish the relation of creep to crystallographic flow mechanisms in 
hexagonal metals. 
§ 4. CONCLUSIONS 


1. The experimental activation energies found for creep of pure poly- 
crystalline zine above one-half the absolute melting temperature 
tentatively suggest that creep of this material is controlled by two 
thermally activated processes acting additively. 

2. At low temperatures, 0-5 to 0-67), a low activation energy process 
is dominant, with @,=21-0keal per mole. This value suggests that creep 
is controlled by dislocation climb. 

3. A high activation energy process is dominant in the temperature 
range 0-8 to 1:07, with Q,=38-0kcal per mole. It is suggested that 
creep in this temperature range is controlled by crystallographic flow on 
the prism plane. 

4. Recent creep theories do not describe completely the high temperature 
creep of polycrystalline zinc. 


ACKNOWLEDGMENTS 
This work was carried out during the tenure of a University Research 
Fellowship (W.J.T.) and a U.S. National Science Foundation Post- 
Doctoral Fellowship (O.D.S.). The authors wish to thank Professor 
A. G. Quarrell for provision of facilities, and encouragement during the 
work. The authors are also grateful to Imperial Smelting Corporation, 
Avonmouth, for analyses of the zine used in this investigation. 


REFERENCES 


Caun, R. W., Buar, I. J., and Bett, R. L., 1953-54, J. Inst. Metals, 82, 481. 

Corrret, A. H., and Avtexin, V., 1950, J. Inst. Metals, 77, 389. 

FRENKEL, R. E., Suersy, O. D., and Dorn, J. E., 1955, Acta Met., 3, 470. 

GERTSRIKEN, S. D., and PryantsHnikov, M. P., 1957, All Union Conference on 
the Application of Radioactive and Stable Isotopes in the National 
Economy and Science, Moscow Edition of the Academy of Sciences, 


SSS R. 


1296 On the Creep of Polycrystalline Zinc at High Temperature 


Giuman, J. J., 1956, Z'rans. Amer. Inst. min. (metall.) Engrs, 206, 1326. 

Jaumor, Jr., F. E., and Smrrn, R. L., 1956, Trans. Amer. Inst. min. (metall.) 
Engrs, 206, 137. 

Krpson, G. V., and Ross, R., 1957, International Conference on Radioisotopes, 
Paris. 

Lav, T., and Drickamer, H. G., 1954, J. chem. Phys, 22, 312. 

Lyrron, J. L., Sauparp, L. A., and Dorn, J. E., 1958, Trans. Amer. Inst. min. 
(metall.) Engrs, 212, 220. 

Miuuer, P. H., and Banks, F. L., 1942, Phys. Rev., 61, 548. 

Pomp, A., and Lanes, W., 1936, Mitt.-K.-Wilh.-Inst. Hisenforsch. 18, 51. 

Ramsey, J. A., 1951-1952, J. Inst. Metals, 80, 167. 

Scumip, E., 1931, Z. Hlektrochem., 37, 447. 

SHmrsy, O. D., 1958, Trans. Amer. Inst. min. (metall.) Engrs (to be published). 

SHERBY, O. D., and Frenxen, R. E., 1954, Trans. Amer. Inst. min. (metall.) 
Engrs, 200, 690. 

SHeRBY, O. D., and Lyrron, J. L., 1956, Trans. Amer. Inst. min. (metall.) 
Enyrs, 206, 928. 

SHERBY, O. D., Orr, R. L., and Dorn, J. E., 1954, Trans. Amer. Inst. min. 
(metall.) Engrs, 200, 71. 

Surrey, G. A., Waspa, F.S., and Huntineton, H. B., 1953, Acta Met., 1, 513. 

Trerz, T. E., and Dorn, J. E., 1956, Trans. Amer. Inst. min. (metall.) Engrs, 
206, 156. 

VINAVER, W., and DruzELLE, P., 1955, Rev. Met., 52, 612. 

WeERTMAN, J., 1955, J. appl. Phys., 26, 1213; 1957 a, Ibid., 28, 362; 1957 b) 
Ibid., 28, 1185. 

WEERTMAN, J., and SHaninran, P., 1956, Trans. Amer. Inst. min. (metall.) 
Engrs, 206, 1223. ¥ 


On the Delay Time in Plastic Flow of Indium Antimonidet 


By J. W. ALLEN 


Services Electronics Research Laboratory, Baldock, Herts. 
[Received April 15, 1958; and in revised form June 6, 1958] 


ABSTRACT 


The delay time before the onset of plastic deformation has been studied as a 
function of temperature and stress for indium antimonide in uniaxial com- 
pression. The temperature range was 150-350°c, and the stresses were 
about 10°-107 dynes/em*. An empirical relation of the form 7= const. 
oT exp (H/kT) has been found. A theory of dislocation production is ~ 
advanced to explain these results, 


§ 1. LyrRopuctTIon 


WHEN a small force is applied to a single crystal of indium antimonide 
there is a delay before flow begins (Allen 1957). A similar phenomenon is 
observed in germanium and silicon (Gallagher 1952). Various interpreta- 
tions of this have been suggested, and it was the aim of the present work to 
provide data on the variation of the delay time with stress and temperature 
in order to distinguish between these interpretations. 'The observations 
have, in fact, suggested a new model which seems to be in better accord 
with the experiments than previous ones. 

Since plastic flow is to be interpreted in terms of a dislocation mechanism 
both the number of dislocations and their movement will be involved. 
We have chosen to concentrate on the former in our experiments since the 
onset of plastic flow in these materials is marked by a very rapid increase in 
dislocation density. The etch-pit technique offers a convenient way of 
measuring this density. 

§ 2. EXPERIMENT 


In the experiment a single crystal block was compressed under a given 
stress at a given temperature and the time under these conditions which 
elapsed before the dislocation density increased was measured. This 
was repeated for a range of temperatures and stresses. 

The specimens were rectangular blocks, 2-8 mm x 0-9 mm x 0-6 mm, 
with a (111) face on one of the long narrow sides. They contained an 
average of fifteen dislocations each (determined by etch-pit count) and 
were cut from a region of the bulk crystal where the dislocation density 
was reasonably uniform. They were compressed in a piston and cylinder 
arrangement in air. Under the conditions of the experiment only very 
light oxidation occurred even at the highest temperatures used. In 
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order that the stress should be a pure linear compression the piston was 
made to a sliding fit in its cylinder, and soft lead platens wrapped in thin 
pure aluminium foil were used to distribute the stress uniformly and to 
prevent frictional forces along the surface. (The aluminium wrapping 
was necessary to prevent alloying between the lead and the indium 
antimonide: for one reading at 350°C pure aluminium was substituted 
for the lead.) 

The cylinder and piston and the specimen were heated separately in a 
furnace to the required temperature, after which the specimen was placed 
between its platens in the cylinder and a load was applied. After a suitable 
time the load was removed and the specimen was taken from the furnace, 
thereby quenching it to room temperature in a few seconds. It was 
polished chemically and etched to show dislocations (for details of these 
procedures and of 111 etch see Allen (1957)) and examined under the 
microscope. If flow had not begun it was replaced in the furnace and 
loaded for a further time, and the sequence repeated until flow was observed 
to have begun. The delay time (also known as the incubation period) 
was taken to be the total time under load before the onset of flow. 

Some points here require further attention. If the specimen is heated 
and then etched with 111 etch without polishing then etch-pits are 
produced which are irrelevant to our present purpose, i.e. they do not 
mark dislocations which penetrate into the bulk crystal. It might be 
objected that dislocation loops are formed at the surface and would in 
course of time move through the crystal. Our results would then be 
misleading, since such loops are removed at each polishing. However, it 
is found that within the limits of the experiment the delay time is inde- 
pendent of the number of periods of heating (each followed by etching) 
of which it is the sum: we therefore conclude that in this experiment at 
least surface sources of the type described should be negligible. These 
remarks only apply to specimens with carefully prepared surfaces and 
with a reasonably high number of dislocations present in the bulk. With 
a smaller initial number of dislocations the results are governed by the 
behaviour of surface sources. In practice, after a few preliminary experi- 
ments the delay time expected under a given set of conditions could be 
estimated, and the heating periods used were one of about two-thirds 
the estimated time followed by further periods of about a sixth of that 
estimated, these being repeated until flow began. In this way about three 
or four etchings were used for each sample whatever the actual length of the 
delay time, and this meant that the decrease in size of the specimens due 
to etching was approximately the same for all of them. 

In the initial stages of deformation dislocations were not produced 
uniformly throughout the crystal, but appeared in localized groups which 
for reasons which will become apparent later we call avalanches. This is 
shown in. fig. 3, Pl. 88. On tracing the dislocations through the bulk 
material it is found that considerable cross-slip has taken place. Linear 
arrays in well-defined glide planes did not occur until a much later stage. 
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Itis thought that this effect is characteristic of the deformation mechanism 
and is not due to a lack of homogeneity in the stress: this is, of course 

difficult to prove, but all possible precautions were taken to ae the See 
homogeneous, and the avalanches did not seem to be particularly associated 
with the upper and lower surfaces of the specimen, where stress concentra- 
tions may be expected to arise. 


Fig. 1 


ISOO 


1000 


NUMBER OF DISLOCATIONS 
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i5 30 AS 60 
Time (MINUTES) ——- 
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Etch pit count as a function of time of compression at 300°C for a series of InSb 
specimens. The solid line is that of eqn. (2). 


In order to find how the dislocation density increases with time a 
subsidiary experiment was performed in which a series of InSb blocks were 
held at the same temperature under the same stress for different times. The 
etch-pit count as a function of time of heating is given in fig. 1 and it will 
be seen that the onset of flow is reasonably well defined. In the main 
experiment flow was deemed to have occurred when the dislocation density 
had increased by a factor of about four as judged by eye (corresponding to 
the appearance of an average of two avalanches). This somewhat 
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arbitrary criterion leads to consistent results, since a factor of two in the 
dislocation density when flow has begun corresponds to only a small change 
in the time, as may be seen from fig. 1. 


$3. RESULTS AND ERRORS 


The values of the delay time as a function of temperature and load are 
presented in the table. The load given therein is the actual weight used 
in the compression : the orientation of the crystal was such that 100 gm wt 
gives a shear stress of 1-6 x 10° dynes/cm? on the (111) plane most favourably 
oriented for slip. The stress is proportional to the load within the accuracy 
set by the scatter in specimen area. The variation in this area was + 6% 
initially, and the area decreased by 1% on each etching. However, since 
the number of etches used was approximately the same for each sample, 
and since they were applied at roughly the same fraction of the delay time, 
the total inaccuracy in the stress, taken as proportional to load, should 
not be greater than +8%. The temperatures were measured with a 
calibrated chromel—alumel thermocouple inserted in the cylinder and 
were controlled to +5°c. Most of the temperature variation was in the 
first few minutes after loading the specimen. The method of estimating 
the delay time described above should give values with an uncertainty 
not greater than + 15% and the scatter of values obtained when a number 
of samples are compressed under the same conditions is within these 
limits. 


Delay Time in Minutes as a Function of Applied Load and Temperature 


Load (gm wt) 


Temperature 
(°c) 47 66 160 260 455 900 
150 == 210 60 
200 — — 200 — 30 = 
250 210 70 25 10 D — 
300 40 — 6 — 0-75 — 
350 — — 1 — — — 


§ 4. THEORY 
The localized arrays of dislocations formed in the initial stages of 
deformation suggest an avalanche mechanism, i.e. that dislocation sources 
are produced by the interaction of dislocations, and that these sources 
produce further dislocations which interact to produce more sources, and 
soon. Ifsis the number of sources per cm?, n is the density of dislocations 


per cm? and v is their velocity, then in the model the rate of production of 
sources follows a bimolecular law 
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where 8 is some constant. As a rough approximation we assume that each 
source gives off a constant number y of dislocation loops before these 
produce a new source and then work-hardens. Observation of the 
dislocation arrays shows that this number is of the order of 10 in these 
experiments. Under this assumption the rate of increase of dislocations 
will be proportional to the rate of increase of sources, that is dn/dt= avn. 
One then has a 1 1 

- tee komt st. een) 


ee cl CUNgt 1—t/7, 
where m9 is the initial dislocation density and 7,=(avn,)-!. Equation 
(2) is plotted as the full curve in fig. 1 and is seen to predict a rapid increase 
of dislocation density after an initial delay period which is a fraction 
nearly equal to unity of 7,. 

The mathematical description contained in eqns. (1) and (2) is open to 
some criticism. Equation (1) is statistical in nature and therefore cannot 
obtain in the early stages of the compression, for the dislocations initially 
are so few and so far apart that some time must elapse before any inter- 
actions occur at all. It is also implicit in eqn. (1) that all the dislocations 
are able to move and that their velocities are randomly distributed in 
space. If those dislocations lying near one preferred glide plane may 
move more rapidly than others then the number of interactions is more 
nearly proportional to the first power of n, giving an exponential increase 
with time in place of the expression in eqn. (2). Again, at a later stage 
of the deformation work-hardening will set in and this will prevent the 
divergence to infinity of n as t approaches 7,. However, the treatment 
should provide a useful approximation to the actual behaviour of the 
system. 

In the experimental part (§ 2) the incubation period was taken to be 
the time at which n~4n,, so from eqn. (2) one finds 7 to be #7,. Since 
this is inversely proportional to the dislocation velocity we require for the 
interpretation of our results this velocity as a function of stress and 
temperature. If the dislocations are surrounded by an impurity atmo- 
sphere, as seems reasonably certain in these materials, the velocity will be 
limited by the drag of this atmosphere. Cottrell and Jaswon (1949) 
have given a theoretical treatment of the effect and have formally solved 
the problem. They express the impurity atmosphere distribution as a 
series of hyperbolic Mathieu functions, and show that the force on the 
dislocation can be expressed as the gradient of the potential energy of an 
impurity atom, integrated over all the impurity atmosphere. Given the 
parameters of interest, namely the temperature 7’, the diffusion coefficient 
of impurities D, the dislocation velocity v and the coefficient A which 
occurs in the expression for the interaction energy V’ of an impurity atom 
in the strain field of a dislocation, namely 


Ve AVSir Oey ee en es, ete ey (8) 


one can find by numerical methods the force F required to maintain the 
dislocation in a state of steady motion, However, their treatment is not 
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suited to the inverse problem of determining the parameters from the 
experimental results, nor does it lead to explicit relations between the 
parameters. 

We have therefore taken the basic equations of Cottrell and Jaswon’s 
formulation, namely the differential equations describing the movement of 
impurities and the conservation of their number together with the integral 
which relates the force on the dislocation to the perturbation of the 
impurity atmosphere, and by suitable transformation of variables we 
have derived the following functional relation : 


vA fF 
Sa -1(=5) : . . ° . . . (4) 


100 


IDT (ARBITRARY UNITS) — 


x 

| 4 
SO Tete) 500 1000 
APPISIE Dea LOAD (ams. wt.) = 


Docantinnn eho ioe See : : : 
Logarithmic plot of -D7 against applied load, showing the functional relation- 
ship between these quantities, 
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in which , is the number of solvent atoms in unit volume, and Co is the 
impurity concentration far from a dislocation. These two factors and 
the quantity A vary only slightly with temperature and may therefore 
be taken as constants in our experiment. Remembering that the incuba- 
tion period + is inversely proportional to the dislocation velocity, one 
finds 

TtDT =g(F). RE eke ae) 


§ 5. COMPARISON OF THEORY AND EXPERIMENT 


It is expected that the diffusion constant of impurities will be of the 
form Dp, exp (—H/kT’) where D, and E£ are independent of temperature. 
From eqn. (5) it follows that H can be found from the slope of a plot of 
In r7’ against 7! at constant stress F. This procedure yields a value of 
0-78eyv. Ifone then plots 7DT against F, using this value for EZ, the results 
should lie on a single curve. All the results given in the table are plotted in 
this manner in fig. 2. It is seen that a good representation of them is 
given by the equation 

TD aCOnst Nena eee er ee (0) 


in accordance with the functional requirement of eqn. (5). 


§ 6. Discussion 


Not much significance is to be attached to the exact value of the 
exponent 1-74 in eqn. (6) since this is a power-law approximation to an 
expression which is expected to be much more complex. In fact we have 
performed a similar analysis of some results of Treuting (1955) for 
germanium and have found an expression of the type of eqn. (6) with an 
exponent 1-37. One point, however, deserves attention and this is that 
previous treatments of this type of microcreep have been based on the 
assumption that the velocity is proportional to the stress, in accordance 
with Einstein’s diffusion relation (Cottrell 1953). But the number of 
impurities bound near the dislocation might be expected to decrease as 
the dislocation velocity increases: pictorially, the moving lattice washes 
away some of the impurity atmosphere. This gives a velocity increasing 
with stress faster than linearly and this agrees with our result. 

So far the drag on the dislocation has been attributed to impurities 
but it may well be that the production of interstitials and vacancies is of 
greater effect. There are then two processes to consider, namely the 
production of defects and their diffusion away. If the former is the 
limiting process than Friedel (1956) has shown that the velocity should 
be linearly related to the force by the Einstein equation, while if the 
latter process is limiting then the treatment we have given should still 
be correct with defects taking the place of impurities. It may be noted 
that the diffusion energy for self-diffusion of antimony in InSb is 0-75 ev 
(Boltaks and Kulikoy 1957), not significantly different from the value 
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0:78 ev which appearsinourresults. Thisis almost certainly a coincidence, 
since the same correlation between the energies does not hold for germanium. 

Seitz (1952) has proposed a rather different explanation of the delay 
time in which he assumes the existence of Frank—Read sources which 
initially cannot multiply because the active length is locked by impurities. 
If a stress is applied then thermal vibrations enable the source to free itself 
from the impurities and slip begins. In this model the activation energy 
derived from a plot of In r against 7! is interpreted as the energy binding 
impurities to a dislocation. Although the experimental values, of the 
order of an electron-volt, are rather high for this interpretation this may 
be due to the addition of a bonding energy to the elastic interaction. 
Seitz treats only a single source and does not consider how the effects 
of a number of sources will add to give the behaviour of the whole crystal. 
Unless all the sources acting are roughly of the same length the sharp 
onset of flow due to a single source will be smeared out. This objection 
can be met by assuming that the probability of occurrence of a Frank—Read 
source is a rapidly decreasing function of its length, for example, 
exp (—1/b). Then since at a given stress there will be a minimum length 
for an active source the active source lengths will have a distribution 
sharply peaked at a certain value, and so all sources will begin to multiply 
at roughly the same time. A more serious objection to Seitz’ theory 
is that it would seem to predict the appearance of well defined slip lines in 
the initial stages of deformation followed later by interaction configura- 
tions, contrary to the observations. 

The crucial experiment which would distinguish clearly between the 
mechanism of Seitz and the avalanche mechanism would be that of 
determining dislocation velocities as a function of stress and temperature, 
rather than the delay time. This has not been found possible so far 
because the presence of impurities at a dislocation causes its etch-pit to 
be rather diffuse, and this renders the unambiguous identification of 
individual dislocations in our crystals as grown extremely difficult. The 
dislocations introduced by deformation do not have dense impurity 
atmospheres, of course, and so can be identified easily. A related problem 
prevents the determination of the relation between delay time and 
dislocation density predicted by eqn. (2), for until the locking impurities 
are identified it will not be possible to produce crystals with different 
dislocation densities but the same impurity atmosphere concentration. 
Correlation between 7 and mn, can therefore only be qualitative at this 
stage, but crystals with high dislocation density do seem to have short 
delay times. 

No hypothesis has been advanced above concerning the detailed 
mechanism of dislocation multiplication in the avalanche, i.e. as to how 
interacting dislocations can produce sources so oriented that the newly 
produced dislocations can also interact. A number of such mechanisms 
are possible in principle, but more detailed information on dislocation 
dispositions is needed to determine which are acting in practice, 
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ABSTRACT 


Oriented deposits of hexagonal ice have been formed by condensing water 
vapour on to silver iodide substrates at temperatures of about — 100°c and 
examined by electron diffraction methods. A feature of the diffraction 
patterns, obtained from the oriented hexagonal ice, is the appearance of 
strong spots in the ‘forbidden’ hk.4 positions. At lower substrate 
temperatures only unoriented deposits of cubic ice were obtained. 


§ 1. INTRODUCTION 


Some ten years ago Vonnegut (1947) discovered that, in suitable 
circumstances, small particles of silver iodide were highly effective as 
nuclei for the formation of ice crystals. When particles of silver iodide 
were introduced into a supercooled cloud in laboratory experiments, ice 
crystals began to form at a temperature of about —5°c and at about — 15°c 
most of the particles of silver iodide appeared to serve as ice nuclei. (If no 
nuclei are present, the spontaneous formation of the ice crystals does not 
occur until the temperature falls below —40°c.) Vonnegut’s discovery 
has led to several large-scale efforts to produce artificial rainfall by 
‘seeding’ suitable supercooled clouds with smokes of silver iodide. 
Although it is not possible to say whether these attempts have been 
successful, the ice nucleating properties of the silver iodide are beyond 
doubt and it is of interest, therefore, to investigate the manner in which 
ice is formed on a single crystal of this substance. 

Vonnegut chose silver iodide for his nucleation experiments with ice 
because of the geometrical similarity of the hexagonal lattice of silver 
iodide and the hexagonal lattice of ice. If the hydrogen atoms of the ice 
lattice are ignored, the atomic arrangement consists of two interlocking 
‘close-packed’ hexagonal arrays of oxygen atoms, the dimensions of the 
unit cell being at 0°c, 

G@=4:52 4; CS TBI A. 
The hexagonal lattice of silver iodide also consists of two similarly 
interlocking ‘close-packed’ arrays, one formed by the silver atoms and 
the other by the iodine atoms. At room temperature the unit cell 
dimensions for hexagonal silver iodide are 


a= 4-59 A; c=T7-51 A, 


{+ Communicated by the Author, 
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The dimensions quoted above have been obtained from the experimental 
results of Megaw (1934) for ice, and Wilman (1940) for silver iodide, by 
multiplying their values by 1-002 to convert from x.v. to &. 

The fact that the geometry of the crystal lattices of ice and silver iodide 
are so closely related suggests that the first step towards establishing the 
mechanism of nucleation is to conduct experiments to see whether there 
is any tendency for ice to grow asa single crystal, oriented in some particular 
manner, on a single crystal of silver iodide. Jaffray and Montmory (1957) 
have observed oriented ice deposits on silver iodide by means of an optical 
microscope and electron diffraction methods have been used, for the first 
time it is believed, in obtaining the results given below. These electron 
diffraction experiments were conducted at temperatures below —90°c; 
it is impossible to maintain thin deposits of ice above this temperature 
as they evaporate very rapidly under the vacuum conditions necessary 
to operate an electron diffraction camera. 


§ 2. PREPARATION OF SILVER [ODIDE SUBSTRATES AND EXPERIMENTAL 
PROCEDURE 


Because it is difficult, if not impossible, to obtain suitable crystals of 
silver iodide prepared by the usual methods of crystal growth, recourse 
to the following procedure was made. Silver was evaporated from a hot 
tungsten wire in vacuo and allowed to deposit on to thin pieces of freshly 
cleaved mica which were maintained at a temperature of about 270°C 
during the deposition process. This produced monocrystalline and 
apparently atomically flat layers of silver having their (111) faces parallel 
to the cleavage face of the mica (see Newman and Pashley 1955). The 
pieces of mica carrying the layers of silver were then subjected to chemical 
attack by suspending them in an atmosphere in equilibrium with a weak 
solution of iodine in water. By this means the oriented silver layers 
were converted to oriented layers of silver iodide. It was necessary to 
conduct the exposure of the silver to the iodine vapour in darkness, as 
it was found that daylight, particularly when bright, brought about the 
decomposition of the silver iodide, leaving a mainly randomly arranged 
deposit of silver crystallites as the final product. 

The chief orientations of silver iodide (which frequently occur as a 
mixture of two forms—one cubic and the other hexagonal—at room 
temperature), observed by means of electron diffraction on the silver 
layers, were as follows: 

(a) (111) cubic AgI parallel to Ag (111), 
[110] or 110 parallel to Ag [211]. 
(b) (00.1) hexagonal AgI parallel to Ag (111), 
[10.0] parallel to Ag [211]. 
These orientations and also others which are less common have been 
examined by Pashley (1951). Figure 1, Pl. 89 shows an electron diffraction 
pattern obtained from an oriented layer of silver iodide. This pattern 
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is common to both cubic—[2II] azimuth—and hexagonal—[11.0] 
azimuth—forms of silver iodide. Not infrequently the cubic form 
predominated although occasionally patterns consistent with a largely 
hexagonal orientation were obtained. A pattern from a cubic deposit— 
[110] and [110] azimuths—is shown in fig. 2, Pi 29, | 
Specimens of silver iodide so prepared were mounted in a special low 
temperature unit which enabled the silver iodide to be maintained at any 
desired temperature down to about —170°c. The low temperature unit 
was then placed in the electron diffraction camera and deposits of ice 
were formed from water vapour injected on to the silver iodide substrate 
from a glass apparatus connected to the camera. Ice, contained within 
the glass apparatus and maintained at a suitably low temperature, 
constituted the source of water vapour. The complete apparatus and 
method of obtaining such deposits of ice are described more fully by 
Blackman and Lisgarten (1957). 


§ 3. EXPERIMENTAL RESULTS AND DISCUSSION 


It is convenient to consider the results of the ice deposition experiments 
in three parts. The first deals with deposits on the silver iodide substrate 
at an initial temperature of about —100°c or a little higher. As soon as 
the deposits were formed the temperature of the substrate was allowed 
to fall as rapidly as possible to about —130°c so that the diffraction 
patterns could be observed and recorded without fear of the ice evaporating. 
In this series of experiments orientation of the ice was obtained. A pattern 
recorded from such a deposit is illustrated by fig. 3, Pl. 90 which shows 
diffraction spots due to hexagonal ice—[11.0] azimuth—superimposed 
on the corresponding array of diffraction spots due to the silver iodide. 
(Figure 4 gives the indices of the diffraction spots due to the ice.) In this 
particular experiment, the substrate consisted of a largely cubic orientation 
of silver iodide and the identity of the hexagonal ice was easily established 
by examining the diffraction pattern associated with the other principal 
azimuth. Thus, several of the diffraction spots due to the oriented ice 
—ice [10.0] azimuth—cannot be associated with corresponding diffraction 
spots due to the cubic silver iodide. When the temperature of the substrate 
was raised, it was possible to observe the disappearance of the diffraction 
spots due to the hexagonal ice as the ice evaporated. The remaining 
pattern due to the cubic silver iodide was generally similar to that shown in 
fig. 2. Another example of a pattern from an oriented deposit of hexagonal 
ice is shown in fig. 5, Pl. 90. Here the deposit was somewhat thicker 
than that of the deposit illustrated by the pattern of fig. 3. The diffraction 
spots as observed on the fluorescent screen of the camera were relatively 
bright and the ring system which can also be seen in fig. 5 is consistent 
with a random array of hexagonal ice crystallites. On warming, the 
bright diffraction spots and rings due to the ice were observed to disappear 
leaving a poor pattern of weak spots from the silver iodide together with 
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some rather diffuse rings indicative of the formation of silver ; this resulted 
from the decomposition of the silver iodide under the action of the electron 
beam. 

The second type of experiment was conducted with the substrate at an 
initial temperature of about —130°o, in which region it was expected that 
cubic ice would be formed. Only unoriented deposits of crystallites of 
cubic ice were observed on the silver iodide in this case. For the third 
kind of experiment, the substrate temperature was lowered to about 
—165°c and in this region ‘amorphous’ ice was obtained. Deposits of 
this type were then heated slowly until the cubic form appeared. No 
orientation of the ice was observed in these experiments. 


Fig. 4 
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Analysis of the pattern shown in fig. 3, Pl. 90. The indices given refer to the 


spots due to the hexagonal ice. 


The fact that there appears to be a temperature limit (at about — 100 °c) 
below which orientation of the ice on the silver iodide cannot be obtained 
is not entirely unexpected and analogous cases are known. For example, 
silver forms an unoriented deposit ou a mica base at room temperature, 
but forms excellent monocrystalline layers on mica which has been heated 
to about 270°c. What is perhaps surprising is that orientation of the ice 
on the silver iodide does occur at such low temperatures ; it seems reasonable 
to conclude, therefore, that at higher temperatures the silver iodide is a 
very efficient ‘orienting agent’. 
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A prominent feature of the diffraction patterns obtained from the 
oriented hexagonal ice is the occurrence of fairly strong diffraction spots 
in the ‘forbidden’ hk.4 positions. These can be observed in figs. 3 and 5 
and their positions are indicated by the diagram in fig. 4. In the case 
of fig. 3, the hk.4 ice spots are just separated from the silver iodide spots 
and since the cell dimensions of the ice are slightly smaller than the corres- 
ponding cell dimensions of the silver iodide, the spots due to the ice 
appear at slightly greater distances from the central spot of the diffraction 
pattern than do those due to the silver iodide. As in the case of the 
‘forbidden’ rings (e.g. hexagonal 00.4 and cubic 222) of the electron 
diffraction patterns obtained from unoriented aggregates of ice crystallites 
(Honjo et al. 1956, Blackman and Lisgarten 1957) there arises the problem 
of accounting for the corresponding ‘forbidden,’ diffraction spots in these 
hk.4 positions in the patterns from the oriented layers of hexagonal ice. 
The 00.4 spots are the hexagonal counterparts of the 222 reflections which 
have been observed in the case of single crystals of silicon (see Heidenreich 
1950) as well as in the case of diamond itself. The 222 spots observed in 
the reflection patterns from silicon and diamond are almost certainly due 
to dynamical scattering effects. In the case of ice one must expect on the 
basis of the Pauling (1935) ice model that the hydrogen atoms play some 
part in producing these hk.4 spots. It is unlikely that the hydrogen 
atoms are solely responsible for the effects observed in the reflection 
patterns from the oriented ice deposits. The hk.4 spots are almost of 
the same intensity as the others, whereas if one examines a transmission 
pattern from unoriented hexagonal ice of roughly the same thickness, 
the hk.4 rings are very weak compared with other rings. If the relatively 
strong hk.4 spots were produced solely by diffraction from hydrogen 
atoms, there would be no reason to expect the corresponding rings from 
the transmission patterns to be relatively so weak. It is therefore 
concluded that the hk.4 spots in the patterns from the oriented ice are 
produced mainly by dynamical interaction or multiple diffraction from 
oriented but different crystals. 

Measurements made on reflection patterns of the type shown in fig. 3 
indicate that, at a temperature of about —130°c, the dimensions of the 
lattice of silver iodide are between 2 and 3° greater than dimensions of the 
ice lattice. According to Blackman and Lisgarten (1957), the lattice 
dimensions for ice at —110°o are 


a=4:49348: c=7-337 4. 


These values are about 24° less than the room temperature lattice 
dimensions of silver iodide (Wilman 1940) already quoted. It follows, 
then, that the lattice dimensions of silver iodide do not change appreciably 
between room temperature and —110°c. To check this point, ring patterns 
from silver iodide were recorded at the above two temperatures and the 
diameters of the rings compared. There was no difference between the 


two within the errors of observation, which were estimated to be about 
0+2.0/ 
ie 
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ABSTRACT 


In the etching of recrystallized aluminium (both of 99-99% purity, and 
zone-refined material) in Lacombe and Beaujard’s etchant, the grain 
boundaries are found to resist attack immediately after electropolishing, 
but to become susceptible to severe attack after periods of a few days at 
room temperature. It is suggested that this behaviour is due to diffusion 
of impurities along the grain boundaries, into the oxide film. 


§ 1. INTRODUCTION 


Tr is generally agreed that etch pits occur on the surface of aluminium at 
places where there is a change, either physical or chemical, in the composi- 
tion of the oxide film. Such changes may be associated with the points 
of emergence of dislocations; the presence of impurity precipitates; or 
more probably, combinations of the two. Clear evidence that dissolution 
can occur where dislocations meet the surface is given by the alignment of 
pits along the traces of sub-boundaries on aluminium crystals in the work 
of Lacombe and Beaujard (1948), Lacombe (1947) and Cahn (1949). The 
pits composing these alignments are sometimes irregularly spaced, which 
suggests that the presence of a dislocation alone is insufficient to promote 
etching. Instead, it is more reasonable to suppose that an etch pit will 
form only on a dislocation upon which impurities have precipitated. 
This conclusion, that regions of the surface susceptible to attack are 
produced by the combination of impurities and dislocations was reached 
by Wyon and Crussard (1951) and also by Forty and Frank (1955). 
When aluminium polycrystals are etched with Lacombe and Beaujard’s 
etchant soon after polishing, pits rarely occur along the grain boundaries. 
Specimens grown from the melt show characteristic ‘clear lanes’ along 
the grain boundaries; these lanes being bounded by high density regions 
of etch pits on rapidly grown polycrystals (Forty and Frank). This does 
not occur when the samples are prepared by recrystallization ; instead, the 
pits are fairly uniformly distributed over the surface, different grains 
being discernible by the variations in the orientation of the etch figures 
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and the appearance of some grain boundaries as lines on the surface. 
In the above treatments grain boundaries are usually not severely attacked. 
Lacombe and Beaujard (1948) observed that intercrystalline attack only 
developed on boundaries between crystals having large differences in 
orientation. Lacombe (1947) has produced marked dissolution along 
grain boundaries after prolonged immersion in a 10°, solution of hydro- 
chloric acid. Again Wyon et al. (1956), using Lacombe and Beaujard’s 
etchant have demonstrated the dependence of boundary etching on 
purity and quenching treatments. 

In the course of the work described here, it was found that by giving 
polycrystalline aluminium suitable room temperature anneals between 
polishing and etching with Lacombe and Beaujard’s etchant the grain 
boundaries were made susceptible to attack so that they were marked by 
lines of pits usually overlapping one another. This result suggests that 
during the interval between polishing and etching, impurities were able to 
diffuse along the grain boundaries to the surface, in sufficient quantities 
to change its constitution enough to allow it fo be readily etched. 


§ 2. EXPERIMENTAL PROCEDURE 


The observations were made partly on aluminium of 99-99°% purity 
but mainly on the same material further purified by the method of zone 
refining. The aluminium was refined in the same manner as described by 
Braun et al. (1958), although in this case the treatment was carried to a 
further stage, eighteen passes being given to the ingot. After the third 
and tenth pass a portion of the impure end was cropped and replaced by an 
equal amount of the starting material, so that the length of the rod remained 
the same. 

This treatment produced aluminium of sufticient purity to enable it to 
recrystallize at —45°c in six to seven days after a 95° compressional 
deformation. Three single crystal specimens from different parts of an 
ingot prepared in the above way, showed residual resistances of 1/3000, 
1/3500 and 1/3900 of the room temperature value. (These ratios. were 
measured by Dr. D. Shoenberg at Cambridge.) 

To facilitate the production of etch figures, portions of the metal were 
cold rolled and cut into plates of dimensions about 24x 1x0-2cm. These 
plates were then divided into three groups; groups A and B consisting 
of refined material and group C of 99-99% aluminium. Batch A was 
allowed to recrystallize at room temperature for a day, whereas B and © 
were annealed at 600°c for six days and then cooled in calm air. 

After these recrystallization treatments all the specimens were carefully 
polished, first on a rotating flat wheel covered with selvyt cloth soaked in 
Brasso and then electrolytically in a bath of 66% acetic anhydride and 
33% perchloric acid, cooled in ice. The surfaces were finally washed 
under a jet of water, rinsed in alcohol and dried with ether. 

Etching was carried out in Lacombe and Beauj ard’s etchant maintained 
below 10°c, and followed by the same washing procedure as after polishing. 
The manner in which the etching was performed was somewhat complicated, 
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each specimen being exposed to the etchant several times. Between each 
etch the surfaces were not repolished, but the aluminium was stored in a 
desiccator at room temperature. Again, some specimens received their 
first etch immediately after polishing, while others were stored for several 
days before treatment. This procedure made it possible to examine the 
etch patterns produced in the normal way, 1.e. on etching immediately 
after polishing ; and also to follow the changes produced by ageing at room 
temperature. 
§ 3. RESULTS 


The results obtained fall into three categories, namely : 


(a) Observations on surfaces etched immediately after polishing. 

(b) Observations on the same surfaces as (a) when etched again after 
various times of storage at room temperature. 

(c) Observations on specimens which received their first etch after 
ageing at room temperature. 


(a) Etching Carried Out Immediately after Polishing 


For the 99-99% pure specimens of group C, an etch of 5-10 seconds 
duration was enough to produce a high density of small pits, fairly uniformly 
distributed over the surface. The presence of grain boundaries had little 
marked effect on the distribution of pits; although very occasionally some 
pits lay along the boundaries. Usually, however, the boundaries were 
seen by the variations in the orientations of the pits or when the surface 
was viewed slightly out of focus. Figure 1+ shows a region in which the 
grain boundary pitting is a maximum for such a specimen. The surface 
is slightly out of focus to show up the grain boundaries. A feature of the 
photograph is the tendency for many of the pits to form lines; this is 
probably due to the fact that this specimen was not electropolished 
sufficiently long to remove all the traces from mechanical polishing. 

A marked change in the configurations was observed on the pure material 
comprising groups A and B compared to those seen in the case of group C. 
The general density of pits was much smaller, a fact which allowed the 
etching times to be considerably lengthened. In addition to the expected 
increase in resistance to dissolution, it was found that the distribution of 
pits was not uniform, but varied from one grain to another depending on 
orientation as found by Lacombe and Beaujard (1948). The least 
densities of pits occurred in grains whose surfaces were close to the (100) 
orientation, while the highest densities were found on those surfaces close 
to the (111) orientation, the orientations being estimated from the shapes 
of the pits, which tend to be bounded by cube faces with the etchant used. 
A striking example of this is shown in fig. 2, which shows a region of the 
surface of a specimen of group B (which had been annealed for six days 
at 600°C) after an etch of 60 seconds. Figure 3 shows a part of the severely 
attacked grain at high magnification, and reveals the orientation as fairly 
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close to (111). Two pits found on a neighbouring grain which was only 
slightly pitted, are shown in fig. 4; these pits being square and flat bottomed 
indicate that the orientation is near (100). 

In general, most of the grain boundaries on the surfaces of the pure 
material in groups A and B were visible as lines. Lacombe and Beaujard 
suggested that these traces were not steps due to unequal dissolution of 
neighbouring grains during polishing or etching, as they were unable to 
produce a relief effect using oblique illumination. In our case such steps 
were present on some boundaries after electropolishing, as was shown by 
interference fringes (fig. 5). Other boundaries, however, did not affect the 
fringes in any discernible way. Etching did not seem to affect the situation 
greatly, and no grooves along the boundaries could be detected. 

There were usually no pits formed along the grain boundaries of group A 
specimens (refined material recrystallized at room temperature), but 
occasionally small portions were attacked, this tendency being even less 
for group B aluminium which had been annealed at 600°c. 


(0) Further Etching after Various Room Temperature Ageing 


It has been mentioned that all the specimens when etched immediately 
after polishing showed a small amount of grain boundary attack. This 
only happened at a few places on the surfaces examined, the majority of 
boundaries being unaffected. The extent of this attack could not be 
increased by lengthening the etching time, as this merely produced deeper 
troughs at places where pitting did occur. However, when the surfaces 
were given another etch after suitable room temperature anneals, the 
majority of grain boundaries were delineated by troughs of overlapping 
etch pits. The periods of ageing required were shorter for the aluminium 
of group A, which had been recrystallized at room temperature, than for 
groups B and C which had been annealed at 600°C. 

In the case of the group A specimens, a time of six days between the 
etches was sufficient to produce widespread grain boundary attack. 
Figure 6 shows a region of maximum grain boundary ‘troughing’ produced 
by the first etch ; the same field of view after a second etch six days later is 
shown in fig. 7. Another typical result obtained after the second etch is 
shown by fig. 8. 

Group B specimens (refined and annealed at 600°c) showed little grain 
boundary attack when etched for the second time, six days after the 
first etch. ‘Troughing’ was however quite extensive following a third 
etch two weeks later (fig. 9). 

The behaviour of the impure aluminium comprising group C was more 
or less the same as that of group B with regard to the grain boundary 
attack. The first and second etches, separated by a room temperature 
anneal of six days, yielded little boundary pitting, while most grain 
boundaries showed continuous troughs after the third etch a fortnight 
later. Figure 10 shows a surface which was etched twice, once immediately 
after polishing which yielded a configuration similar to fig. 1, and again 


after three weeks ageing. 
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(c) Ageing between Polishing and the Furst Ktch 


A number of specimens from group A, which had been recrystallized at 
room temperature, were stored in a desiccator for six days after polishing. 
When they were etched at the end of this interval, it was found that the 
majority of the boundaries were pitted (fig. 11) but the extent of the 
attack was not as great as when the ageing period had been preceded by an 
etch immediately after polishing. It was not possible to increase the 
extent of the troughing by etching again without delay as this only deepened 
the troughs which already existed. However when a further ageing period 
of three days was allowed to elapse before etching dissolution was wide- 
spread along the grain boundaries, e.g. fig. 12 which shows the same field 
of view as fig. 11. : 


§ 4. INTERPRETATION 


If an etch pit is formed on the surface of aluminium where there is a 
singularity in the constitution of the oxide film, such as an impurity 
precipitate associated with the point of emergence of a dislocation, then we 
can explain the resistance to solution (when etched immediately after 
polishing) where a grain boundary intersects the surface by suggesting that 
this region does not possess sufficient impurities to render it susceptible 
to attack. Suitable room temperature anneals after polishing enable the 
traces of the grain boundaries on the surface to be readily etched. This 
must mean that during the ageing treatments the impurity content of the 
oxide film along the grain boundaries increases. Moreover, the only 
conceivable way in which the impurities could accumulate at the surface 
is by diffusion in the grain boundaries. 

Experiments on the self-diffusion of silver by Hoffman and Turnbull 
(1954) have shown that the rate of diffusion along the grain boundaries is 
appreciably greater than in the bulk of the crystals for temperatures well 
below the melting point. Achter and Smoluchowski (1951) also observed 
greater diffusivity along grain boundaries for the diffusion of silver into 
copper. Itseems reasonable to suppose that the impurities in the aluminium 
would, at room temperature, be able to diffuse more readily in the grain 
boundaries than elsewhere, and provided they collected in the oxide film, 
this would lead to selective attack along the grain boundaries. 

The fact that the duration of the room temperature anneal required to 
produce extensive grain boundary pitting in group A (recrystallized at 
room temperature) was much less than for B and C (annealed six days at 
600°C) is probably explained by the different treatments given to the 
groups before polishing. 

During the process of recrystallization impurities have ample oppor- 
tunities to occupy many of the available sites on the grain boundaries. 
An impurity atom on the grain boundary will remain there, provided its 
adsorption energy to the boundary is great enough to prevent it diffusing 
into the bulk of the metal. This is more likely to be true at room tempera- 
ture than at 600°C. Consequently it is possible that aluminium recrystal - 
lized at room temperature (group A) would possess a higher concentration 
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of impurities on the grain boundaries than aluminium which had been 
annealed at 600°c for six days (group B). If this is true then group A 
aluminium would require a shorter ageing time (to render the boundary 
traces on the surfaces attackable) than group B. 

In the case of the aluminium which was allowed to recrystallize at room 
temperature it might be expected that in addition to a general high 
concentration of impurities on the grain boundaries some portions of large 
angle boundaries might possess large enough accumulations of impurities 
to be attacked immediately after polishing, provided these regions are 
intersected by the plane of polish. A small amount of this was in fact 
observed and, moreover, such an occurrence on the surfaces of group B 
was very rare. 

It might, at first, be expected that the material comprising group C 
(99-99%) pure would exhibit marked grain boundary attack after shorter 
periods of room temperature ageing than in the case of group B, because of 
its greater impurity content. No marked evidence of this was observed, 
the behaviour of the two groups being more or less the same. 

A possible explanation is that the impurities in the material, as received, 
were in such stable precipitates that very little change occurred during 
annealing at 600°o, with the result that only small quantities were available 
for diffusion at room temperature. Even if the impurities were all in 
solution at 600°c they would tend to form precipitates on cooling ; leaving 
only a small amount still dissolved (in addition to small precipitates 
which may possibly redissolve) to take part in the diffusion processes we 
are considering. 

If such a situation is true, it may well happen that the times of room 
temperature anneal required to stimulate grain boundary etching would 
be approximately the same for groups B and C, as is observed in practice. 
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ABSTRACT 


The conditions for producing runaway electrons in a fully ionized gas are 
considered, using Chandrasekhar’s coefficient of dynamic friction. The 
runaway currents emitted continuously from linear plasmas are then esti- 
mated and are shown to be in agreement within an order of magnitude with 
the results from some preliminary experimental work. The rate of increase 
of runaway currents in closed circuits, such as toroidal plasmas, is estimated 
and is found tobe proportional to ¢? (¢{=time). From some recently published 
results with toroidal plasmas it is shown that the runaway current, provided it 
is contained and is not drastically affected by the magnetic fields, can rapidly 
become comparable in value with the conduction current. In the experi- 
mental work with linear plasmas, X-rays were observed having an energy 
greater than the applied potential difference to the apparatus. 


§ 1. [yTRODUCTION 


Ir is readily shown that in a fully ionized plasma there must always be some 
fraction of the electrons which will be decoupled by an electric field, no 
matter how small this field might be, and will then be continuously 
accelerated. Thus, in addition to the conduction current there will 
always be a ‘runaway’ current. Hoyle (1949) and Parker (1957) have 
invoked the idea of runaway particles in connection with the physics of 
solar flares; Gibson (1957) and Linhart (unpublished) have discussed the 
subject from the point of view of toroidal discharges. 

In the present paper the subject is treated briefly and the runaway 
currents emitted continuously from linear plasmas are estimated; 
similarly, runaway currents in closed circuits such as toroidal plasma 
are also estimated. As one might expect, the value of these currents is 
critically dependent on the strength of the electric field. 

Given conditions which are otherwise steady, a runaway current con- 
strained to flow in a closed circuit increases at a rate of approximately 2, 
where ¢ represents time. A runaway current in many cases can therefore 
become comparable in value with the conduction current and have an 
appreciable effect in the hydromagnetic behaviour of plasmas. 

The results of some preliminary experimental work with neutralized 
runaway electron beams, emitted endwise from linear plasmas, are shown 
to be within an order of magnitude agreement with the theory. 
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§ 2. Toe Runaway Errecr iw a Fuuiy Ionizmp Plasma 


Let an electron of mass m, and charge e,, having an initial velocity »,, 
move through a plasma containing particles of mass m,, charge e,, density 
n, and temperature 7’. Chandrasekhar (1942, 1943) has investigated a 
similar problem of a star moving through a star field having a Maxwell— 
Boltzmann velocity distribution. On making suitable substitutions in 
Chandrasekhar’s expressions one finds that the mean rate of increase 
Av in the direction of the initial velocity, due to encounters producing 
small deflections, is 


(Anr= — Se (14 St) (1-2) nah) 


Bay 2 
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in which dn; is the number of the ith particles having speeds v,>v,, and 
D;,=47e,"e? In A, where A is the mean ratio of the Debye shielding distance 
to the impact parameter for large angle scattering (Spitzer 1956). In a 
plasma consisting of electrons (‘=1) and positive ions (‘=2), with 
N,~N =n, the coefficient of dynamic friction given by eqn. (1) becomes 


(Av) = — 


sa a (20 — 200, 0n,) 9, 2 = (2) 
Poe 

When an electric field #, having an angle @ to the velocity v,, produces a 
rate of change of momentum opposite to and greater than the dynamic 
friction, the mean electron velocity increases. Thus, provided the 6n, 
are small, there is a critical speed v, given by He, cos@+m,(Av,) =0, 
such that all electrons having v >v, and a velocity component in the same 
direction as # gain energy continuously and contribute tc a runaway flux. 
To a sufficient approximation it will be assumed that the fraction of 
electrons qualifying for the runaway state is 6n,/2n, and the critical speed is 
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always provided it exceeds the most probable random speed of the 
electrons (and therefore also of the ions) and thus allows the non-dominant 
terms of eqn. (2) to be neglected. Therefore, the assumption made is 
that x,> 1, where 

eee (oD Zuerl a0 a 6 (4) 


= 0, 


and m,=m, e,=e, T,=T'; or, alternatively, the present treatment applies 
only for an electric field satisfying # <H,,, in which 


He 60105? aff volts ems?) ys.) a= (5) 
and In A~13. 

So far only electrons have been considered. However, a necessary 
condition for a runaway ion is that it has initially an energy at least as 
great as the most probable electron speed, and therefore for most purposes 
such ions can be neglected, 
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§ 3. EsTrIMATION OF RUNAWAY CURRENTS 


3.1. General Considerations 
On integrating the Boltzmann equation for the runaway electrons, one 
obtains 


ae +Sdivj,= | Ah dtu! Wisse eee 
where 7, is the number density, j, the current density, and Af,(v) is the 
rate at which electrons qualify for the runaway state per unit volume of 
phase space at velocity v. In this treatment the effect of a magnetic field is 
not considered; the presence of a magnetic field parallel to H, however, 
should not appreciably affect the result. When x,> 1 it can be assumed as 
a first approximation that for a zero electric field the integral expression on 
the right-hand side of eqn. (6) is equal to one half the rate at which electrons 
are scattered out of a v>v, spherical region of velocity space. Hence, 


[afc dy =4 I. Anv? f(v)r— dv 


where f(v) is the Maxwell—Boltzmann distribution function and 
7~v/{Av,) ~v?m?/3Dn. Thus, eqn. (6) is now very approximately 


7] 
St +o div j,= Y(«,) MR ce, 
; 3n?2D exp (—2%,”) 
Y'(a,) & ————_—_ 2 
(e) ere 
and numerically ‘Y'(x,) = 107n? exp (—x,?)/T3?z,?. 


3.2. Continuous Runaway Currents in Linear Plasmas 


For a uniform electric field in a steady-state plasma of length z, 

én,/0t=0 and oe (7) can be integrated to give 

see! (x,\ie 

=1-7'x 10*'n*2 exp (— 2,7)/T?"22 amp emis. ee eee oe 
For example, when n=10'%cm-8, 7’'=105°o, H=2-voltcm-, then 
«,?=H,,/H=3, and the runaway current amounts to 0-9 amp cm-? per 
unit length of plasma. The extreme limit for the validity of the present 
treatment corresponds to x,= 1, or H=6 volts cm~ in the present example, 
giving a current density of 20 amp cm per em of plasma. For com- 
parison, the random current density amounts to 80 amp cm-2. For an 
applied electric field greater than H,,, as given by eqn. (5), the velocity 
distribution will tend to become appreciably asymmetric and it is therefore 
unlikely in a steady state condition that a substantially greater current 
will be achieved. Owing to the rapid way in which ¥(a,) varies, appreciable 
runaway current in linear laboratory discharges is only realized when the 
electric field is at least of the same order of magnitude as ee 
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Tt V=Ez is the electric potential difference across the plasma column 
an alternative expression for eqn. (8) is 


y 


dt 
Jr=In Tp OXP (— He") eer, me eee ae) 


for a random current density of Jpr=ne(kT /27m)'2/c, 


3.3. Runaway Currents in Toroidal Plasmas 


Consider now a fully ionized uniform plasma in which the runaway 
current can form a closed circuit, as for example, a toroidal shaped plasma. 
Let this current be contained in the plasma for a time ¢ seconds. For a 
uniform plasma div j,=0, and eqn. (7) may be integrated to give a runaway 
electron density of 

n,=t¥(x,) Peep. oie es ier LO} 


provided that in all other respects the plasma is in a steady state condition. 
A steady state condition implies either that the random electrons are being 
replenished, or that the loss of random electrons is small in the time ¢. 
If, however, the initial random electron number density n is not a constant, 
then 

SE RI CORN East LORY ST Pooh ae a Sa) 


In a similar discharge (Thonemann et al. 1958) with n~10!4em-, 
T ~5 x 10° °c one finds n,/n ~ 8 sec from eqn. (10) when the electric field 
is Lvoltem™. For an electric field of 2 volt cm, n,/n increases initially 
at arate of 6 x 10?sect; the runaway flux if it is totally contained achieves 
a value of 1/e of the initial random electron density in a time of the order 
100 usec. 

Since 0j,/0t=n,e?H/mc, one finds from eqn. (10) (again assuming that 
other conditions are steady) that the circulating runaway current is 


9, =e ¥(x,)/2me 
=1-5x 10-(nt)? Hexp(—2,2)/T?"a2ampem . . (12) 
for a containment time ¢. Or, using eqn. (11) 
j= ne Eft — nin (1+t¥(a,)/n)/P(x,)}/me. pte ee) 


On inserting the previous values for n and 7’ one finds that in an electric 
field of E =0-8 volt cm-! the runaway current density is 7 x 10°¢?amp cm? 
and in 100 psec achieves a value approaching the conduction current density 
which is of the order 100ampem~. It must be emphasized that this 
method of estimating runaway fluxes has neglected the effect of magnetic 
fields, and the true results may therefore be considerably modified for 
example by a transverse field of the kind occurring in a pinched discharge. 


+ In the present context these values of n and 7’ are used purely as an example 
and the precise values are not therefore important. 
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It is interesting to notice that in an idealized closed runaway circuit the 
current density increases as /? and, unlike the linear plasma case, is inde- 
pendent of the linear dimensions of the circuit. In astrophysics, where 
the time scale is large, it would appear that runaway currents can become 
important if not predominant even when 4,> 1. 


§ 4, DESCRIPTION OF APPARATUS 


The apparatus that was used in an attempt to produce and study low 
and high energy runaway electron beams is substantially as shown diagram - 
matically in fig. 1, and in some respects is described elsewhere (Harrison 
and Dawton 1958). Section (a) consists of a water-cooled low voltage 
are operating continuously at a gas pressure of 50-100 wHg. A plasma 
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Schematic diagram of apparatus for producing runaway electron beams. 


beam is drawn through an aperture of radius r, into the next section (bd), 
by an applied electric field, and is prevented from diverging by an axially 
directed magnetic field of approximately 2000 oersteds. The magnetic 
field is continued into the succeeding sections, although in (d) and (e) it 
is reduced in value to around 300 oersteds. To some extent the pressures 
in (a) and (6) can be independently controlled. A neutralized beam 
consisting largely of runaway electrons enters (c), where the apparatus 
is pumped, and is guided by the magnetic field through (d) and (e) and 
finally strikes a water-cooled calorimeter target (Harrison 1957). The 


Saas returns along six parallel conductors equally spaced round the 
eam. 
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§ 5. Low Enerey Runaway Execrron Bram 


For continuous low energy beams of a few hundred ev the first aperture 
had a radius of r,=0-25-0-50 mm in a tungsten dise. Argon was used in 
all the experiments; attempts to produce low energy neutralized beams 
with hydrogen were unsuccessful. Figure 2, Pl. 97 shows a runaway 
beam of electrons produced from an argon discharge. 

With applied voltages varying from 200 to 450 volts, beam currents 
between 0-1 to 1-0 amperes were achieved without any great difficulty. 
Ion bombardment caused considerable erosion of the tungsten aperture 
dise, and it was therefore necessary to change this electrode frequently. 
A typical result was as follows ; a beam of 0-7 amperes, having an apparent 
radius of 1-1-5mm was produced with 320 volts. The pressures in 
sections (a), (b) and (c) were 70, 1 and 0-2 wHg; respectively. From the 
temperature of the calorimeter-target the mean energy of the electrons 
was found to be 173 ev. On biasing the target over a voltage range of 
+30 to —50 with respect to the aperture-dise electrode immediately 
preceding it, the temperature was found to vary linearly and it was 
deduced that the beam contained relatively few electrons with an energy 
less than 50 ey. 

The beams contained very considerable ‘hash’, and occasionally 
fluctuations similar to relaxation oscillations were observed (Cutler 1956). 
An examination of the hash revealed a low frequency, large amplitude 
component in the neighbourhood of 1 Mc/s which agreed moderately 
well with the expected frequency of the ion plasma oscillations (Pierce 
1948). 

From probe measurements made with plasma beams in weak magnetic 
fields, it was estimated that within a factor of 2 the electron density and 
temperature were n~ 10% cm? and 7'~10°°c. For a plasma column of 
length z= 20 cm the maximum beam current, neglecting the effect of the 
magnetic field, is 3 amperes from eqn. (8), which agrees within an order 
of magnitude with the experimental results. The minimum electric 
field necessary for this current, by eqn. (5), is 6 volts cm~ giving a total 
potential drop across the plasma column in section (b) of at least 120 volts. 
No doubt there is some residual electric field in the sections (c), (d) and (e), 
but it is unlikely that these regions will contribute appreciably to the 
runaway flux owing to the low random electron current density. 
This point of view receives some confirmation from the fact that when 
the calorimeter target was biased the beam power varied propor- 
tionately. 


§ 6. HicH Enercy Runaway ELECTRON BEAMS 


To produce higher energy beams the apparatus was connected to a 10 uF 
condenser and pulsed at a rate of 1 sec! at voltages up to 10 kv. Across 
the low voltage arc of section (a), which was kept continuously running, 
was connected a 100 pF condenser. The radius of the first aperture was 
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increased to 1-2mm. ‘The results with hydrogen were substantially the 
same as those with argon except that slightly higher pressures were 
necessary. With argon the approximate pressures in sections (a), (b) and 
(c) were 50, 5 and 0-5 pHg, respectively. 

So far the experimental work is incomplete and only preliminary results 
can be given. It was found that the performance was critically dependent 
on the pressure in (b). When the pressure was too high the discharge in (0) 
of about 104 amperes was accompanied by a loud noise; as the pressure 
was reduced the discharge became silent. Infrequently, these ‘silent 
discharges’ were accompanied by the emission of sparks from the target. 
The discharge, including the beam current reaching the target, usually 
lasted for approximately 1 msec. As with the low energy beams, con- 
siderable hash was observed. A large amplitude component of around 
10 Mc/s was readily identified. For optimum conditions the current at 
the target, measured with a current transformer, was in the region of 
200 amperes. 

The target current was disappointingly small. The results, however 
lend support to the idea that substantially larger currents cannot be 
achieved by using electric fields of H>H,,. From the increased pressure 
in (b) one would then expect the current to be of the order of magnitude 
given when some allowance is made for the increase in the electron 
temperature. The mean temperature of the calorimeter-target, together 
with the current measurements, confirmed that the mean kinetic 
energy of the beam in electron volts was close to that given by the supply 
voltage. 

Measurements of the x-ray emission from the target, using a pro- 
portional counter and aluminium filters, revealed the interesting fact 
that the x-rays had a mean energy of 25-35 kev, although the beam energy 
had a maximum value of 10 kev. One possibility is that some of the 
electrons are accelerated owing to variations in the beam current. The 
maximum possible energy gained by this method is readily estimated. 
Thus, if eV, is the initial energy of a beam of intensity 4/,(1+ cos wf) 
and L is the inductance of the system, the maximum energy possible is 
given by 

V=Vy+3Lwly. 


For L~2x 10~¢ H and w/27~ 10’, this gives a value of 22 kv. Since the 
beam hash contained higher frequencies it would appear that the x-ray 
emission can be explained by an acceleration mechanism of this kind. 
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On the Temperature Dependence of Yield Stress for Aluminium 


By W. D. Syztwestrowicz and D. F. Grspons 
Bell Telephone Laboratories, Inc., Murray Hill, New Jersey 


[Received July 31, 1958] 


TuE temperature dependence of yield stress for aluminium has been 
investigated by Cottrell and Stokes (1955) for single crystals and for poly- 
crystalline material by Sylwestrowicz (1958). In both cases, the yield 
stress shows a strong dependence on temperature. A question has been 
raised, would the ratio of the yield stress to elastic modulus, o/G, also 
depend on temperature? The answer to this question is of significant 
importance for the theoretical consideration of the mechanism of strain 
hardening. In the majority of the proposed theories, the only factor 
dependent on temperature is the elastic modulus and this dependence is 
linear. Therefore, it should be expected that the ratio o/G should be 
independent of temperature. Using the values of Young’s modulus of 


Temperature (1 Temperature 


: G-1 

in °K in °K 

293°3 4-071 173-2 3°817 
273°3 4-024 153-2 3-782 
253-2 3:979 143-2 3°764 
233-2 3-937 133-2 3:747 
213-2 3:°894 717-2 3:663 
193-2 3°853 


Koster (1948), in the range from liquid nitrogen to room temperature, this 
ratio, o/H, is still strongly temperature dependent for polycrystalline 
material (Sylwestrowicz 1958). For single crystals, the temperature 
dependence of o/E is also strong at low temperatures but decreases near 
room temperature, Cottrell and Stokes (1955). The next question arises, 
what would be the result if the values of macroscopic elastic moduli are 
replaced by the elastic moduli measured in a single crystal in the plane and 
direction of the primary slip ? Since all theories of deformation are based 
on phenomena occurring on atomic dimensions, it is appropriate to use 
the value of the elastic modulus in the plane of slip. The shear modulus @ 


Correspondence 1327 


in the slip plane, {111}, can be expressed in terms of the three elastic 
moduli for crystals of cubic symmetry and is 


G-t=4 [84,4 4(8,,—Sj,)]. 


The isothermal elastic moduli, 8,,, S,., 84, have been measured as func- 
tions of temperature down to 60°K (Sutton 1953) and the values for @ 
have been calculated as a function of temperature according to the above 
expression and are given in the table. 

Although these values apply to the primary slip plane, the temperature 
dependence for the secondary slip planes would not be expected to differ 
appreciably. For room temperature and that of liquid nitrogen, the ratio 
Gy /G,,y. ‘is equal to 1-11 and is again smaller than would be needed 
to make the ratio o/G independent of temperature, as the ratio of 
Cp p/o,,.x. for aluminium is about 0-74. 
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The Evolution of Extragalatic Nebulae and the Origin of 
Metagalactic Radio Emission} 


By G. R. BuRBIDGE 


Kellogg Radiation Laboratory, California Institute of Technology 
and Yerkes Observatory, University of Chicago 


[Received June 23, 1958] 


Iy a paper bearing this title, Bruce (1958) has attempted to argue that the 
jet in the elliptical galaxy NGC 4486 (Messier 87) is due to an electric 
discharge mechanism. His last reference to the observations is to those of 
Baade and Minkowski (1954). However, it should be pointed out that in 
1955 Shklovsky (1955) and others predicted that the jet in M 87 was due 
to synchrotron radiation in the optical frequency range, in which case it 
was possible that it would have a high degree of polarization. Following 
this Baade (1956) with the 200-inch telescope at Palomar detected 
polarization in the jet amounting in parts to some 40%. Then the author 
(Burbidge 1956) was able to calculate that the total energy of high energy 
particles and magnetic field necessary to account for the jet in this way was 
NE 

+ Supported in part by the joint programme of the Office of Naval Research 
and the U.S. Atomic Energy Commission. 
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at least of the order of 10° ergs. This theory is the only one able to 
explain easily the high degree of polarization. Itis also generally accepted 
at the present time that the synchrotron mechanism is the most plausible 
explanation of the non-thermal radio sources, of which M 87 is one, both 


inside and outside our Galaxy. 
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The Evolution of Extragalactic Nebulae and the Origin of 
Metagalactic Radio Emission 


By C. E. R. Bruce 


Electrical Research Association Laboratory 
[Received July 8, 1958] 


In comment on Dr. Burbidge’s letter (Burbidge 1958) I should like to 
emphasize that my sole purpose in writing my original letter was to point 
out that the discovery of the ‘jet’ in NGC 4486 and of its being a radio 
source, fulfilled a prediction made in 1944 that such ‘jets’ should exist 
in nebulae of the type of NGC 4486. I was certainly well aware of 
Shklovsky’s suggestions and of earlier rather similar work of Alfven 
(Alfven and Herlofson 1950). However, Dr. Burbidge will agree that the 
great difficulty confronting all such theories, is, as Cowling (1957) and 
others have emphasized, to account for the magnetic field. As I have 
pointed out in the E.R.A. Report mentioned in my earlier letter 
(Bruce 1958), this difficulty would appear to be obviated by the electrical 
discharge theory. 

I should also like to refer to still later work of Woltjer (1958) in Leiden, 
which has appeared since my E.R.A. Report was written. Shklovsky and 
many others have drawn attention to the close similarity between the 
radiation emitted by NGC 4486 and that emitted by the Crab Nebula. 
In the latter, however, a much more detailed study of the polarization of 
its light is possible. As a result of such a study Woltjer has shown that 
the varying directions of the polarization can be accounted for if electric 
currents flow along the filaments. So that another fourteen year old 
prediction of mine would seem to be fulfilled. 

Furthermore in recent communications (Bruce 1957) I have shown that 
the magnetic pinch effect in these electrical discharge channels will give 
rise to relatively high pressures therein and to jets of gas along their 
length. In the long-period variables and combination-spectra stars the 
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observed velocities of the hot gas in which the emission lines originate 
agree with those calculated on this view from the observed temperature 
of the gas in these stars. The conditions in these two types of stars cover a 
range of 10 to I in velocity, and 100 to 1 in temperature. In accordance 
with this theory, Woltjer has observed that the filaments of the Crab 
Nebula are also the sites of jets of gas reaching around 1000 km per second. 
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Measurements of Unit Cell and Physical Dimension Changes 
of Molybdenum after Neutron Irradiation 


By J. Apam and D. G. Martin 
Metallurgy Division, Atomic Energy Research Establishment, Harwell 


[Received July 22, 1958] 


Unit cell and physical dimension changes of molybdenum have been 
measured after irradiation at 30°c and also after subsequent annealing. 
The physical dimension measurements were made on randomly oriented 
polycrystalline bars + in. x ~in.x lin. which had been annealed at 
750°c ; the end faces were ground and lapped parallel to one microinch. 
Changes in sample length were measured to an accuracy of 2 parts per 
million in a temperature controlled room on a calibrated mechanical 
comparator (magnification 4500), using unirradiated samples as standards. 
Specimens were supplied by Johnson Matthey Ltd. and Murex Ltd. 
Both batches had 0-01°% of iron as principal impurity and showed no 
significant difference in behaviour on irradiation and subsequent annealing. 
Unit cell measurements were made on photographs taken in a precision 
X-ray powder camera (Adam 1954) temperature controlled at 26-D°G. 
Errors due to non-uniformity of film shrinkage were eliminated by 
printing an accurate scale on each film before processing. The final 
value of the unit cell size, which was accurate to about 15 parts per million, 
was derived using the standard Nelson and Riley (1945) extrapolation. 
For the powder used (Hopkin and Williams Ltd.}) a=3-14725 + 0-00005 A 
assuming the wavelength of copper Ka, radiation equal to 1:54050 A. 


. i. 
+ Irradiation experiments on bars and powders were started independently 
and materials from different sources were used. No impurities were detectable 
spectroscopically in the powder. 
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- All samples were irradiated in vacuo inside a hollow uranium cylinder in 
the Harwell reactor BEpo. The irradiation position and the method of 
neutron dose measurements were the same as Kinchin and Thompson’s 
(1958). In this position the fast and slow neutron fluxes are approxi- 
mately equal. ; 

Figure 1 illustrates the fractional expansion in unit cell and physical 
dimensions as functions of the measured slow neutron dose. Hach 
irradiated sample was annealed for one hour at successively higher 
temperatures and its dimensions measured. The fractional decrease in 
length obtained from a typical specimen, showing two principal annealing 
stages at 180°c and 550°c is given in fig. 2. 40% of the original elongation 
anneals during the first stage and is independent of neutron dose. The 
fractional decreases in the unit cell expansion measured on a powder 
annealed for one hour at 200°c and 700°C are approximately the same. 
Unfortunately the x-ray measurements are an order of magnitude less 
accurate than those of physical dimensions and therefore the shape 
of the annealing curve cannot be obtained from them. 

An annealing peak at 180°c has been observed by Kinchin and 
Thompson (1958) who measured the electrical resistivity and stored 
energy after neutron irradiation, and also during the recovery of the cold 
worked metal (Martin 1957). 

It is possible to estimate the concentration of defects from the observed 
strain. Tewordt (1958) has shown that the changes in unit cell size of 
copper for 1° of interstitials and vacancies are +0-6 and —0-2% 
respectively. Since a similar value is obtained for vacancies by combining 
Girifalco and Streetman’s (1958) calculations of the relaxation of a vacancy 
in molybdenum with elasticity theory (Eshelby 1954), it will be assumed 
that Tewordt’s results apply also to molybdenum. The concentration of 
defects derived from theoretical considerations of radiation damage 
(Kinchin and Pease 1955) is of the order of a hundred larger than that 
estimated from the observed expansion after irradiation ; this suggests 
that at least one fundamental defect is mobile at 30°c. Although caution 
must be exercised in comparing quantitatively the unit cell and physical 
dimension changes because of the difference in purity of samples, it can 
be accepted that the increase in physical dimensions is larger. Since any 
uniform distribution of defects within the lattice would affect cell 
dimension and macroscopic changes similarly, to account for the larger 
macroscopic expansion we must postulate that new lattice sites were 
formed during irradiation. This could occur if single interstitial atoms 
possessed a lower activation energy of migration than vacancies. The 
defects remaining in the lattice after irradiation at 30°c producing unit 
cell expansion are thus either interstitials trapped by impurity atoms or 
clusters of interstitials. The depths of the traps or the sizes of clusters 
are probably not identical because of the wide range of temperature over 
which annealing takes place. The annealing peak at 180°C can be 
attributed to a particular trap depth or cluster size distribution or 
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Fig. 1 
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alternatively the annealing process may be accelerated by vacancy 
movements if these become appreciable at about 180°c. 
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The Properties of Beryllium-11 


By D. E. AtBurGcER and D. H. WILKINSON{ 
Brookhaven National Laboratory, Upton, New York 


[Received August 21, 1958] 


AJZENBERG AND LAURITSEN (1955) suggested that the energy difference 
between the ground states of “Be and ™B should be about 11-5 Mev and 
that “Be should therefore be stable by about 0-5Mev against neutron 
emission. In this case the Q-value for the reaction |'B(n, p)''Be would be 
about —10-7mev. Nurmia and Fink (1958) discovered a beta-activity of 
half-life 14-1 + 0-3 seconds induced in boron by the 14Mev DT neutrons 
and ascribed it to4Be. We have studied this body, made in the same way. 
We have confirmed its assignment to Be and have determined certain 
properties of the decay. The work was carried out, using Nal(Tl) and 
plastic scintillation detectors both singly and in coincidence and the 
measurements are presented briefly in the table. 

The energies of the "B states in the first column are taken from the 
literature (Ajzenberg and Lauritsen 1955). The entry “ % betas’ gives 
the percentage of all 1'Be beta-decays that take place to the "B state in 
question. “ % gammas ” gives the abundance of the various gamma-rays 
per “Be decay. The half-life entered for the ground state was measured. 
on beta-particles of energy greater than 3-0mev. That for the 2-14 Mev 
state was measured on the gamma-rays from that state. The third 
half-life quoted was measured on all gamma-rays of greater than 3-0 Mev. 
The closeness of these half-lives and the agreement between our measured 
gamma-ray energies (fifth column) and the known energies of the various 
states of '\B confirm the assignment of the activity to “Be. We combine 
these measurements to give an average half-life of 13-57+0-l5sec. Our 


{ Work performed under the auspices of the U.S. Atomic Energy Commission, 


{ Permanent address : Clarendon Laboratory, Oxford University, Oxford. 
England. , s 
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various measured beta -particle transition energies (second column) may be 
combined with the known energies of the “B_ states to yield : 
"Be—'B=11-48+0-15Mev which corresponds to a mass defect for 
NBe: M—A=23-39 + 0-15 Mev. 

The log ft values of the various beta-ray decay branches are rather 
large but those leading to the 6-76 and 7-99 Mev states probably indicate 
allowed transitions. The parity of the former state is probably odd (Ccx 
and Williamson 1957) which gives odd parity for “Be if indeed the transi- 
tion leads to that state but of this we cannot be sure from our data. The 
beta-decay to the 3— state at 4:46 Mev is strongly forbidden while transi- 
tions to the }— ground state and }— state at 2-14 Mev are considerably 
faster although still surprisingly slow for allowed transitions. This suggests 
J=%}-— for Be which is as expected from its position in the 1p-shell. 


3 state 0 -f 
E OF cians vi Half-life 
MeV ) piten) mp notes log ft pion) gammas (sec) 
11-48 + 0-15 61 6-77 — -— 13-68 + 0-15 
138 9-32 + 0-15 29 6-63 2-121 + 0-010 oe 13-47 + 0-20 
459 = =—0-2 25°2 — —) 
034 = <0-2 338-2 ce a | 
6- 5: 


- 13-48 + 0-20 


Although the speed of the two beta transitions of highest energy may 
suggest first-forbidden transitions (even parity for ''Be) the strong relative 
discouragement of the transition to the 4-46 Mev state argues against this. 
It seems most likely that Be is J=}— and that we have two rather 
retarded allowed transitions whose theoretical analysis should be very 
interesting in view of the fact that the independent-particle model wave 
functions are fairly well known in this region (Kurath 1956). 

We have considered the possibility that the 13-6sec period belongs to 
an isomeric state in “Be but we have found no evidence for this in a 
search for low-energy gamma-rays accompanying the decay. 

We would like to thank Mr. Stuart Crampton for his help. 
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REVIEWS OF BOOKS 


Biophysical Chemistry. Volume 1. ° Thermodynamics, Electrostatics and the 
biological significance of the properties of matter.’ By J. T. Epsaun and 
J. Wyman. (Academic Press 1958.) £5. 


A TexTBooK such as the present one must first be judged by the quality of its 
treatment of certain of the basic principles of physical chemistry. Edsall and 
Wyman are to be congratulated on producing a book which could well serve as 
a standard text for advanced undergraduates on thermodynamics and the 
properties of ionic solutions. The greatest care has been taken in presenting 
the physical ideas which motivate the choice of mathematical methods. Only 
in the use of Caratheodory’s rather than a more elementary treatment of 
the second law of thermodynamics is the reader’s understanding of abstract 
mathematics tested more than is absolutely necessary. 

By the clever choice of examples the authors have succeeded in introducing 
so much material dealing with the properties amino acids, peptides, poly- 
phosphates, etc., that the sections on the physical chemistry of electrolytes 
also form a detailed introduction to this aspect of biophysics. This material 
taken together with the chapters on the constitution and structure of biological 
materials make up an excellent introduction to the more biological parts of the 
subject. 

This book can be wholeheartedly recommended both to students and to 
biologists who wish to brush up their physical chemistry. LEO: 


Electric Conduction in Semiconductors and Metals. By W. EHRENBERG. 


(Oxford: Clarendon Press). [Pp.389.] 63s. 


TuIs is mainly about semiconductors. The first half derives the electrical 
transport properties of solids, first in the mobility approximation and then 
more rigorously from the Boltzmann equation. The properties of actual 
semiconductors are then summarized. The final chapters give the theory of 
rectification and transistor action. Although there is a good deal of mathematics 
it is all relatively elementary, and should easily be followed by physicists 
working on semiconductor devices. Perhaps Chapter 7 is a little too heavy; 
the brute-force solution of the Bloch integral equation is extremely complicated 
to explain, and has little direct physical significance. Also, it is surprising 
that the concept of ‘deformation potential’ is not mentioned. The book 
as a whole will be useful for the careful detailed derivation from first principles 
of a number of standard formulae, and for the discussion and clarification of 
various difficult points which are not adequately treated elsewhere. J. M. Z. 


A History of Mathematics. By J. F. Scorr. (London: Taylor & Francis, 
Ltd.) [Pp. x+266, with six plates.] 3 guineas. 

Dr. Scort now follows his specialist works on seventeenth century mathematics 

with a general survey of the history of mathematics from its beginnings to 

the early nineteenth century. 

_ For such a task the first requisite of the author is a sound knowledge of the 
literature, particularly as it affects the periods with which he is least familiar. 
Unfortunately, this knowledge Dr. Scott does not always possess. His 
bibliography includes the great classics of Montucla and Cantor, but it omits 
most of the important modern works on the history of mathematics. In some 
fields these modern works have compelled a revision of earlier views, and 
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accounts based on earlier books are no longer acceptable, This is particularly 
the case with Babylonian mathematics, a field in which exciting developments 
have taken place in recent decades. Dr. Scott takes the surprising course of 
quoting the work of the man chiefly responsible as reported in a textbook for 
undergraduate mathematicians, and he pays the penalty, 

On Dr. Scott’s home ground of the seventeenth century, and on topics where 
recent research has added little of importance, this book is of value. Elsewhere 
it must be used with caution. M.A. H. 


Dislocations and Mechanical Properties of Crystals. By J. C. Fisoer, W. G. 
JoHNsTON, R. THomson and T. VrEELAND Jr. (New York: John Wiley & 
Sons, Inc.; London: Chapman & Hall.) [Pp. xiv+634.] 120s. 


TuHis book describes the proceedings of an International Conference on 
dislocations held at Lake Placid in September 1956. The Conference was 
organized by the U.S. Air Force Office of Scientific Research, the Air Research 
and Development Command and the General Electric Research Laboratory, 
and covered eight main fields. 

it is not possible in a short review to do more than indicate a few general 
impressions. One notes, on one hand, the elegant and beautiful experimental 
techniques that have been developed in recent years to observe directly the 
presence of dislocations; on the other, the underlying theoretical uncertainty 
as to whether these techniques do, in fact, always reveal the most relevant 
dislocations. Secondly one is impressed by new theoretical mechanisms that 
continue to be elaborated in an attempt to explain more adequately such 
problems as work-hardening, yield-point phenomena and the effect of reversed 
stresses. Thirdly one is struck by those areas in which dislocation theory 
appears to provide a remarkably successful synoptic explanation of a whole 
class of experimental phenomena, such as Liicke and Granato’s work on 
internal friction. But chiefly there remains the impression of a Conference 
in which the participants discussed freely and openly their own and each 
other’s ideas. If for no other reason the editors are to be congratulated on 
conveying so effectively an atmosphere of vigorous and thought-provoking 
discussions. As a result they have produced an account of recent advanced 
experimental and theoretical studies in the field of dislocations that is both 
authoritative and stimulating. 10, A, 


Mass Spectroscopy. By Hmnry E. DuckwortH. (Cambridge: University 
Press, 1958.) [Pp.ix+206.] 35s. 


In the past decade mass-spectrometers have been incorporated in many kinds 
of apparatus used to study phenomena which depend on atomic or molecular 
mass. In nuclear physics, for example, the masses of stable and radioactive 
isotopes, neutron absorption cross sections, fission yields and nuclear moments 
have been studied. Chemical and biological processes are more easily 
interpreted when isotopic abundances are measured, and these also play an 
important réle in geological dating. The possibility of mass analysis of reaction 
products has led to great advances in molecular physics and the understanding 
of dissociation processes and chemical reactions. 

This monograph contains a concise description of two aspects of mass- 
spectrometry, firstly, the design and operation of different types of instrument, 
and secondly, the application to various fields of research and measurement. 
The first six chapters are concerned with the design principles, the ways of 
obtaining sufficient resolution and intensity, and the techniques used for the 
production and detection of the ion beams. No attempt is made to cover the 
engineering, electronic and vacuum technology. The remaining five chapters 
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are devoted to the determination of atomic masses and isotopic abundances, 
applications in nuclear physics, the study of ionization and dissociation of 
molecules under electron bombardment and finally, applications to geology. 
Tables of isotopic abundances are included, and there is a generous bibliography. 
The book gives a very good indication of the sort of research which is possible 
when a mass-spectrometer is included in the apparatus and typical modern 
instruments are described in ample detail. K.S. 


Atomic Physics and Human Knowledge. By Nrmts Boner. (New York: John 

Wiley & Sons, Inc.) [Pp.101.] $3.95. 

Tuts collection of articles, extending over the years 1932-57, is a ‘ theme 
with variations ’ on the principle of complementarity. Although only one of 
the essays attempts anything like a systematic treatment (Bohr’s discussions 
with Einstein) their main interest lies perhaps in the evidence they yield of the 
repeated attempts of the author to arrive at a clearer statement of the com- 
plementary idea (cf. pp. 6, 19, 25, 39, 60, 72, 90,99). This is said to “ sym- 
bolise the fundamental limitation, met with in atomic physics, of the objective 
existence of phenomena independent of the means of their observation ”’ (p. 7). 
Elsewhere it is a “ consistent generalisation of the idea of causality ”’ (p. 27); 
or again, it ‘‘ characterises ” the “‘ mutual exclusion ” of certain experimental 
procedures designed, e.g. to measure spatio-temporal and momentum or energy 
aspects, (p. 61). The principle, then, draws our attention to certain incom- 
patibilities. Does it do more than this? By 1955 a still sharper definition is 
adopted: the mutual exclusiveness of certain experimental situations stems 
from the “ very definition ” of the clocks and measuring rods adopted (p. 90). 
However, even here the complementary principle only tells us that this 
exclusiveness “‘ must be termed complementary ....’’ (/bid., my italics). 

Bohr also attempts to apply the complementarity idea to the spheres of 
biology, psychology and cultural history. Here he seems less sure, and it is 
indeed questionable whether analogies suggested by the term complementarity 
are relevant to the extension of the concept beyond the ground of quantum 
mechanics. Thus, in the free-will controversy, man as actor and as object are 
hardly complementary aspects in the quantum sense. And in the biological 
sphere, the theory that complementarity applies to the organism as a living 
and as a dead object of study is based on the weak analogy of some sort of 
limitation of investigation as a feature shared with the physical situation. 

GERD BUCHDAHL. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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Fig. 1 
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(a), (b): opposed faces of a ‘ blue-white * octahedral diamond; 
(c) irataneont region within this diamond, seen through face (a). 
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Fig. 2 


Central area of fig. 1 (2)—normal photomicrograph. 


Fig. 3 


Central area of fig. 1 (a)—‘ schlieren ’ illumination. 


Original magnification x 160; these two figures are reduced by one-half. 
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(a) Below upper face. 


(b) Originating from lower edge. (The reflections of tunnels can be seen in the 
adjacent octahedral face.) 

Tunnels (indicated by arrows) penetrating into a diamond etched in wet 
hydrogen at 820°c. x 20. 
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(¢) (d) 


(a)—(c) Diamond etched in chlorine with trace of oxygen at 800°C. (a) Pyra- 
midal and flat-bottomed pits. 1700. (b) and (c) Apparently non- 
crystallographic features. 1700. (d) A similar feature to (b) and (c) 
observed on an unetched octahedron. Electron micrograph. x 7500. 
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Pits in trigon orientation etched in blue ground at about 1450°c. (Surface 
Oo o 
layer of graphite also present.) x 50. 


Fig. 4 
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Pits in trigon orientation etched in blue ground at about 1450°c. x 420. 
Dark ground illumination. 
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Bottom of deep pits etched in blue ground. x280. Dark ground illumination. 
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A typical avalanche in InSb (111) faces <130: 
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Fig. 1 


Pattern from silver iodide. This pattern is common to both cubic—{211] 
azimuth—and hexagonal—{11.0] azimuth—silver iodide. 


Fig. 2 


Pattern from a mainly cubic deposit of oriented silver iodide, 
[110] and [110] azimuths. 
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Fig. 3 


Pattern from hexagonal ice—[11.0] azimuth—oriented on silver iodide. 
Glancing angle of incidence of electron beam, 1-6°. 


Fig. 5 


Pattern from hexagonal ice, [11.0] azimuth. The underlying silver iodide has 


been obscured by the ice deposit. Glancing angle of incidence, 0-¢ 
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Fig. 5 
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Fig. 7 


X 


x 50 


|. BRAUN, F. C. FRANK and G. MEYRICK Phils Magy ser<8, Volo37 Piezo: 
Fig. 9 
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Fig. 11 
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Photograph of a continuously running low energy beam. 


